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ABSTRACT 


A  system  of  n  first  order  linear  differential  equations  in 


matrix  form 


(*) 


dW 

dz 


=  z 


s-1 


k=o 


z  k  +  a  (z)  z 
s  ' 


-  s 


W  (s  =  0, 1,2, • • • ) 


where  h  is  an  integer,  A  (z)  =  A(z),  each  A  (k  =  0,1,2, ••• )  is  a 

O  K 

constant  n  x  n  matrix,  the  n2  elements  of  A(z)  (s  =  0, 1,2, • • • )  are 

s 

holomorphic  in  some  domain  £  extending  to  infinity  and  for  z  in 

lim  |(A  («))  .  .  |  =  K .  ,  <  oo  ,  has  formal  solutions  which  are  asymptotic 

z  — *•  oo 

expansions  (as  |z|  -*oo)  of  actual  solutions  in  certain  sectors  extending 
to  infinity  in  S. 


In  this  thesis  known  results  of  error  bounds  for  partial  sum 
approximations  of  formal  solutions  of  second  order  equations  to  actual 
solutions  of  second  order  equations  are  extended  to  the  general  system  (*). 
This  is  achieved  by  first  transforming  the  system  (*)  to  canonical  form. 
Formal  solutions  are  then  explicitly  constructed  for  the  transformed  system 
(henceforth  to  be  referred  to  as  (**))  and  error  bounds  are  obtained  for 
partial  sum  approximations  of  these  formal  solutions  to  actual  solutions 
of  (**). 


The  treatment  of  the  general  case  (*)  is  preceded  by  the  treatment 
of  the  special  cases  h  =  -2  (which  includes  h  <  -2),  h  =  -1,  and 
arbitrary  non-negative  h  with  the  eigenvalues  of  Aq  distinct.  The  case 
h  =  -2  is  ready  for  formal  solution  as  it  stands.  Transformations  are 
explicitly  constructed  which  transform  the  particular  h  >  0  case  when  the 


,n  " ■  i  .■ 


. 


eigenvalues  of  are  distinct  and  the  h  =  -1  case  to  canonical  form. 

A  known  existence  theorem  concerning  transformation  of  the  general  h  >  0 
case  is  stated  and  a  modification  of  this  theorem  corresponding  to  our 
modification  of  the  usual  h  =  -1  canonical  form  is  proved. 

Assuming  that  we  require  actual  solutions  of  (**)  for  which  the 
formal  solutions  are  asymptotic  expansions  (as  |z|  ->oo  ),  and  assuming 
the  domain  j£)  is  sufficiently  large  it  is  always  possible  in  the  special 
cases  h  =  -2  and  h  =  -1  to  bound  the  error  vector  by  use  of  a  single 
Volterra  vector  integral  equation.  A  simple  geometric  criterion  is  given 
by  which  one  can  determine  whether  a  single,  or  a  simultaneous  pair  of 
Volterra  vector  integral  equations  are  required  to  express  the  error  in 
the  case  when  h  >  0  and  the  eigenvalues  of  A^  are  distinct.  In  all 
these  special  cases  it  is  possible  to  choose  the  end-points  of  integration 
at  infinity.  This  is  no  longer  always  possible  in  the  general  case,  even 
if  is  a  neighborhood  of  infinity.  Although  it  may  be  necessary  to  choose 

one  finite  end-point  of  integration  it  is  nevertheless  always  possible  to 
express  the  error  vector  by  at  most  a  simultaneous  pair  of  Volterra  vector 
integral  equations. 

As  an  example,  the  general  n  =  2  case  is  completely  solved, 
including  a  flow  chart  illustration  of  the  transformation  of  this  case  to 
canonical  form. 

Real  variable  inequalities  useful  for  bounding  solutions  of  integral 


equations  are  obtained. 
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CHAPTER  I 


INTRODUCTION 


This  introductory  chapter  is  intended  to  facilitate  the 
understanding  of  the  thesis. 

Some  definitions  and  notations  are  established. 

A  brief  review  of  the  literature  leading  to  the  problem  of  this 
thesis  is  given. 

The  problem,  namely  that  of  obtaining  a  bound  on  the  difference 
between  actual  solutions  and  partial  sums  of  formal  solutions  is  described, 
and  the  need  for  solution  of  this  problem  is  emphasized. 

A  preview  of  the  organization  of  the  remainder  of  the  thesis 

is  given. 


1.1  Matrix  Notations  and  Norms 


We  shall  use  the  following  notations  to  refer  to  an  tn  x  n 


matrix  A; 


(1.1.1) 


A  = 


all 

a  o « o 

12 

ai 

In 

a  , 

a  o  o  o 

a 

21 

22 

2n 

a  , 

O  6  0 

3,  O  o  o 

a 

ml 

m2 

mn 

= 


To  refer  to  the  (isj)  "th  element  of  a  matrix  A,  we  shall  use 


(1.1.2)  (Aj.j  =  . 

If  the  matrix  A  is  square  (m  =  n)  we  designate  the  unit 
matrix  consisting  of  ones  along  the  diagonal  and  zeros  elsewhere  by  I. 
The  diagonal  part  of  the  square  matrix  A  with  diagonal  elements 

all’a229  ann  iS  desiSnated  bY  [an»a22>  °  *  ° »  annh  or  simply 


fc 


2  - 


by  D  when  it  is  convenient  not  to  illustrate  the  diagonal  elements. 

Furthermore  by  A  we  shall  mean  A  -  f-all’a22>  **’»  ann^  where  A  =  [a^]. 


We  shall  also  use  the  notation 


(1.1.3) 


'  [AU] 


i>  j  =  1,2. 


where  A„ . 

ij 


Ml  +  ^2  + 


is  a  x  pu  matrix  for  some  positive  integers 
•  +  \x&  =  n. 


u.  and 
i 


f 


In  carrying  out  error  analysis  we  shall  often  replace  an  m  x  n 
matrix  A  by  its  "comparison  matrix"  | A |  where 

(1.1.4)  | A i  =  [|{A)tJ|] 

i.e.  J  A  j  is  the  matrix  of  absolute  values  of  the  elements  of  A.  This 
applies  equally  to  the  case  when  m  is  an  arbitrary  positive  integer 
and  n  is  18  i.e.  to  the  case  when  A  is  a  vector. 


Given  two  ra  x  n  real  matrices  A  and  B  of  the  same 

dimension,  we  interpret  A  <  B  to  mean  (A) . .  <  (b) . .  for  each 

^  J  ^  J 

(i,j)(i  =  ls28 • ° » sm i  j  =  192s • • • ,n) „  We  note  that  this  definition  includes 
the  possibility  that  A  and  B  are  vectors. 


It  will  often  be  convenient  to  use  norms  of  vectors  and  matrices 


We  assume  A  =  [aij3  be  an  m  x  n  matrix.  Possible  matrix  norms  ||A|[ 


for  A  ares 


(a)  norm  ]|A||  =  max 

j 

(1.1.5)  (b)  i  norm  ||A||  =  max 

00  i 

m 

(c)  norm  ||A||  =  QT 

xiii 

We  again  note  that  the  definitions 
a  vector. 

Each  definition  of  ||A|| 
properties s 

(a)  ||A||  >  0  and  ||A||  = 

(h)  ||cA||  =  | c |  ||A||,  c 

(1.1.6) 

(c)  ||A  +  B||  <  || A ||  +  |[B ||  s 

(d)  ||AB||  <  ||A||  ||B||,  A 


3  - 


m 


(1.1.5)  include  the  case  of  A  being 


in  (1.1.5)  satisfies  the  following 


0  if  and  only  if  A  =  0 
is  a  complex  number 

A  and  B  both  m  x  n 
is  m  x  kj,  B  is  k  x  n. 


The  definitions  (a)  and  (b)  of  ||A||  in  (1.1. 5)  are  readily 
shown  to  satisfy  (1.1.6).  Case  (c)  of  (1.1. 5)  clearly  also  satisfies 
(a)8  (b)  and  (c)  of  (1.1.6).  With  B  =  [b^  ]5  the  proof  for  (d) 
of  (1.1.6)  is  given  in  [12],  page  228. 


Combining  the  definitions  of  (1.1.4)  and  (1.1.5),  it  follows 


that  (1.1.7)  |||A|||  =  ||A||. 


-  4  - 


1.2  Some  Transformations  on  a  2*2  Matrix 


A  general  theory  of  functions  of  matrices  is  presented  in  the 
first  volume  of  [5].  In  this  section  we  shall  consider  some  functions 
of  the  matrix 


(1.2.1) 


where  a*  J3,  7  and  6  are  non-negative  numbers  „  Suppose  this  matrix 
has  distinct  eigenvalues.  Then  there  exists  a  non-singular  matrix  T 


(1.2.2) 


such  that 


(1.2.5) 


where 


and  A 
2 


are  the  eigenvalues  of  A. 


Substituting  (1.2.2)  into  (1.2.3)  an<*  equating  off-diagonal 
elements  in  the  resulting  matrix  to  zero,  we  find  that  a  solution  for  T 
is  obtained  if  we  choose  u  and  v  as  follows? 


(1.2.4) 


U  =  2/  ’  6  +  +  } 
v=--^-{a-6  +  \T(  a  -  6)s  +  4(3/)  . 

^ p 


We  find  moreover  by  back  substitution  that 

Af  =  -|  {a  +  5  -  \T(a  -  S)2  +  437} 

Ag  =  i  (a  +  6  +  -JJa  -  6)2  +  Up7)  . 


(1.2.5) 


n9%iB  3onkJ%lb  earf 


-  5  - 


Let  f(A)  be  defined  by 


00 


(1.2.6) 


f(A)  =  £ 


Ck  A 


00 


where  ^  je^j  is  convergent  for  |A|  <  A^.  Then  using 


-1 


f(A)  =  T  f(TAT  )T  or  more  directly  using  the  Lagrange-Sylvester 
interpolation  formula  (see  e.g.  [15]  P.  77 )  we  find  that 


(1.2.7) 


f  (A)  = 


(A  -  AJ  f(Aj  -  (A  -  A  )  f(A1) 


A  -  A 
2  1 


If  the  eigenvalues  of  A  are  not  distinct,  an  arbitrarily 
small  change  in  the  elements  of  A  can  make  them  distinct.  Hence  by 
continuity  of  ( 1.2.7)  equation  (1.2.7)  holds  whether  the  eigenvalues  of 
A  are  distinct  or  not.  For  if  they  are  not  distinct  we  can  treat  (1.2. 7) 
as  if  they  were  and  write  this  equation  in  the  form 


(1.2.8) 


f (A)  =  (A  -  Ax) 


(\)  ~  f(A. 


A1  '  A2 


4 


-  f(A2)  • 


Upon  noting  that  f(A)  is  an  analytic  function  of  the  complex  variable 
A  we  obtain,  taking  the  limit  A^  -*.A^S 

f(A)  =  (A  -  A^l)  f'fAp  -  I  f(Ax)  . 


(1.2.9) 


.JdAl^aib  9^bi  neo  A  '  in.'-  >JU>  nil  i  ?.  -r  H^r 


-  6  - 


1.3  Variations  on  a  Contour 

Let  f(z)  be  an  analytic  function  of  z  and  let  IP  be  a 
contour  such  that  the  integral  on  the  right  of  (I.3.I)  exists.  We 
define 

(1.3.1)  Y7(f(z)'')=  f  |  j-  f(z)  dz|  =  f  |d  f(z)| 

Jp  up 

where  the  contour  IP  lies  in  the  domain  of  analyticity  of  f(z). 

1.4  Formal  Series  and  Asymptotic  Expansions 

We  shall  use  the  definition  of  formal  series  given  in  [4] 

(pages  (115)  and  (l4l)). 

Poincare ' s [38]  definition  of  an  asymptotic  expansion  will  be 
employed.  More  general  definitions  are  given  in  [25,28].  A  formal  series 

00 

~  r-1  -v 

(1.4.1)  F(z)  =  y  F^  z  >  an  mxn  matrix  (k  =  0,1,***) 

k=o 

is  called  an  asymptotic  expansion  of  the  mxn  matrix  function  F(z) 
valid  in  a  given  range  of  values  of  arg  z,  if  for  every  fixed  non¬ 
negative  integer  s,  the  expression 

s 

(1.4.2)  ZS  |f(z)  -  y  Fk  z"k} 

k=o 


The  symbol  ”="  will  here  and  henceforth  mean  a  formal  equality. 


i  ' 


-  7  - 


tends  to  the  mx  n  zero  matrix  as  |z|  -►oo  while  arg  z  remains  in 
the  given  range.  When  this  is  the  case  we  shall  write 

00 

(1.4.3)  F(z)  ~  Fk  z  or  F(z)  F(z) 

k=o 

More  generally,  if  H(z)  is  a  non-singular  m  x  m  matrix  and 

(1.4.4)  H( z)~ 1  F(z)  ~  Fq  +  F1  z"1  +  F2  z’2  + 
we  shall  write 

(1.4.5)  F(z)  H( z)  (Fq  +  Fx  z'1  +  F2  z'2  +  •••). 

The  term  H(z)  Fq  is  called  the  dominant  term,  or  dominant  matrix  of 
the  asymptotic  representation  (1.4.5)  of  F(z). 

The  Borel-Ritt  theorem  (see  e.g.  [13]  page  BIO7)  states  that 
to  any  given  formal  series  and  any  sector  in  the  z-plane  there  exists  a 
function  f(z)  holomorphic  in  the  given  sector  such  that  the  given 
formal  series  is  the  asymptotic  expansion  of  f(z)  as  | z |  — ►  00  in  the 
sector. 

1.5  The  System  of  Differential  Equations 
In  this  thesis  we  are  concerned  with  the  system  of  differential 

equations 

(1»P»1)  4^  «  A(»)W  9  r  integer 


from  the  p@int  el  view  §f  asymptotic  solutions.  The  matri*  A(g)  is 


o  J* 


bnB  xliJj  n  m  x  m  r  JugniB-noo  a  al 


s.tJttw  J  i»  ®v 


' 


° 


. (s )  o  (  .4.-0  rre  i  ff  v trfd 
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n  x  n  and  holomorphic  in  a  domain  JS  extending  to  infinity  in  the 


z-plane.  Moreover,  for  any  positive  integer  s  we  have 


s- 1 


00 


(1.5.2) 


A(z)  = 


V  A  -k 

L  ^  z 

k=o 


+  Ag(z)  z s  /x,  y 


Ak z 


-k 


k=o 


where  the  elements  of  A  (z)  are  bounded  for  all  sufficiently  large 

s 

|z|,  z  in  ,  That  is,  the  matrix  A(z)  need  merely  have  an 
asymptotic  expansion  for  z  in  cQ  ;  we  do  not  require  it  to  have  a 


convergent  Taylor  series  expansion  about  z  =  oo. 


If  r  >  0  equation  (1.5. 1  is  said  to  have  an  irregular  singularity 
of  rank  r  +  1  at  infinity;  if  r  >  -1  equation  (1.5. l)  will  be 
referred  to  as  the  singular  equation. 


An  actual  solution  W(z)  of  (l.5.l)  is  any  non-singular 
n  x  n  matrix  of  holomorphic  functions  which  satisfies  this  equation  for 
all  z  in  JQ  .  If  W(z)  is  an  actual  solution  of  (1.5. l)  we  shall  also 
say  that  W(z)  constitutes  a  full  set,  that  is,  a  linearly  independent 
set,  of  (vector)  solutions . 

'•V 

A  formal  solution  W(z)  of  (1.5-1)  is  a  formal  series  of  the 
type  (1.4.3)  or  (1.4.5)  with  m  =  n  which,  when  substituted  into 
(1.5.1)  as  if  it  were  a  convergent  series,  enables  us  to  uniquely  deter¬ 
mine  all  the  coefficients  of  the  formal  series  when  we  equate  equal 
powers  of  z  ^ .  If  the  dominant  matrix  of  the  formal  series  obtained 
in  this  manner  is  non-singular,  we  shall  say  that  W(z)  constitutes  a 
linearly  independent  set  or  a  full  set  of  n  formal  (vector )  solutions . 


. 


_ 
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1.6  A  Review  of  the  Literature  Concerning;  Formal  and  Actual  Solutions 

Only  a  brief  review  of  the  literature  relevant  to  the  present 
work  will  be  given. 

✓ 

In  1877  Thome  showed  that  the  equation 

y"  +  p(z)y'  +  q(z)y  =  0 

( 1 .6. 3.)  »  » 

p(z)  =  ^  Pk  z"  ,  q(z)  =  y  qk  z"k,  | z |  >  p  >  0 
k=o  k=o 


has  formal  solutions 

00 

'V  W .  Z  p  r-1  . 

(1.6.2)  yj  =  e  J  -  X  «Jk  «*  ;  %  /  0,  j  =  1,2  . 

k=o 

The  formal  series  on  the  right  of  (1.6.2)  do  not  generally  converge. 


In  1885  Fabry  [48]  showed  that  there  exists  a  linearly 

til 

independent  set  of  n  formal  solutions  for  the  n  order  equation 


(1.6.3) 


jn 

p  4-*  +  p 

n  ,  n  n-1  .  n-1 

dz  dz 


,n- 1 
d y 


+  •••  +  pi  dt  +  V  =  0 


where  (k  =  0,l,***,n)  is  an  arbitrary  polynomial  in  z. 


In  1886  Poincare  [38],  using  Laplace  transformations,  showed 
that,  if  the  degree  of  none  of  the  polynomials  P^  (k  =  0,l,*»*,n)  exceeds 
that  of  P^  in  (1.6. 3),  (i.e.  the  rank  1  case;  for  proof  see  e.g.  [4] 
page  I69),  then  the  n  formal  solutions  of  (1.6.3)  obtained  by  Fabry 
are  asymptotic  expansions  of  actual  solutions  of  this  equation. 


■  i\  J  *  :  - 
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Horn  [39]  extended  the  method  of  Poincare  to  the  case  where 
the  coefficients  P^  (k  =  0,l,***,n)  in  (1.6. 3)  have  a  convergent 
power  series  expansion  of  the  form 


(1.6.4) 


kr 

z 


k  =  0, 1,  •  •  •  ,n 


about  z  =  00.  He  showed  moreover  that  in  the  case  when  thfe  roots  of 
the  characteristic  equation 

(1.6.5)  ^  +  an_1  ^  1  +  +  aQ  =  0  (a^  =  Pk  k  =  0,1,  •  ••  ,n) 


are  distinct  the  system  (1.6. 3)  with  P^  (k  =  0,1, •••,n)  given  by  (1.6.4) 
has  a  linearly  independent  set  of  factorial  series  solutions  of  the  form 

00 

1 1.6.6)  C'V‘A) 

k=o 

where  7  and  (k  =  0,1,2,**»)  are  constants.  In  his  treatment  of 
the  equation  (1.6.1)  in  [40]  Horn  uses  the  leading  terms  or  partial 
sums  of  formal  solutions  to  construct  a  differential  equation  which  is 
in  a  certain  sense  close  to  the  given  equation  when  |z|  is  large, 
and  to  prove,  with  the  aid  of  Vol terra  integral  equations,  the  existence 
of  actual  solutions  of  (1.6. l)  for  which  the  formal  solutions  (1.6.2)  are 


* 


This  paper  contains  other  references  of  the  works  of  Horn. 


3  C  '•  (f  .(  )•  I  ‘  0  ->  3 


o  ooi  f»-  ,  i r 
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asymptotic  expansions.* 

The  existence  of  n  linearly  independent  actual  solutions  of 
(I.6.3)  f°r  which  the  formal  solutions  of  this  equation  are  asymptotic 
expansions,  in  the  case  when  z  is  real  and  the  coefficients 
I?k  (k  =  0,l,**#,n)  merely  have  asymptotic  expansions,  was  proved  by 
Sternberg  [44]  in  1920. 

In  1909  [^5]  and  again  in  1913  [50  ]  G.D.  Birkhoff  attempted  to 

reduce  the  equation  (1,5.1)  in  the  case  when  A(z)  has  a  convergent 

) 

power  series  expansion  andt’r  >  <D  >  to  the  form 

s 

(1.6.7)  Y'  =  (**  ^  c 

j=o 

with  b  =  r  +  1,  by  making  the  transformation  X  =  B(z)Y,  where  the 
n* 2  elements  of  B(z)  are  all  regular  for  |z|  sufficiently  large  with 
the  determinant  B(oo)  not  zero.  His  proof  is  correct  for  the  case  when 
Aq  in  (1.5.1)  has  distinct  eigenvalues  but  it  was  pointed  out  by  F.R. 
Gantmacher  [49]  that  Birkhoff' s  Theorem  (reproduced  with  proof  in  Ince 
[23]  page  470)  is  false  and  that  it  will  in  general  be  necessary  to  take 
the  integer  sy  (  (  (l„. 6. 7)  above)  -  larger  than  r  +  1. 

In  1933  W.J.  Trjitzinsky  treated  the  equation  (1.6. 3)  f°r  the 
case  when  the  coefficients  are  holomorphic  and  need  merely  have  an 

* 

For  an  earlier  result  along  this  line  see  Liouville,  J.  de  Math  (l) 

2  (1838)  p.  19  and  (1)  3  (1838),  p.  565. 
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asymptotic  expansion  (as  |z|  -*-oo)  in  some  domain  extending  to  infinity. 
He  established  the  existence  of  a  full  set  of  actual  solutions  for  which 
the  formal  solutions  of  this  equation  are  asymptotic  expansions.  His 
theorems  of  a  similar  nature  concerning  the  more  general  system  (1.5.1) 
are  also  correct,  although  the  existence  of  a  linearly  independent  formal 
series  solution  of  (1.5.1)  had  at  this  time  not  yet  been  fully  established. 
Trjitzinsky  [ 3 1 ]  (1937)  also  extended  the  work  of  Horn  [39]  in  obtaining 
convergent  factorial  series  solutions  in  certain  cases  when  the  roots  of 
the  characteristic  equation  (1.6. 5)  are  not  necessarily  distinct. 

In  1935  Hokuhara  proved  the  existence  of  a  full  set  of  formal 
solutions  of  (1.5.1)  (Ill.  [22]). 

The  paper  of  H.L.  Turrittin  [32]  (1955)  contains  the  most  com¬ 
plete  treatment  to  date  of  the  reduction  of  the  general  system  (1.5. l) 
to  canonical  form.  This  result  is  explicitly  stated  in  Chapter  III  of 
this  thesis.  In  his  paper  of  1955*  Turrittin  assumed  the  coefficient 
matrix  in  (1.5.1)  to  have  a  convergent  power  series  expansion  about 
z  =  00;  this  enabled  him  to  replace  certain  of  the  divergent  asymptotic 
series  which  represent  solutions  of  (I.5.I)  by  convergent  factorial  series. 
He  thereby  extended  the  work  of  Trjitzinsky  [31].  Two  other  papers  of 
Turrittin  extend  the  work  of  Birkhoff  [46,45].  These  concern  (a)  [34]  the 
reduction  of  the  system  (1.5. l)  to  the  Birkhoff  canonical  form  -  in 
this  paper  he  proved  the  existence  of  the  integer  s  (possibly  larger 
than  r  +  l)  in  (1.6. 7);  and  (b)  [33]  he  replaced  the  n  x  n  system 
(1.5.1)  rank  (r  +  l)  by  an  n(r  +  l)  x  n(r  +  l)  system  of  rank  1. 
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The  book  [47]  edited  by  C.H.  Wilcox  contains  a  more  complete 
list  of  references  than  those  mentioned  above. 

1.7  A  Brief  Review  of  the  Literature  Concerning  Error  Bounds 

In  1890  H.  Weber  obtained  error  bounds  for  asymptotic  approxi¬ 
mations  to  the  Hankel  functions.  By  his  method  (summarized  in  [30]  page 
194)  he  used  an  integral  representation  for  the  Hankel  function,  together 
with  solving  the  differential  equation  for  the  partial  sum  (which  is  used 
to  approximate  the  Hankel  function)  of  an  asymptotic  expansion. 

In  1912  0.  Blumenthal  [21]  obtained  error  bounds  for  the 
Liouvi lie -Green  (or  WKB)  approximation 

(1.7.1)  W  =  A  q  4  exp  ^P  u  q2  dx^+  B  q  u  exp  ( -u  P  q2  dx) 

for  the  solution  to  the  equation 

(1.7*2)  w"  =  u2  q(x)  w 

where  u  is  a  large  parameter.  His  results  are  valid  for  x  in  an 
arbitrary  finite  interval. 

In  I96I  F.W.J.  Olver  [14]  sharpened  the  results  of  Blumenthal, 
and  generalized  them  for  x  a  complex  variable  in  some  (possibly 
infinite)  complex  domain.  In  I963-64  Olver  [26,27]  found  error  bounds  for 
first  approximations  and  asymptotic  expansions  in  turning  point  problems. 
In  1964  he  [35]  wrote  a  general  paper  on  error  bounds  for  asymptotic 
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expansions,  with  application  to  Bessel  functions  of  large  argument.  In 
1964  he  [6]  wrote  a  paper  on  the  asymptotic  solution  of  second  order 
differential  equations  having  an  irregular  singularity  of  rank  1  with  an 
application  to  Whittaker  functions  of  large  argument.  In  1964  F.W.J, 
Olver  and  F.  Stenger  [7]  found  error  bounds  for  asymptotic  solutions 
of  second  order  differential  equations  having  an  irregular  singularity  of 
arbitrary  rank. 


In  all  the  above  cases,  Olver' s  procedure  is  to  express  the 
error  term  (i.e.  the  difference  between  an  actual  solution  and  the 
partial  sum  of  a  formal  solution)  in  the  form  of  a  non -homogeneous  second 
order  differential  equation  and  then  to  obtain  accurate  bounds  for  the 
solutions  of  this  equation. 

In  his  book  of  I962  [29]  H.  Jeffreys  gave  a  method  of  obtaining 
error  terms  for  asymptotic  expansions  obtainable  by  use  of  Watson’s 
Lemma  ([29]  page  14). 


Richard  [42]  in  1963  found  error  bounds  for  the  approximations 
sin  p(x),  cos  p(x)  ^3(x)  =  \l  1  +  ^(£)  to  the  solutions  of  the 

pOO 

equation  y"  +  [l  +  ^(x)  +  r(x))  y  =  0  where  /  j^'(x)|dx  <  00, 


J 


poo 

/  |r(x)|  dx  <  00,  and  lim  [1  +  ^(x)]  >  0. 


x  — 00 


Useful  results  pertaining  to  error  bounds  for  approximations 
to  solutions  of  differential  equations  arise  out  of  stability  theory. 


' 
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In  particular  we  mention  here  Gronwalls  [ 17 ]  (or  Bellman's  [2])  inequality. 
The  paper  of  Li  Yue  Shing  contains  (Bernoulli  equation)  norm  bounds  for 
the  solution  of  =  A(z)X  +  where  ||^(z,X)||  <  P(z)  ||X||  +  Q(z)  ||X||a, 

a  >  0.  The  paper  of  Willett  and  Wong  [44]  contains  generalizations  of 
Gronwall's  inequality;  this  paper  also  contains  a  wealth  of  references  on 
this  inequality. 


1.8  The  Problem 


An  asymptotic  expansion  of  a  function  f(z)  with  respect  to 
large  values  of  the  argument  z  can  usually  be  shown  to  hold  uniformly 
in  some  sector  a  +  5  <  arg  z  <  3  -  6,  for  some  positive  number  6  in 
0  <  6  <  ■§■(  3  -  a) .  This  implies  for  example  that  for  all  z  sufficiently 
large  and  in  the  sector  a  +  5  <  arg  z  <  3  -  6  we  have 

(1.8.1)  |f(z)/g(z)  -  aQ  -  z'1-  ...  -  ak  z"k|  <  K  |z|"k_1 


for  some  constant  K  depending  only  on  6.  Such  a  bound  gives  one  a 

very  limited  idea  concerning  the  actual  magnitude  of  error.  It  is 

known  for  example  in  the  case  of  the  Bessel  functions  that  the  expansion  [35] 


(1.8.2) 


J  (z) 
v  ' 


'hi  {■ 
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sin  cp  (  — 


!s 


for  fixed  v  ,  where 
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holds  uniformly  for  large  z  in  |arg  z|  <  Jt  -  0  <  5  <  jt. 

The  accuracy  yielded  by  using  a  finite  number  of  terms  of  each  of  the 
series  on  the  right  of  (1.8.2)  decreases  steadily  as  arg  z  ->t  it. 

That  is,  the  constant  K  on  the  right  of  (1.8.1)  becomes  increasingly 
large  as  S  -*■  0 . 

Now  in  the  case  when  a  computer  evaluates  a  function  by  its 
Taylor  expansion  he  is  warned  of  inaccuracies  incurred  near  the  boundary 
of  the  region  of  convergence  since  the  magnitude  of  successive  terms 
of  the  series  decrease  less  slowly  as  the  radius  of  convergence  is 
approached.  There  is  no  such  warning  in  the  case  of  an  asymptotic 
expansion  because  the  rate  of  reduction  of  the  magnitude  of  successive 
terms  in  the  series  is  independent  of  the  proximity  of  the  boundary. 

In  deriving  asymptotic  expansions,  "exponentially  small" 
contributions  to  the  final  asymptotic  series  are  discarded.  For 
example  [35]  it  is  readily  shown  by  repeated  integration  by  parts  that 
the  integral 

(1.8.1*)  I(n)  =  f*  £2|J2_£  dt 

^  O 

* 

has  the  Poincare  expansion 

(1.8.5)  l(n)  ^  (-l)n  '•*')  • 

^  n  n  n 

Evaluating  (1.8.4)  by  direct  numerical  quadrature  for  n  =  10  we  obtain 
I(10)  =  -  0.0004558  to  seven  decimals.  On  the  other  hand  using  (1.8.5) 
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for  n  =  10  we  obtain  I(10)  ^  (0.0005518  -  0.0000048  +  0.00000  01  -  •••) 

=  -0.000527¥\,  an  incorrect  result.  If  however  we  write  (1.8.4)  in  the 
form 


(1.8.6) 


r>° 0  .  pco 

l(n)  =  T  £2|_n_E  dt  .  f 
Jo  *  +  1  Jn 


cos  n  t 

"FTT 


dt,  we  get 


1 


jr 


-n 

e  + 


(-1  r1 


6 


n 


) 


where  the  first  integral  on  the  right  of  (1.8.6)  equals  ^  jt  e  n  and 
where  the  second  has  the  same  asymptotic  expansion  as  (1,8.6).  Using 
the  "complete”  expansion  on  the  right  of  (1.8.6)  to  evaluate  1(10),  we 
obtain  the  correct  result  -  0.0004558.  The  point  n  =  10  is  not  a 
special  case;  (1.8.6)  gives  a  better  result  than  (1.8.5)  for  all  finite 
n.  There  is  however  no  self-evident  reason  why  this  should  be  so. 


Asymptotic  analysts  have  long  recognized  that  in  general  it 
is  possible  to  achieve  only  limited  (though  often  very  high)  accuracy 
when  evaluating  a  function  by  its  asymptotic  expansion.  In  attempting 
to  overcome  this  difficulty  they  have  tried  summing  divergent  series  and 
obtaining  "complete"  asymptotic  expansions.  The  problem  of  finding 
convergent  sums  of  divergent  series  solutions  has  not  yet  been  completely 
solved  even  for  the  case  where  A(z)  in  (1.5. l)  has  a  convergent  power 
series  expansion  about  z  =  00.  Finding  convergent  sums  of  divergent 
series  solutions  when  A(z)  merely  has  an  asymptotic  expansion  in  certain 
.  i: v ■  I,- : .■  <•'£  i.n  i.  i. r»  .■  Ly  O'  1  i  >.  ■  ■ . 
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In  the  sense  of  Miller  and  Watson  [24]. 
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sectors  of  infinity  has  not  at  all  been  attempted. 

If  however  we  demand  somewhat  less,  namely,  a  bound  on  the 
difference  between  an  actual  solution  and  the  partial  sum  of  an  asymptotic 
expansion  for  the  actual  solution,  we  obtain  a  problem  which  does  not 
seem  to  be  as  difficult  to  solve.  Moreover  since  in  most  cases  it  is 
much  easier  to  evaluate  a  partial  sum  of  an  asymptotic  expansion  than  the 
actual  solution,  or  a  convergent  factorial  series  solution,  such  bounds 
are  of  considerable  value. 


1.9  Layout  of  the  Thesis 

In  this  thesis  the  method  of  Olver  (see  e.g.  [6])  is  extended 
to  obtaining  error  bounds  for  a  system  of  a  first  order  equation  of 
arbitrary  rank.  The  method  in  [6]  (or  [7])  does  not  carry  over  to  the 
arbitrary  system  directly,  since  a  Liouville-Green  (or  WKB)  transformation 
is  not  known  for  the  general  system  of  n  first  order  equations.  Other 
difficulties  arise  because  of  the  appearance  of  more  than  two  eigenvalues 
in  the  lead  matrix  of  the  general  system,  and  because  of  two  or  more  of 
these  eigenvalues  being  the  same. 

As  we  have  already  mentioned,  there  exists  a  transformation  which 
transforms  the  general  system  (1.5.1)  to  a  canonical  form  (see  Theorem 
2.4.1).  We  shall  find  solutions  and  error  bounds  for  solutions  of  the  trans¬ 
formed  system.  It  is  possible  that  this  transformed  system  has  a  regular 
point,  a  regular  singular  point  or  an  irregular  singular  point  at 

# 

See  however,  the  paper  of  H.L.  Turrittin  in  [47]  for  a  different 
approach  to  this  problem. 
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infinity.  As  an  aid  in  overcoming  difficulties  we  have  started  with 
the  simplest  case  and  progressively  proceeded  to  the  most  difficult  and 
complex  case.  Thus  we  considered  first  the  case  of  a  regular  point. 

The  case  of  a  regular  singular  point  is  more  difficult  in  that 
we  need  to  choose  the  canonical  form  (Section  2.2)  so  as  to  enable  us  to 
easily  bound  any  particular  j'th  error  vector  solution.  Additional 
restrictions  are  required  on  the  paths  of  integration  in  the  Volterra 
integral  equations  defining  the  solutions. 

When  all  eigenvalues  of  Aq  are  distinct,  it  is  possible,  in 

the  case  of  an  irregular  singular  point  to  find  a  vector  solution  (both 

formal  and  actual)  corresponding  to  each  eigenvalue.  The  analysis  for 

obtaining  error  bounds  is  considerably  simpler  in  the  case  when  the  solution 

* 

corresponds  to  an  eigenvalue  which  is  an  extreme  point  of  the  smallest 
strictly  convex  closed  polygon  containing  all  the  eigenvalues  of  A^,  and 
is  considered  first.  In  this  case  the  error  can  be  expressed  in  the 
form  of  a  single  Volterra  (vector)  integral  equation.  In  the  case  when 
the  vector  solution  corresponds  to  an  eigenvalue  which  is  an  interior 
point  of  the  polygon,  the  error  can  be  expressed  in  the  form  of  a 
pair  of  simultaneous  Volterra  (vector)  integral  equations. 

In  all  of  the  above  cases  we  have  assumed  that  end-points  of 
integration  can  be  chosen  at  infinity.  This  assumption  is  dropped  when 
we  consider  the  general  case,  which  includes  the  above  cases  as  well  as  all 
others . 

By  extreme  point  we  mean  vertex. 
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In  Chapter  II  of  the  thesis  we  transform  the  general  system  to 
canonical  form,  and  explicitly  construct  formal  solutions  for  the 
transformed  system.  We  give  an  explicit  construction  of  the  transformation 
to  canonical  form  for  the  case  of  a  regular  singular  point  and  for  the  irregular 
singular  case  when  all  eigenvalues  of  the  lead  matrix  are  distinct.  For 
Turrittin's  theorem  [32]  concerning  transformation  of  the  general 
system  to  canonical  form  we  state  a  slight  modification  which  enables  us 
to  easily  bound  the  solutions  of  integral  equations. 

In  Chapter  III  of  the  thesis  we  obtain  actual  solutions  and 
error  bounds  for  formal  partial  sum  approximations  to  actual  solutions. 

The  error  bounds  are  given  in  two  forms:  as  a  norm  bound  and  as  a  vector 
bound. 


In  Appendix  A  of  the  thesis  we  solve  the  general  singular  n  =  2 
case  to  illustrate  the  results  of  the  thesis. 

The  proofs  concerning  existence  and  uniqueness  of  solutions  of 
integral  equations  are  contained  in  Appendix  B. 

Appendix  C  contains  some  extensions  of  Gronwalls  [17]  inequality  to 
matrices,  to  Fredholm  integral  equations,  and  to  simultaneous  Volterra 
integral  equations  These  extensions  are  an  aid  in  arriving  at  the  results 
of  theorems  in  Section  3.3.5  and  3.4.5  of  the  thesis. 

Each  chapter  and  each  appendix  begins  with  a  short  synopsis 
which  summarizes  its  results. 


icl  1mm  3aioq  xeXu,  uie  xaj 1/307  a  io  erta  x<  1  ixoSt  I*oiaon«D  ox 


>rv  V  \  :  t  bgflljT, JrtOO  J»XE  SnolXfiijrpfl  I»Xy<>3ni 


l  •  ■  *  •  ;■  ,< 


auooaaJiuaicj  o3  hr it  t«fKvUftft>pe  Ijvt£b*a1  mlcrfb®*!*  03  „  t«dli3«e* 


eXXoasx  aria  3b  gnlvixxfi,  oJ  blt>  an  si*  noataiic  t-  *earfT  iuojlaupa  imsal'ik 


CHAPTER  II 


TRANSFORMATIONS  ON  THE  ORIGINAL  SYSTEM  AND  FORMAL 
SOLUTIONS  FOR  THE  TRANSFORMED  SYSTEM 


The  analysis  is  considerably  simplified  after  transforming  a 
given  system  of  differential  equations  to  canonical  form.  Besides  trans¬ 
forming  systems  of  differential  equations  to  canonical  form  we  also 
obtain  formal  solutions  of  the  transformed  system  in  this  chapter. 

Transformations  which  transform  a  given  system  to  canonical 
form  are  explicitly  described  for  the  case  of  a  regular  singular  point 
and  for  the  case  of  an  irregular  singular  point  when  the  eigenvalues  of 
the  lead  coefficient  matrix  are  distinct.  There  is  as  yet  no  complete 
constructive  proof  of  Turrittin's  theorem  concerning  transformation  of 
the  very  complicated  general  singular  system  to  canonical  form.  We  have 
however  proved  modifications  of  this  theorem  corresponding  to  our 
modifications  of  the  usual  canonical  form  for  the  case  of  a  regular 
singular  point.  These  modifications  are  an  aid  in  bounding  the  solutions 
of  integral  equations  in  the  next  Chapter. 


2 . 1  The  Case  of  a  Regular  Point 

In  this  case  the  given  system  of  differential  equations  takes 
the  special  form 

(2.1.1)  ~  =  z'2  A(z)W 

where  for  each  arbitrary  positive  integer  s 

s-1 

(2.1.2)  A(z)  =  ^  A^  z  k  +  Ag(z)  z  S  ;  A^  constant  (k  =  0,1,...) 

k=o 

and  where  the  n2  elements  of  A  (z)  are  bounded  for  z  in  some  domain 

sv  ' 

(&  extending  to  infinity.  That  is,  we  assume  that  the  elements  of  A  (z) 

s 

remain  bounded  as  |z|  -*•  00  in  o 9  ,  and  that  A(z)  has  the  properties 


described  in  Section  1.5. 


'  li<  I  1  98l:o  ->rf  l  l<  t SStoh  /ill  >1J  :X9  dlB  Mo'! 


( . . .  tI  ,0 


22 


In  the  case  of  a  regular  point  we  need  not  make  any  preliminary 
transformations ;  the  system  (2. 1 . 1 )-(2. 1 .2)  is  ready  for  formal  soltuion 
as  it  stands. 


To  obtain  a  formal  solution  we  substitute 


00 


(2.1.3) 


W(z)  =  ^  \  z'\  u0  =  1 

k=o 


into  (2.1.1)  with  the  coefficient  matrix  A(z)  formally  expanded;  i.e. 


00 


00 


(2.1.4) 


y  <.-)k 


-k+1 


-k 


or 


k=o 


k=o  i-o 


00 


(2.1.5) 


y  h + ^  Vi + y  ^ 


k=o 


-k  /v 
z  =  0 


i-o 


On  equating  equal  powers  of  z  to  zero  in  (2.1.5)  we  obtain 


Theorem  2.1.1  The  system  (2. 1 . 1 )- (2. 1 .2)  has  a  formal  independent 
series  solution  (2.1.3)  with  Uq  =  I,  and 

k 


(2.1.6) 


X  Ai  Vi!  /  <k  +  !)•  k  =  0,1,-..  . 

1-0 


It  is  of  course  well  known  (see  e.g.  [4]  Chapter  III)  that  in 
the  case  when  A(z)  has  a  convergent  power  series  expansion  in  a  neighbor¬ 
hood  of  z  =  oo  then  the  formal  series  (2.1. 3)  converges  in  this  same 

*  /V 

The  sign  ”  =  "  indicates  a  formal  equality. 
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neighborhood  of  infinity  and  is  an  actual  solution  of  this  equation. 

2 . 2  The  Case  of  a  Regular  Singular  Point 

The  given  system  now  takes  the  form 

(2.2.1)  g  =  z'1  A(z)  X 

where  A(z)  is  given  by  (2.1.2)  and  has  the  properties  described  in  the 
previous  section. 

We  shall  reduce  the  system  (2.2.1)  to  canonical  form  by  a  method 
which  closely  resembles  that  of  H.L.  Turrit in  [32]  (pp.  Jl,  32  and  37) 
but  with  a  slight  modification  allowing  for  actually  carrying  out  the 
computations . 

The  normalizing  transformation 

(2.2.2)  X  =  T  W 

where  T  is  a  constant  non-singular  matrix  is  made  in  (2.2.1),  changing 
(2.2.1)  to  the  form 

(2.2.3)  =  z  1  B(z)W  where  B(z)  =  T  ^  A(z)T, 

so  that  for  each  arbitrary  positive  integer  s 

s-1 

(2.2. It)  B(z)  =  Bk  Z'k  +  Bs(z)  z's 

k=o 


n.  to  )rU  ■■  ?>)(  3  won  1  r  iv  1  f*  #T 
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where  B,  =  T  1  A  T  (k  =  0,1,***)>  B  (z)  =  T  ^  A  (z)  T.  It  follows 
K  K  s  s 

that  since  the  elements  of  ^s(z)  are  bounded  for  all  z  in  a  particular 

domain  extending  to  infinity,  the  elements  of  B  (z)  are  similarly 

s 

bounded  for  z  in  this  domain.  It  is  further  presumed  that  T  has 

been  chosen  so  that  B  takes  on  the  classical  Jordan  canonical  form 

o 


(2.2.5) 


B  = 
o 


M, 


0 


0 


0 


M„ 


0 


0 


0 


0 


M. 


I 


with  (i  =  l,2,***,i)  a  p^  x  p^  M-  =  n>)  matr:*-x  t*ie  f°rm 

i=l 


(2.2.6) 


M.  =  A. I  +  J  . 
i  i  pt  |it 


In  (2.2.6)  A.  is  a  complex  number,  I  is  a  p.  x  p, 
l  ^i  i  i 


unit  matrix 


and  J  is  a  p„  x  p„  matrix  with  zeros  everywhere  except  in  its  first 

Pi  11 

subdiagonal  where  it  may  have  ones  or  a  mixture  of  zeros  and  ones. 


If  0  <  Re  X,  <  1,  -1  <  Re  A,  <  0  and  A,  -  A  .  J.  1 

=  i  =  j  i  j  / 

(i  =  l,2,#,,,i,  i  ^  j)  and  if  the  transformation  W  =  T(z)X  where  T(z) 

is  a  diagonal  matrix  with  i'th  diagonal  block  I  (i  ^  j)  but  with 

-1 

j ! th  diagonal  block  z  I  completely  fills  B  with  zeros  we  can 
then  proceed  directly  to  Theorem  2.2.1.  If  these  conditions  are  not 


satisfied  we  may  assume  without  loss  of  generality  that 
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Draw  vertical  lines  through  the  points  X  =  k  (k  an  integer) 
in  the  complex  A  plane.  If  A^  lies  in  the  strip  k.  <  Re  A  <  k.  +  1 
we  associate  the  integer  k^  with  A^  (i  =  1,2,  •••,.£)•  It  is  clear  that 

k!  i  k2  i  ***  5 


With  'k  =  k^  -  k^,  we  next  make  the  sequence  of  zero  inducing 
transformations 


(2.2.7) 


\-l  =  (I  + 


z"k  T)  2  (k  =  1,2, •••,»);  Z  =  W, 


on  (2.2.5).  Thus  for  example,  if 


(2.2.8) 


dZ 


k-1  -1 

=  z 


dz 


Bo  +  B^1  +  ...  +  Bkz'k  +  Bk+1(z)  z-k'l 


'k-1 


where,  whenever  possible  all  blocks  B^,  B^, 


,  B^  ^  have  already 


been  made  to  vanish,  and  B^+^(z)  is  defined  in  a  domain  (fQ  as  on  page  &/ 

except  that  |z|  may  need  to  be  taken  larger  than  before,  to  exclude 

-k 

the  zeros  of  the  determinant  det(l  +  z  T^)*  tken 


dZ*  -1  /  -1 

<2-2-9)  ZT  =  Z  (Bo  +  V  + 


1-k  n  -k  *  .  »  -k-l\ 

+  Bk-iz  +  V  +  Vi(z)  z  K 


where 


(2.2.10) 


C.  =  B.  +  B  T,  -  T.  B  +  k  T. 
k  k  o  k  k  o  k 


and  where  the  elements  of  B^  ^(z)  are  defined  similarly  as  those  of 
Bk+i(z>- 


' 


1 


.. 
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We  divide  the  constant  matrices  C,  ,  B.  and  T,  into  smaller 

k  k  lc 

blocks  in  the  same  way  that  is  subdivided  in  (2.2.5).  After  this 

subdivision  is  made  we  denote  the  i'th  row  of  blocks  and  the  j ' th  column 
of  blocks  respectively  by  C  Bijk  anc^  ^ijk  =  1>2,  .  .  .  ,£ ) . 

Equation  (2.2.10)  for  C..,  then  becomes 

1  j  K 

(2.2.11)  C.  ..  =  B.  ..  +  (A.  I  +  J  )  T  -  T  (A  .1  +  J  )  +  k  T.  ; 

ljk  ijk  l  ni  (ii/  ijk  ljk  j  }i  1 1/  ijk 

if  A.  -  A.  J.  k,  elements  of  T...  can  be  computed  (for  further  details 
j  l  '  ijk  r  v 

see  [52]  p.  33)  starting  with  the  right  hand  vector  and  calculating  from 

top  to  bottom  and  in  this  manner  successively  computing  the  vectors  from 

right  to  left  so  that  C..,  =0.  If  A.  -  A,  =  k,  C  cannot  in 

ijk  j  i  ijk 

general  be  made  to  vanish,  and  we  set  T..,  =  0. 

ijk 

Proceeding  in  this  manner  to  k  =  k,  we  fill  all  blocks  of 
Bn ,  B  ,  ...,  B  with  zeros,  except  those  (i,j)'th  (i  <  j)  for  which 
Aj  -  A^  =  k  =  kj  -  k^.  We  emphasize  in  particular  that  all  blocks  on 
and  below  the  diagonal  of  B  ,  B  ,  . . . ,  B  can  be  filled  with  zeros  in 

1  c  K 

this  manner. 


The  root  equalizing  transformation 


(2.2.12) 


Z  = 

K 


ki  _ 

Z  I 


L 


Hi 


0 


\ 
z  I 


H, 


0 


'  . 
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applied  to  (2.2.9)  with  k  =  k  changes  this  equation  to 


(2.2.13) 


dY  -1  ,  *  **  .  x  -l. 

—  =  z  (B  +  B.  (z)  z  l)Y 
dz  v  o  1  ' 


* 


where  B  is  zero  below  the  diagonal.  The  matrices  in  the  (i,i)'th 
o 

block  of  (i.e.  with  respect  to  B^  in  (2.2.5)~(2.2.6) )  are  changed 

to  (Ay  -  k.)l  +  J  (i  -  1,2, *  * ' , i)  while  a  non-zero  (i,j)'th  block 

1  1  |J._.  u. 

x.  1 

(i  <  j)  will  be  (with  respect  to  equation  (2.2.9)  with  k  =  k)  if 

ljK 

and  only  if  A.  -  A.  =  k,  a  positive  integer.  Since  the  transformation 

J  1 

(2.2.12)  merely  multiplies  each  block  of  B(z)  by  some  integral  power 


of  z 


B. 


** 


which  coincides  with 


^  (z)  is  holomorphic  in  a  domain  ofy 
the  domain  in  which  B^(z)  (page  £l)  is  holomorphic  for  all  |zj 

larger  than  the  largest  zero  of  any  of  the  determinants  of  the  sequence 
of  zero  inducing  transformations  used.  Moreover  B^  (z)  has  an  expansion 
exactly  similar  to  B^(z)  defined  in  equation  (2. 2. A). 


If  some  of  the  differences  A^  -  A^s  (i,j  =  l,2,***,i5)  were 
integers  at  the  outset  another  normalizing  transformation  is  applied  to 
(2.2.13)  to  bring  (2.2. 13)  to  the  form  (2.2.3)  where  Bq  has  the  form 
(2.2.5)  with  the  integer  &  smaller  than  before.  The  new  system, 
regardless  of  whether  or  not  A^  -  A^  was  an  integer  at  the  outset, 
will  have  all  its  eigenvalues  in  0  <  Re  A.  <  1, 


For  later  convenience  we  apply  another  zero  inducing  trans¬ 
formation  which  makes  zero,  and  lastly,  we  follow  this  up  by  a  root 

equalizing  transformation  which  puts  the  real  part  of  a  particular  A 
in  -1  <  Re  A^  <  0  while  for  i^j,  0  <  Re  A^  <  1,  This  last  root 


it  >loo  ''  e  IqlSli  1  \'  ...  u  t  i  (  .S.S) 
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equalizing  transformation  takes  the  form  W  =  T(z)X  where  T(z)  has 
the  form  (2.2.12)  with  =  0  (i  =  l,2,...,i5;  i  ^  j)  and  k  =  1. 

Theorem  2.2.1  The  system  (2.2.1)  can  be  transformed  to  the  system 
(2.2.3)  with  B(z)  as  in  (2.2.4)  and  Bq  a_s  in  (2.2.5)-(2.2.6)  by  a 
finite  sequence  of  normalizing,  zero  inducing  and  root  equalizing  trans¬ 
formations  as  described  above .  Moreover,  corresponding  to  any  fixed  integer 

j  in  [ 1 , Z  ]  we  can  achieve  0  <  Re  A  <  1 ,  -1  <  Re  A,  <  0  and 

^  —  J 

A ^  -  A ^  jL  1  (i  =  1 , 2 , . . . ,  i ;  i  ^  j).  The  elements  of  B ( z )  are  holo- 
morphic  and  bounded  in  a  domain  r$  which  coincides  with  the  domain  in  which 
A(z)  is  holomorphic  for  all  |z)  sufficiently  large. 

Theorem  2.2.2  The  system  (2.2.3)  with  B(z)  as  described  in  the 
above  theorem  has  a  formal  independent  series  solution 

00 

(2.2.14)  W( z)  =  (7  z  ^  exP  [B0  l°g  Z1 

k=o 

with  Uq  =  I  and  given  by  the  unique  solution  of 

k 

(2.2.13)  kU.  +  BU.-U.B  =  •  )  B .  U.  . 

kokko  i  k-i 

i= 1 

Proof,  Substituting  (2.2.14)  into  (2.2.3),  cancelling  exp  [Bq  log  z] 

on  each  side,  using  (2.2.4),  expanding  and  equating  equal  powers  of  z  we 

obtain  (2.2.15).  It  is  proved  in  detail  by  Turrittin  in  [32]  p.  33  that 

with  U  =  I,  (2.2.15)  uniquely  defines  U  ,  k  >  1.  Moreover,  the  j ' th 
o  Ic  — 

column  of  blocks  may  be  solved  independently  of  the  other  blocks;  the 
successive  solution  of  columns  in  a  particular  column  of  blocks  being 
carried  out  from  right  to  left,  and  the  elements  of  a  particular  column 


X 


on*o  «(<  .S.2)  o3iii  (iI.S.2)  8«J 
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bein^  solved  from  top  down. 

It  is  well  known  (see  e.g.  [4]  Chapter  iv)  that  if  A(z)  has  a 
series  expansion  convergent  in  a  neighborhood  of  z  =  oo,  then  the  series 
in  round  brackets  on  the  right  of  (2.2.14)  also  converges  in  this  same 
neighborhood  of  z  =  oo,  and  hence  in  this  case  the  formal  solution 
(2.2.14)  is  an  actual  solution  of  (2.2.3). 


We  take  this  opportunity  to  make  a  remark  concerning  a  possible 
economizing  in  computing  successive  coefficients.  Suppose  we  are 
interested  in  knowing  the  coefficients  of  a  particular  vector  solution  - 

r\j  r\j 


say  the 


v'th  column 


in  the  j  '  th  column  of  blocks  W  .  (z)  ()f  W(z)  . 


Let  J  be  a  [i  x  |i  matrix  with  zeros  everywhere  except  in 
its  first  sub-diagonal  where  it  has  a  full  set  of  ones.  Then  by  [5]  page  32 

0  0  •  •  •  0 

1  0  •  •  •  0 

log  z  1  ...  0 

(log  z)11'1  (log  z)u"2  ...  x 

.  (U  -  1):  (n  -  2)1 


(2.2.16)  exp[J  log  z]  = 


log  z 

U°g  ZY 

2: 


Hence  if 
v  +  1 '  th. 


has  a  one  in  its  v'th  column  as  well  as  in  its 
(p  -  l)’th  but  not  in  its  p'th,  where  p  < 


then 


(2.2.17) 


ONlS 


. 
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S 

where  U.,  is  the  s'th  column  vector  in  the  i 1 th  column  of  blocks  of 
jk 

U^.  Thus  it  is  only  necessary  to  compute  the  v' th,  v+  l'th,  •••,  n'th 
column  at  each  stage,  when  computing  the  formal  solution. 

2.3  The  Irregular  Singular  Case  When  the  Eigenvalues 

of  A  are  Distinct 
—  o  - 

The  given  system  now  takes  the  form 

(2.3.1)  =  zr  A(z)  X 

where  r  is  a  non-negative  integer  and  where  A(z)  has  the  properties 
described  at  the  beginning  of  Section  2.1. 

Theorem  2,3.1  If  Aq  has  n  distinct  eigenvalues  then  there 
exist  constant  matrices  T  *  T_ ,  • • • ,  T  .  such  that  when  the  trans- 

-  o>  r+i  - 

formation 

r+1 

(2.3.2)  X  =  T(z)  W  =  (Y  Tkz'k)w 

k=o 

is  applied  to  (2.3. l)  we  obtain 

(2.3.3)  ^|=zrC(z)W 

where  for  each  arbitrary  positive  integer  s 

r+1  s-1 

(2.3.4)  C(z)  =  zr+1  y  Ck  z"k  +  z“2  Y  Bk  z'k  +  Bs(z)  z‘S  , 

k=o  k=o 

and  where  the  elements  of  B  (z)  are  holomorphic  and  bounded  for  all 
| z |  sufficiently  large,  z  €  . 


i{  ,*»•  +  v  ,/fJV  &il3  fiJirqtoc.D  03  yn<8  assn  \lno  it  it  *urfT 

f  .S 


' 

■ 

♦  I.S  nolJD^c  lo  ^nlanlsaii  s.-fj  SB  badJtioeab 


j  w  a  a » 


^5£H  24i  22S!2  1£*  42H5  u*  .•••  .1  <0T  aa^a*  Zi^lSUZt  2Sh 52 


■ 


-  31  - 


In  (2.3.4)  each  matrix  ( k  =  0 , 1  *  •  •  •  , r  +  1 )  _is  diagonal . 


Proof.  We  shall  prove  this  important  theorem  here  (for  the 

original  proof  see  Birkhoff  [46])  in  order  to  furnish  an  explicit 

algorithm  for  computing  the  coefficients  (k  =  0,  1,  •  •  •  ,r  +  l) .  By 

the  analysis  of  the  previous  section  we  suppose  that  the  transformation 

(2.2.2)  has  already  been  made  and  that  the  lead  matrix  A  =  C  is 

00 

in  diagonal  form  with  distinct  diagonal  elements.  In  this  case  we  may 
take  T  =  I  in  (2,3.2).  On  making  the  transformation  (2.3*2)  in  (2.3. l) 
we  obtain 


(2.3.5) 


We  next  choose  the  r  +  1  T^5s(k=l,2,*»*,r  +  l)  so  that 

the  matrix 

(2.3.6)  C(z)  =  T"1(z) 

has  the  form  (2.3.4).  For  this  purpose  we  observe  that  in  order  to 
solve  for  the  (r  +  l)  T  “  s  we  need  only  consider  the  coefficients  of 

IV 

-1  “2  -r-1 

z  ,  z  ,  •••,  z  in  the  equation 

(2.3.7)  T(z)  C(z)  =  A(z)  T(z)  -  z"r  . 

/  \  -1 
Expanding  (2.3.7)  an<l  equating  equal  powers  of  z  we  obtain 

l 

(2.3.8)  7  (T^k  Ck  -  Ak  T^_k)  =  °,  (i  =  +  1)  , 

k=o 


A(z)  T(z)  -  z 


-r 


d  T(z) 


dz 
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With  i  =  0,  (2.3*8)  gives  us  Cq  =  Aq.  With  Z  =  1  we  have 

(2.3.9)  Tx  Co  -  Ao  Tx  +  Cx  -  Ax  =  0. 

Now  for  arbitrary  T^,  the  diagonal  part  of  Cq  -  AqT^  =  T^  Aq  -  Aq 

is  zero.  Thus  =  diag  (A^).  The  off-diagonal  elements  of  T^  can 

be  solved  for  from  {T-}..  ({A  )..  -  {A  )..)  =  [A  }  under  the 

1'ij  oJjj  v  oJii'  1  ij 

assumption  that  the  diagonal  elements  of  A^  are  distinct.  We  set 
the  undeterminable  diagonal  elements  of  T^  to  zero. 

Suppose  that  we  have  solved  for  C  ,  C. ,  C .  and  the 

O  1  H/ 

off-diagonal  elements  of  T^,  ••*,  T^  and  that  we  have  set  the  diagonal 
elements  of  T^,  •  •*,  T^  equal  to  zero.  Then  from  (2.3*8)  we  have 

Z 

(2.3*10)  TJ+1  Aq  -  Aq  Ti+1  +  X  ( Vl_k  Ck  -  Ak  Ti+1_k)  +  CM  -  A^+1  =  0. 

k=l 


Again,  the  diagonal  part  of  T.  ,  A  -  A  T„  ,  is  zero.  Thus 

or  i+1  o  o  i+1 


the  diagonal  part  of 


Z 


(2.3*11) 


I  (Ti+l-k  Ck  "  Ak  Ti+l-k^  +  Ci+ 1  "  AZ+ 1 


k=l 


is  also  zero.  This  enables  us  to  determine  C.  ..  The  off-diagonal 

Z+ 1 


elements  of  T„  .  are  determinable  from 

i+1 

Z 

(2.3.12)  {Vlly  (f.}u  -  {*.}J  *{I 


T„  ,  .  C,  -  A,  T„  , 
Z+ 1-k  k  k  i+1 


-k)\j = 


0, 


(i  /  j) 


assuming  the  diagonal  elements  of  A  to  be  distinct,,  The  diagonal 


elements  of  T  are  again  not  determinable  and  we  set  these  equal 

1 


to  zero.  This  proves  by  induction  that  each  and 


-DIP 

—  W  9  -i-  3>c_)ooo 


can  be  determined  as  described  in  the  theorem. 


-  r  - 1  -  k 


On  again  expanding  (2.3«7)  and  equating  coefficients  of 
k  >  1  we  obtain 


min(k,r+l) 


B.  =  A 


k+r+2 


(T,  B 


l  k-£ 


4=1 


A 


k+r+2-4 


r+1 

I 

i=k+l 


Gk+4+2-i 


A.  .  T  )  +  k  T 

k+r+2 -4  l 


k 


for  all  k  >  0  if  we  define  T,  =  0  for  k  >  r  +  1.  Since  the  elements  o 
=  k 

_  i 

A(z),  T(z)j  T  (z)  are  holomorphic  and  bounded  for  all  jzj  sufficiently 

large.,  z  eod,  the  elements  of  Bo(z)  are  similarly  holomorphic  and 

s 

bounded.  This  completes  the  proof  of  Theorem  2.3.1. 


A  desirable  consequence  of  the  previous  theorem  is  that  it 
simplifies  considerably  the  construction  of  a  formal  solution  (e.g.  the 
proof  to  be  given  below  is  much  shorter  than  that  in  Coddington  and 
Levinson,  [4]  Chapter  V). 

Theorem  2.3.2  The  system  (2.3«3)  where  C(z)  has  the  canonical 
form  (2.3°^)  possesses  a  formal  independent  series  solution  matrix  of 


the  form 


' 


.1  owe  o  f  vrs-  n  . 


-  Jk  - 


(2.3.13) 


<\>  <V. 


W(z)  =  U(z)  exp  Q(z) 


where  Q(z)  Is  diagonal  and  has  the  form 


V*  Qk 

(2.3.14)  Q(z)  =  )  r  +™ "V 


r  _  r+l-k 
C,  z 


k=o 


r+l-k  V  k  z 

z  +  Q_.  z  =  J 

k=o 


r+1 


-7-  +  C  .  z 
k  r+1 


and  where 


(2.3.15) 


00 

v\_  -k 

z 


U(z)  =  I  +  /^Uk 
k=l 


Proof.  Substituting  (2.3. 13)  into  (2.3*3)>  cancelling  exp  Q(z), 
and  equating  coefficients  of  equal  powers  of  z  \  we  find,  that  if 
we  define  =  0  for  j  >  r  +  1,  U  =  0  for  j  <  0,  then 


(2.3.16) 


min(k,r+l) 

l 

j=° 


.  ,  Q.  -  C.  U  .)  -  (k  -  r  -  1)  U 

k-j  J  k- j 7 


k-r-2 


k-r-1 


I  Bj  \-r-J-2  =  0  (k  =  0,1.2,- -) 


J=° 


We  shall  prove  that  on  choosing  U  =  I  and  assuming  Q(z)  diagonal , 


o 


equation  (2.3.16)  thereafter  uniquely  determines  each  and  Q^. 


For,  with  k  =  0  we  have 


(2.3.17)  Qq  “  Cq  =  0  which  implies  Qq  =  Cq. 


With  k  =  1, 


U.  C  -  C  U,  +  Q.  -  C. 
1  o  oil  1 


(2.3.18) 


=  0. 
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Now,  if  Uj  isai  arbitrary  n  x  n  matrix,  the  diagonal  part  of 

is  zero.  By  assumption  and  are  diagonal  and 

hence  Q.  -  C  .  Since  the  diagonal  elements  of  C  are  distinct  the 

JL  JL  O 

off-diagonal  elements  of  given  by  the  solution  of 

{U.}  ({C  }  -(C)  )  =  0  are  zero.  The  diagonal  elements  of  U 

ij  °  jj  °  ii  1 

are  at  this  point  not  determinable  and  all  we  can  say  here  is  that 

UL  =  where  is  an  arbitrary  diagonal  matrix. 


Now  suppose  that  for  1  <  k  <  r  we  have  =  C^s 
where  is  an  arbitrary  diagonal  matrix.  Then  solving  for  ^  and 

Qk+i  we  have  from  (2.3*16), 

k 

(2.3.19)  U„  C  -  C  U.  .+)  (D.  .  .  C.  -  C.  D.  .  .)  +  Q.  -  -  C.  .  =  0 

7  k+1  o  o  k+1  ’  /_j  N  k+l-j  j  j  k+ 1  - j 7  k+1  k+1 


or,  since  diagonal  matrices  commute. 


(2.3.20) 


U,  -  C  -  C  U  1  +  Q.  “  A  =  0 
k+1  o  o  k+1  k+1  k+1 


and  on  comparing  equation  ( 2.3.20)  with  (2.3. 18)  it  follows  by  exactly 
similar  arguments  as  those  used  for  (2.3. 1.8)  that 

Uk+i  =  where  ^  is  an  arbitrary  diagonal  matrix.  Thus  we  have 

proved  by  induction  that  for  any  integer  k  in  0  <  k  <  r  +  l  j 
=  l>k  where  are  yet  undetermined  diagonal  matrices. 


With  k  =  r  +  2  we  have  from  (2.3,16)  and  the  above  results. 


r+1 

(2.3*21)  U  .  C  -  C  U  +  /  (D  0  ,  C  -  C.  D  . )  -  B  -  D  = 

7  r+2  o  o  r+2  /J  v  r+2- 3  j  j  r+2-j7  o  1 

j=l 


0 


0  -  tQ  ’  o8  -  (t-s«a  1*  -  I3 &  [  ♦  s*,'  «3  -  3  s+,°  <ls-t-s> 
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y 


or,  since  the  sum  )  on  the  left  is  a  sum  of  commuting  diagonal  matrices. 


(2.3.22) 


U  0  C  -  C  U  0  -  B  -  D  =  0. 
r+2  o  o  r+2  o  1 


The  diagonal  part  ofU^C-CU^is  zero;  hence 
°  r+2  o  o  r+2 

=  -  diag  Bq.  The  off-diagonal  elements  of  U  ^  are  §*ven 


by  the  solution  of 


Kl .  ({co}. .  -  {».}. 

1  ij  jj  ^  J  ii' 


B»L-  The 

L  iJ 


diagonal  part,  D  ^  of  U  cannot  be  determined  at  this  point. 

r+2  r+2 


Now  suppose  that  the  matrices  U  ,  U  ,  • - • ,  U  as  well  as 

12  Xr 

the  off-diagonal  elements  of  the  matrices  U.  , ,  U.  .  • • • s  U  *  ,  have 
b  Z+i  Z+2  r+i+1 

been  uniquely  determined,  while  the  diagonal  matrices 

(Z  +  1  <  k  <  r  +  Z  +  l)  consisting  of  the  diagonals  of 
U^(i+l<k<r+i+l)  are  still  undetermined.  Then,  from  (2.3. 16) 
we  have 


(2.3.23)  U  ,  0  C  -  C  U  tfo  + 

r+^+2  o  o  r+i+2 


r+1 

>1 


u  „  .  c.  -  c.  u  „  _  .) 
r+2+i-j  j  j  r+i+2-j ' 


Z+l 


- 1 V. 

j=l 


^1-j  -  <* +  V  Vi  =  °- 


On  writing  +  D^,  j  >  Z  we  find,  using  the  fact  that  diagonal 

matrices  commute. 


,  3d  c  oq  (  r  J a  brx  ?  Joan  a 


9  t  3  av$e  '  s.’tnrr  -  y  *  P  n*i  n*»d 
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r+1 

(2.3.24)  ur+^+2  CQ  "  c0  ^r+i+2  +  Y  ^r+i+2-j  Cj  ’  °j  ^r+i+2-j^ 

j=l 


i+1 

v 

)  B.  U.  . 
/_  j-1  i+l-j 

j=l 


(i  +  1)(P^+1  +  Vl^ 


and  by  arguments  similar  to  those  above,  U  ,  _  as  well  as  D.  _  is 

r+i+2  .0+1 

uniquely  determined  by  (2.3*24). 


Thus  each  U.,  j  >  1  is  uniquely  determined.  This  completes 


the  proof  of  Theorem  2.3*2. 


Clearly  the  above  procedure  can  be  easily  modified  to  compute 
a  particular  formal  series  solution  vector  independently  of  the  others; 
the  j ' th  vector  series  solution  being 


(2.3.25) 


~  ~  <i,(z) 

«j(z)  =  Uj(z)  e  J 


<\y 

where  XJ^(z)  t^ie  J 5 vector  in  U(z)  and  ( z ) 

later  convenience  we  set 


{Q( z) 3 j j •  For 


(2.3*26) 


qj(z) 


L 

•Y 


r  +  1  -  k 


r+l-k 
z  + 


q.  log  z. 

J»r+1 


k=o 


The  representation  (2.3.25)  illustrates  that  there  is  a  one-to- 

one  correspondence  between  the  eigenvalues  of  A  and  the  formal  vector 

o 


solutions . 


t  vo<^  .  J  or  “te.  f  *  3«  yd  bn  r 


-  38  - 


2.4  The  General  Case 


The  following  theorem  due  to  H.L.  Turrittin  [32]  concerns  trans¬ 
formation  of  the  system  (2.4.1)  with  eigenvalues  of  not  necessarily 

distinct  to  canonical  form.  By  an  exponential  transformation  in  the  state¬ 
ment  of  this  theorem  we  mean  a  transformation  of  the  form 


(2.4.1) 


^  k 

f\Z 

VI  -  e  IX 


where  W  and  X  are  matrices,  A  is  a  constant  and  k  is  an  integer.  The 
shearing  transformation  is  defined  in  [32]  page  39*  We  shall  give  an  example 
of  a  shearing  transformation  for  the  n  =  2  case  in  Appendix  A. 

Theorem  2.4.1  A  given  system  of  differential  equations  of  the  form 


(2.4.2) 


dX 

dt 


A(t)  X 


where  h  is  a  non-negative  integer ,  where  for  each  arbitrary  positive 


integer  s 


s  - 1 


00 


(2.4.3)  A(t)  =  /  Ak  t“*  +  Ag  (t)  t 

k=o 


-k 


-s  5  " 
k=o 


and  where  the  A^ ’ s  are  constant  n  x  n  matrices  and  the  elements  of 
A  (t)  (s  =  1,2,***)  are  holomorphic  and  bounded  in  some  domain  rAl .  extending 

S  1  tl 

to  infinity,  can  be  transformed  to  the  canonical  system 


(2.4.4) 


=  zT  C(z)  W 
dz  ' 


od.1  Jo  noi  JMrtvt&a,  13  >,  «c  r  ir.noartJ  «2ri1  io  Inara 


' 
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where 


r  _ ,  \  r 
z  L  ^ z)  =  z 


(2.4.5) 


z  Cjz)  = 


r+1 


k=o 


1 

1 

c1(z) 

Bll(z)  B12(z)  •”  Bli^z) 

c2(z) 

9 

9 

-2 

+  z 

B2i^ z)  B22(z)  ”*  B2i^Z^ 

•  •  • 

9 

Vz> 

Bn(z)  Vz)  Vz) 

r-k  -1 

c..  z  I  +  z 

ik  fi. 

l 

J 

*i 

s  i,  +  \i2  +  •  •  •  +  [x£  =  n 


B 


ij 


i  =  lj  2,  * 
s-1 

i  (z)  ■  X  Biik  z'k  +  Biis  <*> 


-s 


k=o 


hi  £  finite  sequence  of  normalizing,,  exponential ,  zero  inducing,  root 
equalizing,  and  shearing  transformations .  In  (2,4.4)  r  is.  an  integer. 


and  z  =  t F  for  some  positive  ii 


The  elements  of  the  matrices 


B^(z)  (i,  j  =  1,2, 0  ° ° ,  i)  are  holomorphic  and  bounded  for  z  in  the  domain 

0^9  which  coincides  with  the  domain  under  the  trans format ion  t ^  =  z, 

for  all  |z|  sufficiently  large.  If  r  =  -2  the  matrices  C^(z)  (i  =  1,2, • • • , 

are  identically  zero.  If  r  >  -2,  I  are  u.  x  u.  unit  matrices,  J 
-  - — - L  -  —  1 1±  - -  i  i - (ii 

are  fj.„  x  square  matrices  with  zeros  or  ones  or  a  mixture  of  zeros  and 

ones  on  the  first  subdiagonal  while  all  other  elements  of  are  zero. 

No  two  of  the  matrices  C^(z)  (i  =  1,2, • • • ,i)  in  (2.4.5)  are  identical. 

In  particular  if  i  j  and  c^  =  for  k  =  0,l,**°,r  then  corresponding 

to  any  particular  fixed  j  (j  =  1,2, • • * ,i) ,  we  have  0  <  Re  c^  <  1, 

-1  £  Re  cj,r+l  <  0  3Sd  ci>r+1  -  cj>r+1  /  1  (i  =  !  i  /  j). 


fii  amofa 
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The  above  statement  of  Turrittin' s  theorem  is  in  some  aspects 

different  from  that  given  by  Turrittin  in  [32].  In  the  first  place  Turrittin 

proves  this  theorem  for  the  case  when  A(t)  has  a  convergent  power  series 

expansion  about  t  =  oo  (actually  in  Turritin's  case  A(t)  converges  about 

t  =  0  but  this  difference  is  immaterial).  Secondly,  Turrittin' s  original 

version  of  the  last  statement  above  is?  Ijl  particular  if  i  jl  j  and 

c„,  =  c for  k  =  0.1.....  r  then  c.  ,  -  c.  ,  is  not  only  not  zero 

lk  jk  -  *  *  *  -  i,r+l  j,r+l - - 

but  is  also  not  an  integer . 


We  have  already  seen  that  the  normalizing,  zero  inducing  and  root 

equalizing  transformations  preserve  the  form  of  the  expansion  of  A(t) 

as  well  as  holomorphy  and  boundedness  of  the  elements  of  A(t)  for  |t| 

sufficiently  large,  t  e  In  addition  the  shearing  transformation 

_  1c 

shifts  the  coefficients  of  t  in  a  manner  similar  to  the  root  equalizing 
transformation  and  hence  it  also  preserves  these  properties. 

In  order  to  arrive  at  our  version  of  the  last  statement  of  Turrittin* s 
theorem  we  start  by  assuming  Turrittin 's  version  and  proceed  by  a  method  similar 
to  that  given  on  pages  25  and  26  of  this  chapter.  We  emphasize  that  the 
procedure  which  we  shall  now  describe  is  necessary  only  if  the  last  condition 
in  our  statement  of  Turrittin' s  theorem  is  not  satisfied. 


Consider  the  sub-set  of  integers  (i)  of  the  set  [i  =  1,2,..., 
such  that  with  i  €  [i],  j  €  [i]  and  i  ^  j,  c^  "  c  ^  =  0,  k  =  0,1,.  ..,r. 
Draw  vertical  lines  through  the  points  A  =  k  (k  an  integer)  in  the  complex 
A  plane.  If  c.  ,  lies  in  the  strip  kt  <  Re  A  <  k,  +  1  we  associate 
the  integer  k^  with  c^  r+-p  Without  loss  of  generality  we  may  assume 


(  aJ:  m  ztsh  H  i  (  «  3 


o  *vod6  3  ns  ?3r>:  3i  •  !  W1  io  riolBisiv 


>  i2<!  ££i£  0.  SJii 


i  ^  nl  -  i  rn  >-»rfi  4  f3  n.>^2  ^  »'  o  v 


och  Ho  i  ol  -,/3  svwssiq  amxiJaijJioHami  .nisi  I  .tups 


o 

•  ooi  -•  ■  .‘V  r43  :  :f:  >;•  tf  atis  cj  Ji  J  s»v,  Id! 

.  :‘ (jr.rK  ,i  [J  *1c  ,  >>r  £  2  ftagp.q  i  o  navig  3  >  rt3  oJ 


ino  \fifcar.30dn  i  adiToejb  won  liftrft  aw  r‘.;if fw  aiubaooiq 


35nioq  srfi  (Xyuot  M  »»■•  XX  wind 


that  {i}  =  {1,2, • • • 


and  that  k.  <  k_,  < 
1  =  2  = 


<  k  .  Let  k 
co 


We  shall  next  need  k  +  1  (n  x  n)  matrices  T  (k  =  1,2, • • • ,k+l) 
with  each  being  zero  everywhere  except  in  its  leading  u  x  u  block 

where  we  shall  define  it  presently.  The  sequence  of  zero  inducing  trans¬ 
formations  (2.2.7)  with  to  be  defined  here  and  with  k  (in  (2.2.7)) 

replaced  by  k  (defined  here)  +  1,  is  now  applied  to  (2.4.4).  Thus  for 
example  if 

k-l  rr„  - 1  „  -r- 1  „  -r-2  „  -k-r-1 

dz  ^  o  1  r+1  o  k  k-l 


where  the  leading  w  x  uj  blocks  of  B  ,  B  ,  *  • 0 ,  B  have  already  been 

O  J_  K  ■  JL 

filled  with  zeros  then 


dZ 

IT  =  zr[co  +  ci  z  +••• 


-r-1  -r-2  *  -k-r-1 

+  Cr+1  z  +  Bo  z  +-"]z, 


where 


B  =  B  -  T  C  +  C  .  T.  +  k  T, 
k  k  k  r+l  r+1  k  k 


This  equation  is  similar  to  (2.2,11);  hence  we  can  clearly  choose  T^  so 

* 

that  the  leading  w  x  u  block  of  B^  is  zero.  It  is  thus  clear  that  each 
leading  w  x  U)  block  of  B^  (k  =  0,l,*»«,Kfl)  can  in  this  manner  be  filled 
with  zeros. 


The  root  equalizing  transformation  (2.2.12)  with  i  =  w,  with 
k^  (i  =  1,2,***,W  ;  i  ^  j)  as  defined  here  (above),  but  with  k^  (in  2.2.12) 
replaced  by  k^  (as  defined  here)  +  1,  is  now  applied  to  (*)  (above) 


- 
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changing  ci  r+1  (i  =  1,2,  •••,(*))  to 
of  the  last  statement  of  Turrittin’s 


the  form  as  described  in  our  version 
theorem. 


For  details  of  proof  of  Theorem  2.4.1  we  refer  the  reader  to 
the  paper  cited  above.  We  emphasize  however  that  a  constructive 
(transparent' [47]  page  7)  proof  of  this  extremely  complicated  theorem 
would  be  very  desirable. 


The  following  theorem  is  also  proved  by  Turrittin  [32].  The 
main  purpose  of  presenting  a  proof  here  is  to  obtain  a  recurrence  relation 
for  the  coefficients  of  the  solution  matrix  and  to  furnish  an  explicit 
algorithm  for  computing  these  coefficients. 


Theorem  2.4.2  The  system  of  differential  equations  in  canonical  form 
(2.4.4)-(2.4.5)  possesses  a  formal  independent  series  solution  matrix  of 


the  form 


(2.4.6) 

where 


W(z)  =  U(z)  e 


'J(s) 


00 


00 


”(*)  =  iv*»  =  Z  \  z'k  =  Z  [v]  z 


-k 


k=o 


k=o 


A/ 


U  (z)  is  |~i .  x  |i . ,  U.  =  5  I 
lj  1  J  ijo  ij  u. 


(2.4.7) 


r 

Q(z)  =  [(i^  qt(z)  +  log  z)&ij^  =  y 


xk  r+l-k 

TTTk  z  +  Qr+1  108  Z 


k=l 


‘i(z>  =  z 


qikz 


r+l-k 


r+l-k 


+  q 


.  l°g  z  (l»j  =  1,2,  •••,!). 

i,r+I 


k=o 


. 
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Proof  of  Theorem  2.4.2  Let  us  put 


(2.4.8) 


r+1 

zr[C.(2)Bij]  =  Y  Ck  zr+1"k  =  *  c(z) 

k=o 

s-1 

j  ( z)  ]  =  B(z)  =  y1  Bk  z  k  +  Bg(z)  z"s 

k=o 


Hence  if  we  substitute  (2.4.6)  into  (2.4.4),  use  the  fact  that  Q(z)  and 
Q'(z)  commute  we  obtain 


(2.4.9) 


U* (z)  +  U(z)  Q* (z)  2 


C(z)  + 


U(z) 


powers  of 
(2.4.10) 


Expanding  (2.4.9)*  using  (2,4.7)  and  (2.4.8)  and  equating  equal 
z  we  obtain 

min(p,r+l)  p-r-2 


(ck  Vu 


-  U  Q,  )  +  (n-r-l)O  ,  + 
p-k  k7  ^  7  p-r-1 


B  U 

k  p-r-2-k 


k=o 


k=o 


0 


(p  =  0,1,2s040)  where,  in  order  for  this  equation  to  be  defined  for  each 
non -negative  integer  p  we  put  =  0  if  k  <  0,  Ck=Qk=0  if 
k  >  r  +  1. 


The  proof  of  Theorem  2.4.2  will  be  complete  if  we  prove  the 
following  lemma. 


Lemma  2.4.1  Under  the  assumption  that  Q(z)  has  the  form  given  in 
(2.4.7)  and  UQ  =  I*  equation  (2.4.10)  thereafter  uniquely  determines 
Q  (p  =  0,1,°°°, r+1)  and  U  (p  =  l,2,oo«). 

r1  H’ 


0  0  9 


(8.4.S) 


l 


. 

let,  )  $r;t  -JSI/p.9  >'fj&  (8.-f.S)  b  >r  (  ,  :»  t(9,4.S )  fcoM>«»qpi3 

. 


. 


»v.  OJ  .  j'L  [  1  70  '  ,  X»J>  v»*»tSfIt0  *  m) 


ninroiloi 


^  aj  *  [r  -  ti  j  j  (y.4.s) 


•  I  '  14)  '<  h£H  (  ,J,0  «  m)  0 
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Proof  of  Lemma  2.4,1  The  proof  of  the  lemma  is  by  induction. 

With  p  =  0,  (2.4.10)  gives 

(2.4.11)  C  -  Q  =  0 

which  has  the  unique  solution  Q  =  C  .  With  p  =  1,  (2.4.10)  becomes 

00 

(2.4.12)  C„  -  Q,  +  C  IL  -  U.  C  =0. 

'  Iloilo 


Now  if  is  an  arbitrary  n  x  n  matrix  we  know  from  the 

definition  of  Cq  that  the  (i,i)'th  (i  =  l,2,*°*,i)  blocks  of 

Co  -  Ux  Co  consisting  of  p_^  x  p^  matrices  are  identically  zero. 

Hence  these  blocks  of  C,  -  Q,  must  also  be  zero,  and  therefore  Q,  =  C, . 

Furthermore  it  follows  that  all  (i,j)9th  blocks  of  for  which 

c.  -  c  .  J  0  are  uniquely  determined  and  are  zero.  All  (i,j)'th 
10  jo  '  / 

blocks  of  U,  for  which  c,  =  c  „  cannot  be  solved  for  at  this  point. 

1  10  jo 

In  particular  it  is  not  possible  to  solve  for  U  ^  (i  =  1,2, • * 0 ,i) • 


Proceeding  in  this  manner,  let  us  suppose  that  we  have  just 
completed  solving  (2.4.10)  for  p  =  k  (0  <  k  <  r);  that  we  have  found 


Q  =  C  (p  *  0,1 

p  p 


000 

S>  S> 


k)  and  that  we  have  determined  all  blocks  of  the 


matrices  U  (p  =  0, 1, 0 0 • ,k)  as  far  as  possible.  More  precisely,  we  make 
M> 

the  following  assumptions  for  each  (i,j)9th  block  of  U  (i,j  =  l,2,°**,i 

M' 

p  =  0, 1, 0  •  •  ,k)  .  If  c  -  c  =  0  (p  =  0,1, 0  • 0  ,k~l)  then  assume  U  , 

ill  JM 


Ufj2,  0 0 # »  are  arbitrary  so  far  undetermined  p^  x  p^  matrices.  On 


the  other  hand,  if  p  is  the  smallest  non-negative  integer  in  0  <  p  <  k 


(JJtV 

•  :  ■  '  •  .  •  ,'i  v  .  '  j  »i.  :  rialrfv 

* 


♦ 


* 

t- 


for  which  c  -  c  j/  0  we  assume  U.  .  =  IT  ,,  =  •••  =  U.  .  -  =0 

ip  jp  r  ijo  ljl  ij>k-p 

While  Uij,k+l-p*  Uij,k+2-p>  -"-“ijk  are  arbitrary  so  far 
undetermined  p^,  x  p^,  matrices.  Under  these  assumptions  let  us  proceed 
to  p  =  k  +  1. 

For  p  =  k  +  1  equation  (2.4.10)  becomes 


(2.4.13)  C  U  -  U,  C  + 

'  o  k+1  k+1  o 


+  ck  U1  -  uick +  Ck+1  -  Vi  =  °- 


If  u. ,  U 


1’  '*'2* 


,  Uk+1  are  arbitrary  n  x  n  matrices  the  (i,i)'th 


(i  =  l,2,,*’,i)  blocks  of  C  U.  .  -  U,  _  C  are  zero  (p  =  0,1, •••,k) 

v  *  *  *  '  p  k+l-p  k+l-p  p  7 

Hence  these  blocks  of  C,  ,  -  Q,  -  must  also  be  zero,  and  therefore 

k+1  k+1 

Qk+i  =  Let  us  incorporate  this  result  in  (2.4.13)  and  let  us 

consider  the  (i,j)“th  block  of  (2.4.13)  (i,j  =  1,2, • • • ,-0 .  Here  we  have 


(2.4.110  (c.o  -  cj0)aij(k+1  +  ■  cjl^Uij,k  + 


+  (cik  -  cjk'iUijl 


If  c.  -c.  =0,  p  =  0,1,  •••  sk¬ 

ip  JP 


then  all  we  can  say  is  that 


U..-,,  U.  .  are  so  far  undetermined  arbitrary  matrices.  On 
ljl  ij2  ij,k+l  J 

the  other  hand,  if  p  is  the  smallest  non-negative  integer  in  0  <  p  <  k 
for  which  c^  -  c  ^  0,  then  by  the  above  induction  hypothesis  (2.4.14) 
becomes 


(2.4.15) 


(c  -  c .  )  U . .  ,  .  =0  (p.  x  p.) 

ip  JP  ij, k+l-p  Vhi  "V 


and  so  U..  ,  ,  =0,  under  the  assumption  that  U. .  ,  _  , . 

ij, k+l-p  r  ij ,k+2-p  ij ,k+l 

are  arbitrary  p^  x  p^  matrices. 


■ 


is  uniquely  determined 
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It  is  thus  clear  at  this  point  that  Q 


for  each  p  =  0, 19 • • • »r+l,  and  we  may  write  the  (i,j)'th  block  of 
(2.4.10)  in  the  form 
r+1 

(2.4.16)  y  (qj,  -  q 


k=o 


4-  -  ,  +  j  u.  .  _  -  U.,  -J  +  (p-r-l)U..  . 

ij  |jt±  ij jM--r-l  ij,p-r-l  p  ij,|i-r-l 


|i-r-2  i 

*  I  E 

k=o  s=l 


B  U  oi=0 

isk  sj,p-r-2-k 


For  p  as  r  +  2S  (2.4.16)  becomes 


r+1 

(2.4.17)  y  (qik  -  ^ jic)Ui 
k=o 


..  /  ~  .  +  J  U».-s  ~  U .  .  1 J  +  U .  . ,  +  B  .  .  =  0 

jk  ij  sr+2-k  ljl  ijl  p^  ljl  1  jo 


if  q.  "  q.  =0  (p  =  0,1 

ip  jp 
matrix 


0  0  9 

9  9 


r)  then  for  U. .  _  an  arbitrary  p.  x  p, 

ij  #r+2  1  j 


(2.4.18)  (1  +  q.  .  -  q.  JIT,,  +  J  U.  41  -  U...J  +  B.,  =  0, 

'  i9r+l  j  j.r+1'  ijl  ni  ijl  ijl  p ljo 

where  under  the  assumption  of  Theorem  2.4.1  q.  .  -  q,  1  ^-1  if  i  ^  j . 

In  either  case,  whether  i  =  j  or  i  /  j  equation  (2.4.18)  is  of  the 

same  form  as  equation  (2.2.11)  and  hence  it  uniquely  determines  U 

On  the  other  hand,  suppose  p  is  the  smallest  non-negative  integer  in 

0  <  p  <  r  for  which  q.  -  q.  /  0.  Then  from  the  above  analysis  we  know 

that  U„ .  =  U ...=•••  *  U  .  =0  (p.  x  p„).  Thus  if 

ijo  ijl  ij,r+l-p  Vfi  *j' 

Uij Sr+3“P*  Uij ,r+2  are  arbitrary  V±  x  p,  matrices 


(2.4.19) 


^qiP  qjP^Uij»r+2-p  +  Bij°  ° 


-  ,  iic  :  '■>  »oi 


xl  j3sur< 


»  on- n  t  -m!  i  r  t  q  '  a  ,  na 4  i*riJo  1  nO 

D  ,13 life  TO 5 


which  has  up  to  this 


and  this  equation  uniquely  determines  U. . 

iJ,r+2-p 

point  been  defined  only  to  be  an  arbitrary  p_^  x  p^  matrix.  Thus  we  see 
that  the  p_„  x  p_  matrix  is  at  this  point  uniquely  determined  for 

each  pair  of  integers  (i,j  :  i,j  =  l,2,»**,i)  (unless  of  course,  only  the 
j ' th  column  of  blocks  is  desired,  in  which  case  it  is  unnecessary  to 
compute  any  other  column  of  blocks,  since  a  particular  column  of  blocks 
can  be  computed  independently  of  any  other),  and  hence  has  been 
completely  and  uniquely  determined. 


Proceeding  in  this  manner,  let  us  suppose  that  we  have  just 

completed  solving  (2.4.10)  for  p  =  r  +  2  +  k,  k  >  0;  that  we  have 

uniquely  determined  U^,  •••,  ^  and  that  we  have  determined  all  blocks 

of  the  other  matrices  ^+s»  s>  k  as  far  as  possible.  More  precisely, 

we  make  the  following  assumptions  for  each  (i,j)'th  U  ^  of 

IL  ,  0  <  p  <  r  +  2  +  k.  If  q.  -  q.  =0,  (p  =  0,1, •••,r)  then 

assume  U..,,  •••,  U. .  .  ,  are  arbitrary  so  for  undetermined  p.  x  p, 

ijl  ij  ,r+2+k  17  pi 

matrices.  On  the  other  hand,  if  p  is  the  smallest  non-negative  integer 


in  0  <  p  <  r  for  which  q.  -  q.  4  0  then  assume  U, • 0 • ,  U 
=  -  ip  JP  r  ijl 


have  been  uniquely  determined  while  U. .  ,  ,  •••,  U..  _  , 

ij ,r+3+k-p  lj ,r+2+k 


ij ,r+2+k-p 

are 


arbitrary  so  for  undetermined  p^  x  p^  matrices.  Under  these  assumptions 
we  proceed  to  p  =  r  +  3  +  k. 


With  p  =  r  +  3  +  k  (2.4.16)  is 


r+1 

(s.t.20)  £  (q1,-qJ,)«1J. 

s=o 


-z  1  +>J  U..o1  -  U.. 

r+3+k-s  p_^  ij,2+k  ij,2+k  p 


+  (k  +  2)Uij,2+k  + 


k+1  Z 


S-O  V— 1 


S.  U  .  =  0 

ivs  vj,k+l-s 


,0.-'  :  • 


*•  .cure** 
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If  -  q.  =0  ( s  =  0,l,*#,,r),  then  under  the  assumption  that 

is  j  s 

U.  .  ,  is  an  arbitrary  p  x  p.  matrix,  by  the  above  induction  hypothesis, 

ij>r+5+*:C  i  j 

and  the  assumptions  of  Theorem  2.4.1,  equation  (2.4.20)  is  an  equation  in 

U.,  vo  °f  exactly  the  same  type  as  (2.2.11)  was  in  T...  ,  and  thus 
ij,k+2  ljk 

(2.4.20)  uniquely  determines  U  ,  .  On  the  other  hand,  suppose  p  is 

i j , k+2 

the  smallest  non-negative  integer  in  0  <  p  <  r  for  which  q^  -  q^  ^  0. 

Then  assuming  >r+1)+k.p,  Uij>r+3+k  to  be  arbitrary  n.  x  ^ 

matrices,  equation  (2.4.20)  may  be  written 


(2.4.21)  (q.  -  q.  )U. .  _  . 

ip  jp'  ij,r+3+k-p 


r+1 

7  (qis  -  qjs)Uij, 


r+3+k-s 


s=p+l 


-  (k  +  2)U 


i j  »k+2 


k+1  4 


s=o  v=l 


B.  U  4  1  ! 
ivs  vj,k+l-s 


By  the  induction  hypothesis,  all  matrices  on  the  right  of  (2.4.21)  are 

known,  and  hence  U.  .  ,  .  is  uniquely  determined. 

ij  »r+3+k~P 

Thus  each  U  (p  =  1,2,*»»)  is  uniquely  determined  by  equation 

M* 

(2.4.10).  This  completes  the  proof  of  Lemma  2.4,1  and  also  of  Theorem  2.4.2. 


By  a  discussion  analogous  to  that  leading  to  equation  (2.2.19) 

we  obtain 


(2.4.22) 


k=o  s=v 


0  J  >u-  i  d;  . ,  n  x.  i. 


1  :  ' 


i  >  q  >  0  n.  avian  '>n-no«  3a  >Ii jwj  i  »rf3 


t1**!**  ^,UU  •  “  'q-^.U® 


nsi.-'t'B  d  /.  it  <0S»4.S)  aolJaup*  «e*3*Tj£cx 


v  .-  .2;  aolJttrpo  03  gal£&al  ja/fi  oj  auqgclr./rn  noia  ooelb  a  ^ 


(2S.4.S) 


for  the  v'th  vector  solution  in  the  j  *  th  column  of  blocks .  In  equation 

(2.4.22)  we  have  assumed  that  J  has  a  one  in  its  v’th,  v  +  l'th,*'*, 

s 

[i  -  l'th  columns,  but  not  in  its  jj.  1  th ,  where  p  <  p .  and  U  is  the 

-  j  jk 

s'th  (s  =  l,2,***,p.)  column  vector  in  the  j'th  (j  =  l,2,*#*,i)  column 
of  blocks  of  the  matrix  (k  =  1,2,***;  Theorem  2.4.2).  Thus  to  obtain 

a  formal  solution  it  is  again  only  necessary  to  compute  the  p'th, 


p  -  l'th,  •••,  v'th  column  of  at  each  stage. 


CHAPTER  III 


ERROR  BOUNDS  FOR  FORMAL  PARTIAL  SUM  APPROXIMATIONS 

TO  ACTUAL  SOLUTIONS 


It  is  to  be  expected  that  the  differential  equations  of  the 
previous  chapter  possess  actual  solutions  for  which  the  formal  solutions 
are  asymptotic  expansions  (for  proof  of  this  statement  see  [8]). 

In  this  chapter  we  generalize  Olvers  (see  e.g.  [35])  procedure 
for  the  second  order  equation  in  that  we  solve  the  differential  equation 
for  a  partial  sum  of  a  formal  series  solution.  Thus  in  setting  up  integral 
equations  for  the  error  vector  we  simultaneously  prove  the  existence  of 
actual  solutions  for  which  the  formal  series  solutions  obtained  in  the 
previous  chapter  are  asymptotic  expansions,  as  well  as  obtain  error  bounds 
for  partial  sum  approximations  to  actual  solutions. 

The  error  bounds  are  given  in  two  forms;  as  vector  bounds  and 
as  norm  bounds. 


3. 1  The  Case,  of  a  Regular  Point 

3.1.1  The  Differential  Equation  for  an  Approximation 

We  shall  find  an  actual  solution  vector  corresponding  to  the 
j’th  formal  solution  vector  in  (2.1. 3).  For  this  purpose  we  denote  the 
j'th  column  of  U^  (k  =  0, 1, 0 0 • )  as  defined  by  Theorem  2.1.1  by  U^. 

The  partial  sum 

m-1 

(3.1.1)  V(Z)  =  Ujk  z-k 

k=o 


satisfies  the  differential  equation 


♦  <  i  ii/Io  o:  11  .14  kfxoVK*  i  jb  X*. r *x  q  ioI 
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where  on  partially  expanding  A(z)  as  in  (2.1.2)  and  using  (3.I.I)  and 
(2.1.6)  we  find  that 


m- 1 


(5.1.3) 


R.  (z)  =  m  U.  z 
jm  jm 


-m- 


1  -  I  Vk(z)  v 


-m-2 


k=o 


Hence  if  a  vector  W.  =  W.(z)  satisfies  (2.1.1)  then  the 

J  J 


error  vector 


e.  (z)  =  W  .  ( z )  -  tp.  (z) 
jm  j  jm  ' 


satisfies  the  differential  equation 


(3.1.5) 


T~  e.  (z)  -  z  ?  A(z)  e.  (z)  = 
dz  jm  '  v  '  jm 


-  R.  (z) 

jm 


where  R^(z)  is  given  in  (3.1. 3). 


5.1.2  Matrix  and  Norm  Error  Bounds 

Let  £)$( z,£)  be  the  union  of  a  set  of  points  z  which  can  be 
joined  to  £  by  a  path  P  with  the  following  properties: 

1)  P  (with  the  exception  of  the  point  £  if  £  is  a  boundary 
point  of  0&  )  lies  wholly  in  oO; 

2)  P  consists  of  a  finite  number  of  Jordan  arcs  t  =  t(s) 
such  that  t’(s)  is  continuous  and  non-vanishing; 


3) 


(3.1.6) 


7^  (t_k)  =  f  lk  t'k_1  dt  |  (k  =  1,2,...) 


is  bounded. 


rOJosv  I01*1  » 


<«V 


.. 


5-  ,  . 


.(  i.^)  »*  *  »*J  *  8l  (s  I  w 


. 


;  $0  <Morfw  «  I 
£ 3TE  .  ?o  SJ^r/j:  s  o  wJai  n  o  H  ( S 

- 


Now  if  G.  (z)  is  a  solution  of 
jm 


(3.1.7) 


€  (z)  =  r  [t  2  A(t)  € .  (t)  -  R.  (t)]  dt 


then  g.  (z)  also  satisfies  (3.1.5). 

jm  '  s  1 


By  the  analysis  of  /)  ppfrjdi*  B  q$(z,£)  a  domain,  and  the 
solution  of  (3.I.7)  exists  and  is  a  unique  vector  of  functions  holomorphic 

J9  (*.5). 


m 


By  (3.1*3)  we  have,  with  A. .  (z)  =  {A  ( z) ) .  . ,  U. .  =  (U  ) 

J  w  '  ijs'  '  '  s  lj  ijm  nrij 

and  R..  (z)  the  i’th  element  of  R.  (z),  (i  =  l,2,*,,,n),  that 
ijm  jm'  '  v 


(3.1.8) 


£  I”  iUijml  +  ’’ijm 


where 


m-1  & 


(3.1.9) 


(m  +  1)  7. 


ijm 


sup 
t  Gtf 


X  Z  Ais,m-k  ^  Usjkl 


k-o 


S  =  1 


We  also  define  a  matrix  of  non-negative  elements  B  by 


(3.1.10) 

B  = 

sup 
t  €  g=> 

to  obtain,  on 

denoting  |V^| 

to  be 

elements  V. . 

ij 

of  the  vector 

V., 

J 

ij 


that 


is  s  08 1  »'•  (s)  n)  3 


if  ,r  ’  £  J*  !  ,  1 


.  ,\  (i.  ♦  ■■>)  (•:•■  ) 


f  *8 
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Applying  Lemmas  C.l  and  C.2  of  Appendix  C  to  (3.1.10)  we  obtain 
Theorem  3.1.1  Let  (z,£)  be  defined  as  above.  Corresponding  to 


each  formal  solution  vector  VL(z)  21.  (2.1.3)  there  exists  an  actual 
solution  vector  Wjm(z)  oj.  functions  holomorphic  in  p$L(z,£)  and  depending 
on  £  and  an  arbitrary  positive  integer  m  such  that 


m-1 


(3.1.12) 


W 


(z)  =  y  u 


Jk 

k=o 


-k  (  * 

Z  -r  G  .  (z) 

jmv 


where 


(3.1.13)  |eJn(0l  <  exp  +  7  . 


In  equation  (3. 1.12)  ijs  the  j 3  th  column  vector  of  the  matrix 

\  (k  =  0, 1, • • • )  defined  in  Theorem  2.1.1,  while  in  (3. 1. I3)  the  matrix 
B  and  the  vector  7^  are  defined  by  equations  (3.I.I0)  and  (3*1.9) 
respectively. 


Similarly s  to  obtain  a  norm  bound  for  e^m(z)  above,  let 

m-l 


(3.1.14) 


Oil.  (m  +  l)  7  =  sup  II  y  Am  k(t)U 


B”  /UP,o  T  *'  'jm_  » /_,  “m-k'*'”jk 

t  1  £  f  k=o 


where  any  one  of  the  definitions  of  norms  in  Chapter  I  may  be  used. 


Tl^orem<i^!Sals^  A  corresponding  norm  bound  €^(z)  in  Theorem  3*1.1 
is 


(3.1.15)  He  («)|  <  exp  (B7^(t'1))  {||Ujm||7^(t-m)  +  7.my^(L-m-1)) 


Jm" 


1 
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where  B  and  7 are  defined  by  equation  (3. 1.14). 

3.2  The  Case  of  a  Regular  Singular  Point 

5.2.1  The  Differential  Equation  for  an  Approximation 

With  reference  to  Theorems  2.2.1  and  2.2.2  we  shall  find  error 
bounds  for  an  actual  solution  vector  of  (2.2.3)  corresponding  to  the  v’th 
formal  solution  vector  (2.2. ly)  in  the  j  '  th  column  of  blocks  of  formal 
solutions . 


The  partial  sum 


m-1 


(5.2a)  <pjm(z)  =  ( y  Ujk  z'k)  exp  [(A.  y  +  Jj  log  z] 

k=o  J  2 


satisfies  the  differential  equation 


(3.2.2) 


d  cp. 


im 


dz 


( z )  -  z  1  B(z)  <p.  (z)  =  R .  (z)  exp  [A  .1  +  J  )  log  z]. 


jm  jm 


1  u  .  [i . 

J  J 


On  expanding  the  left  hand  side  of  (3.2.2)  we  obtain 


(3.2.3) 


R.  (z)  =  [U.  (A.I  +  J  )  -  B  U.  ]  z"1 
jm  jo  J  ly  o  jo 


m-1 


+  )  {-k  U,.  +  U,.  (A. I  +  J  )  -  B  U. 

l_j  jk  jk  j  \±  [±  '  o  j 


k=l 


j  'j 


jk  ”  /  , Uk-i 5  Z 

4=1 


-k-1 


m 


m-1 


-  (I  Bi  Vi  )  z'”'1  '  (X  Vi-s(2)ujs)  z 


-m-2 


4=1 


s=o 


•  J  (  .  ,  ,)  :•  rro  uXoa  I*  mol 


(3.S.5) 


i!o  a  ha  »X  »di  *.  ?£bafi<jx»  nO 


-  55  - 


But  by  equation  (2,2.15)  we  obtain 


(3.2.4)  R.  (z)  =  [m  U.  +  B  U.  -  U.  (A. I  +  J  )]  z 
3m  jm  o  jm  jm  j  p  p  ' 

J  J 

m-1 


-m- 1 


-d 


ujs ) 2 


-m-2 


s=o 


Thus,  if  a  block  of  p  vectors  W,(z)  satisfies  (2.2.3)  the  error 

j  J 

block  of  vectors 


(3.2.5) 


Tj.  (z)  =  [W  .  (  Z  )  -  Cp.  (z)] 
3m  J  jm 


satisfies  differential  equation 


(3.2.6)  “  q.  (z)  -  z  1B(z)qc  (z)  =  -  R,  (z)  exp  [(A. I  +  J  )  log  z] 

'  dz  'jm  '  v  '  ‘jm  jm  '  j  (i,  p  ' 


where  R^(z)  is  given  in  (3.2.4)  above. 


3.2.2  Matrix  and  Norm  Error  Bounds 


With  oQ  the  domain  in  which  B(z)  (Theorem  2.2.1)  is  holomorphic, 
let  oQ*  =  o9-  (0)  -  {00}^.  We  define  a  region  J9(z  „£)  by  the  union  of 
the  set  of  points  z  (/  £)  which  can  be  joined  to  £  by  a  path  iP  with 
the  following  properties*. 


dt  ( s) 

ds 


order 


1)  Except  for  £9  ft3  lies  wholly  in  qQ ; 

2)  (P  consists  of  a  finite  number  of  Jordan  arcs 
continuous  and  non-vanishing  for  each; 

3)  r/t  remains  bounded  for  every  pair  of  points 
£,  t,  t,  z  on  P  ; 


t  =  t(s)  with 


t ,  t  in  the 


* 


By  e.  g.  {00}  we  miUln  the  set 


00  e 


:e 


6  real. 


(*)  3  STflrfw 


•  i 
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*o  =  f  ■ 1)  fcT1 2  dtl 

J(p 

is  bounded ,  where  rj  may  be  taken  to  be  zero  if  J  =  0,  but  otherwise 
rj  is  an  arbitrary  positive  number  less  than  1. 

Set  B  .  =  A„I  +  J  .  Equation  (3.2.6)  may  be  written  in  the  form 

oj  J  Pj 

H  1  o  B  log  z 

(5-2-7^  dl  V(z)  ‘  2  Vjm(z)  =  2  "  B1(Z)V(Z)  '  Rjm(z)e 


By  our  choice  of  a  canonical  form  in  Section  2.2.  and  our  assumptions  on 

the  path  above  f'  |p  dtj  exists ,  where  p  is  any  element  of  the 

matrix  Rc  (t)  exp  (B  .  log  t}.  Hence  if  n .  (z)  is  a  solution  of 
jm  7  oj  J  'jm  7 

rz  z  B^ilog  1 

(3.2.8)  T)  (z)  =  /  exp  (B  log  “}(t  B  (t)  q  (t)  -  R  (t)e  J  )  dt 

Jin  J  ^  o  l  i  jm  jm 


where  the  integral  is  taken  along  one  of  the  paths  P  of  the  family 
described  above,  q^(z)  also  satisfies  (3.2.7). 

A. 

Let  us  put  q.  (z)  =  e.  (z)  z  and  let  us  consider  the  v'th 

r  'jmv  7  jm  7  9 

column  vector  (z)  of  e  fz).  With  reference  to  the  discussion 

jm  jm  7 

surrounding  equation  (2.2.17)  the  integral  equation  for  e^(z)  takes 
the  form 


(3- 2. 9)  <5m(z)  =  /  |(Bo  -  A.I)  log  Q}|t-2Bl(t)  ^m(t) 


S”  v 


•  .  • 


. 


1 

A  4. 


■  :rol  '  V.  ’  3  n-  s  be  1  a  •  ,nic>  k^lrv  i'>  *  Mo  a  vfl 

, 


82 


to 


1  ■■  :  (  Y* ''  1  ■  J  ‘  1  •  : '  V 


«oi  e  'iz2b  9/i  ol  ©an©*©!©*  ntoitf  , (a)  3  5o  s  »  to2s  iv  r  wloo 


so1  (v  *  .8i  *“  \  <<•  «-*'<> 
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g 

where  R  (t)  is  the  s'th  column  vector  (s=v,  v+1,  \i)  in  R.  (t). 

J  «3  ^ 


Now 


(3.2.10) 


j  /  -B  log  t  B  .  log  t' 

&  (*  °  >-  • 


-B  log  t 


=  -  e  0  ~  [m  U.  +  B  U,  -  U.  B  .]  t"®"1  e 


B  .  log  t 


jm  o  jm  jm  oj 


thus  the  first  term  of  R.  (z)  in  (3.2.4)  is  in  a  form  suitable  for 

jm  '  ' 

bounding.  The  contribution  to  the  v'th  error  vector  of  the  second  term 
In  (3.2.4)  is 


(3.2.11) 


r 


m-l 


exp  [(Bq 

i 


-V)log®]  7 


b  .  (t)us 

jp  (s-v): 


p=o  s=v 


To  obtain  a  vector  bound  we  choose  a  fixed  number  where 

Tj .  =  0  if  J  =0  (i  =  1,2, • • • ,£)  but  is  otherwise  an  arbitrary  positive 
1  Mi 

number  less  than  1,  and  set 


(3.2.12)  (m  +  1  -  r^) 


=  sup 
t  »  T  €  JP 


t"^  jexp  [  (Bq  -  Vl)log  ] 


m-l 


s- 


B 


p=o  s=v 


(t)us  (-1-og  t} 
rn-l-p  JP  (s  .  v): 


where  t  and  r  are  located  on  "Jp  in  the  order  t,  t,  z,  and  where 


r  *0) 

indicates  the  co'th  element  of  the  group  of  elements  labelled  i 
(corresponding  to  the  i'th  row  of  blocks  of  B(z))  of  the  vector  V  ^ .  We 
define  7J  such  that  {7^  =  7^. 


tc  9 


d/fT  .grtjtbaucd 


■  • •  f  "ttu  •  'it'!.  <  '  'J  ,!..  r:  do  '.*■ 


}M 

tf3  ar  si  i  d 


t,|r  -  I  ♦  ■>  (SI.S.I 


:) 


no  .»oI  »  a i  t  bns  3  sxsrfv 


lo  lool:  5  3  wox  ••/:)'  .  arf ;  03  ;;crlbno<j89xxoo ; 
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Similarly  we  fix  a  number  rj  where  r)  =  0  if  J  =0  but 
otherwise  rj  is  any  positive  number  less  than  1,  and  for  t  and  t 
in  the  order  t,  r,  z  on  we  define 


(3.2.13)  (1  -  n)  b  = 


sup 

■t.T  ey 


-n 


exp  [(BQ-A^l)log  Q]  B1(t)|- 


where  {A}^  is  the  (p,s)'th  element  of  the  (i,j)'th  block  of  A. 


By  Section  B.l  of  Appendix  B  (z,£)  is  a  domain  (actually  a 
Riemann  surface)  and  the  solution  to  (3.2.9)  exists  and  is  a  unique  vector 
of  functions  holomorphic  in$$(z,£).  Applying  the  above  inequalities  to 
(3.2.9)  we  have 

(3.2.14)  l^)l5/fl(i-’i)|t’'-2||B|£;(t)|  +  |i£  ujm  n 


s=v 


+  («+  1  -  /j  |t-m-1+^|]|dt|  . 


Applying  Lemmas  C.l  and  C.2  of  Appendix  C  to  (3.I.I3)  we  obtain 


Theorem  5.2,1  Let  j()  (z,Q  be  defined  as  above.  Corresponding  to 
the  y 1 th  formal  solution  vector  in  the  j '  th  column  of  blocks  of  formal 
solutions  as  described  by  Theorem  2.2.2.  for  the  equation  (2.2.3)"(2.2.4), 
there  exists  an  actual  solution  vector  W^m(z)  of  functions  holomorphic 
in  A  z ,  £  )  depending  on  £  and  an  arbitrary  positive  integer  m,  such 


that 


(3.2.15)  w v  U)  = 

jm 


m-1 


-J 

p=o  s=v 


S  (log  z)  V  -p  V  /  V 
U.  -  Z  r  +  €.  (z) 

jp  (s  -  v): 


a  (r  -  x) 


- 


.0,s)  .  j  ni  olfi  iiomolori  snolJ&atft  io 


- 
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where 


(3.2.16) 


s=v 


(log  t)S"V 

(s  -  v): 


+  7 


for  each  z  in  Jfy ( z ,  £ )t  The  integers  p  and  v,  and  the  vectors 

are  defined  as  in  equation  (2.2.17).  If  each  J  =0  (i=l,2,***,l0 

=  .  Mi 

then  the  numbers  r)  and  may  be  taken  to  be  zero;  if  J  -  0 

then  Tj  may  be  taken  to  be  zero ;  otherwise  r\  and  r)1  are  arbitrary 
positive  numbers  less  than  1.  The  n  x  n  matrix  B  =  B(rj)  and  the 
n  x  1  vector  7^  =  )  are  defined  through  equations  (3*2.13)  and 


(3.2.12)  respectively. 


To  obtain  a  norm  error  bound,  we  first  choose  any  of  the  definitions  of 
norms  given  in  Chapter  I.  Next,  for  t]  and  %  defined  as  for  Theorem 
3*2.1  and  t  and  r  in  the  order  £,  t,  t,  z  on  ^  ,  we  define 


(3.2.17)  (m  +  1  -  rij)  7V 

=  sup  || t  ^  exp  [(B 
t,r  €  p  0 


Vl)log 


m- 1 


p=o 


(t)U 


s  (log  t)S'V 
(s  -  v)! 


and 


(3-2.18) 


B 


sup 
t ,  T  € 


exp 


A.l)log  (f)]  B^t) 


io‘;  ,3X^.1  .1  li  )q;..!0  ni  a»v>n  .mo n 


ban 
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Theorem  5,2.2 

b£ 


A  norm  bound  for 


in  Theorem  3*2.1  jls  given 


(3.2.19) 


s=v 


where 


and  B 


respectively. 


are  positive  numbers  defined  by  (3.2. I7)  and  (3.2. 18) 


3.3  The  Irregular  Singular  Case  when  the  Eigenvalues 

of  Aq  are  Distinct 

5.3.1  The  Differential  Equation  for  an  Approximation 

We  shall  obtain  an  actual  solution  vector  W 

'V 

to  the  j'th  formal  solution  vector  W^(z)  as  defined  by  (2.3.25).  The 
procedure  simultaneously  gives  us  a  bound  on  the  difference  between  an 
actual  solution  and  the  sum  of  the  first  m  terms  of  a  formal  solution. 


j 


(z)  corresponding 


We  denote  the  j'th  column  of  (k  =  0,1,  • 
by  Theorem  2.3*2  by  U^*  l*et  tbe  partial  sum 


)  as  defined 


(3.3.1) 


m-1 

%(z)  ■  (I  V 

k=o 


-k\  qj(z> 
z  e  J 


satisfy  the  differential  equation 


(3.3.2) 


dP,m(z)  r  q.(2) 

- j -  -  z  C(z)  $ .  (z)  =  R.  (z)  e  J 

dz  *3™  jmv 


nartv 


.Yl»Vjt30»q«3l 


C.5 
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where  C(z)  is  defined  by  (2.3.4);  that  is,  we  use  (3.3*2)  to  define  the 

vector  R.  (z). 

Jin 


Expanding  (3.3. l),  using  equations (2. 3.26)  and  (2.3. 16)  we 


obtain 


(3.3.3) 


R.  (z) 

jm 


r+m+1  min(u.,r+l) 

r* — i  r  r — t 

••>.  1>,  (\  -  V  uL-k  +  (^-r-l)uj, 

(j.=o  k=o 


fi-r-1 


p-r-2 

V 

+  )  B,  U. 


k  j ,  fj.-r-2-k 


1  ,»>* 


m- 1 

V1 


/ 


k=o 


k-o 


(Vk<2>Ujk)  2 


-m-2 


where  U  =  -jU  ^  if  0  <  k  <  rn-1;  0  otherwise}*-  . 


Consider  the  i'th  element  R. ,  (z)  of  R.  (z).  If  i  =  j 

ljm  jm 

then,  using  (2.3. 16)  once  more  we  obtain 


(3*3*4)  R.  .  (z)  = 
J  jm 


m- 1 

y  (B  U  )  z-™"1  +  B  (z)  U  z"m"2\ 

x  k  j  ,m-l-k  m-kv  '  jk 

k=o 


where  {V.}.  indicates  the  i'th  element  of  the  vector  V  .  But  on  again 

J  1  j 

using  (2.3.16)  we  have 


m- 1 

(3.3.3)  R  .  .  ( z )  =  m  U  .  .  z'm_1  -  -TV 

'  Jjmv  '  J  jm  j 

k=o 


(b  (z)  u  )  z  m  y 

m-k  jk' 


If  1  /  j  we  define 


q,.j(z)  =  q±(z)  -  q,(=0 


(3-3.6) 


,  t)  /I  101 
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and  again  make  use  of  (2 .3.16)  in  (3.3.3)  to  obtain 


(3-5<7)  Rijm(z)  =  qij(2)Uijm  ^  *  Z[qij(z)  '  (qio  •  qjo)zr]Uijm  Z 


-m-1 


m+r+l 


T. 


(i=m+l 


min([i,r+l) 

I  (qiu  ' 

k=|i+l-m 


q.JU.  .  .  +  (|i-r-l)U.  .  . 


p-r-2 

V" 

v 

k=o 


m-1 


zr'"-  y 

1“  k=o  1 


-m-2 


3.3.2  Two  Possible  Cases 


In  this  section  we  shall  illustrate  that  it  will  in  all  cases 

be  possible  to  express  the  error  vector  by  at  most  a  simultaneous  pair 

* 

of  Volterra  vector  integral  equations  .  The  error  bounds  we  shall  obtain 
are  considerably  sharper  in  the  case  when  we  can  express  the  error  by 
a  single  Volterra  vector  integral  equation  than  in  the  case  when  we 
require  a  simultaneous  pair  of  Volterra  vector  integral  equations  to 
express  the  error  vector.  Thus  it  is  desirable  whenever  possible  to 
use  a  single  Volterra  vector  integral  equation. 


We  have  already  noted  at  the  end  of  Section  2.3  that  by  our 
construction  of  formal  solutions  there  is  a  one-to-one  correspondence 
between  the  eigenvalues  of  Cq  and  the  vector  solutions.  Thus  whenever 
it  is  possible  to  express  the  error  vector  by  a  single  Volterra  vector 


* 

We  assume  throughout  this  section  that  end-points  of  integration  be 
taken  at  infinity. 


eb'i.'rpd  loyi'j  ©rfT  ,  enoJUx,up9  ltt*xg«*raJ  loJoov  mio JloV  3to 
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integral  equation  we  shall  refer  to  the  corresponding  eigenvalue  as  an 
extreme  eigenvalue;  otherwise  we  shall  call  it  an  interior  eigenvalue . 

Consider  the  eigenvalues  of  Cq  to  be  represented  as  points  in 
the  complex  plane,  and  to  be  enclosed  there  by  a  smallest  closed  convex 
polygon  j[.  If  an  eigenvalue  of  C  is  an  extreme  point  (i.e.  a  vertex) 
of  the  polygon  it  will  always  be  an  extreme  eigenvalue;  if  it  is  an 
interior  point  but  not  a  boundary  point  it  will  always  be  an  interior 
eigenvalue;  if  it  is  a  boundary  point  but  not  an  extreme  point,  then  to 
determine  whether  it  is  an  extreme  eigenvalue  or  an  interior  eigenvalue 
it  is  necessary  to  also  examine  the  eigenvalues  of  C^. 


FIGURE  1 


FIGURE  2 
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The  polygon  !]  and  the  differences  between  the  various  eigenvalues 

Azr+1 

of  A  are  left  invariant  under  the  transformations  W  =  X  e  , 

o  * 

W  =  T(z)X  (T(z)  as  in  Theorem  2.3.  l),  z  =  wl|,  u)  a  constant  and 
| co |  =  1.  Thus  we  may  assume  without  loss  of  generality  that  the  vector 
solution  of  interest  corresponds  to  the  zero  eigenvalue.  Moreover,  in 
the  case  of  an  extreme  eigenvalue,  e.g.  when  zero  is  an  extreme  point  of 
the  polygon  TJ”,  we  may  assume  that  all  other  eigenvalues  of  are 

within  the  closed  sector  |arg  A|  <  a  <  ^  jt  and  that  there  is  at  least 
one  eigenvalue  on  each  of  the  rays  arg  A  =  a,  and  arg  A  =  -a.  In  the 
case  of  an  interior  eigenvalue  we  assume  that  no  eigenvalue  (other  than 
the  zero  eigenvalue)  is  located  on  the  imaginary  axis;  that  there  are 
some  eigenvalues  of  Aq  in  each  of  the  sectors  |arg  A|  <  a  and 

|rt  -  arg  A|  <  a,  0  <  a  <  “  it  and  that  there  is  at  least  one  eigenvalue 

on  each  of  the  lines  arg  A  =  a  and  arg  A  =  -a. 

Suppose  for  example  that  we  have  the  extreme  eigenvalue  case 
indicated  in  Figure  1;  that  C(z)  is  regular  (holomorphic )  for  all 
|  z  |  >  p;  that  r  =  0  and  that  =0  (i  =  1,2, 3, 4).  It  is 

then  clear  that  |e  |  (i  =  1,2,3>^)  increases  monotonically  along  any 

path  for  which  |arg  ~|  <  ~  -  a.  Thus  assuming  the  points  A^ 

(i  «  1*2, 3»^)  to  lie  In  the  z-plane  as  indicated  in  Figure  1  we  can  for 

example  connect  any  point  in  the  interior  of  the  overlapping  sector  of 
Figure  1  with  oo  +  0  \T-1  by  a  path  (P  consisting  of  at  most  two  straight 
lines  along  each  of  which  |arg  <  ~  -  a  and  neither  of  which  passes 


through  the  origin. 


n  eaoi  2  mantel  eiitt.1  t  ii  b.  u  *  ulifivai  Jltl  '**& 
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In  the  general  rank  r  +  1  (distinct  eigenvalue)  extreme  eigen¬ 
value  case  we  assume  C(z)  is  holomorphic  in  a  domain  oS  extending  to 
infinity  in  the  direction  arg  z  =  2jtk/(r  +  l)  for  some  integer  k  in 
0  <  k  <  r„  Let  £  be  the  smallest  closed  disc  with  center  at  the  origin 


. . o  ,n;  i  ^  j) ,  and  let 


containing  the  zeros  of  q^'(t)  (i  =  1,2, 
oQ*=o&-£  -  (oo) .  We  next  define  a  region  0&(  z,  to  be  the  union  of 

all  points  £  )  such  that  there  is  a  path  P  connecting  z  €  d$(z»  Ck) 

and  =  00  exP  [2itk  \f- T/(r  +  l)]  satisfying  the  following  extreme  eigen¬ 
value  conditions? 


N* 


1)  Except  for  the  point  P  lies  entirely  in 

2)  P  consists  of  a  finite  number  of  Jordan  arcs  t  =  t(s)  such 
dt  (  s  ) 


that  — ^ is  continuous  and  non- vanishing; 

3)  For  any  two  points  t  and  t  in  the  order  £  ,  t  ,  t 

-L  td  K  JL  cL 

z  on  P  we  have  | exp  [qij(t2)  -  j ( x ^ !  <  1  U  =  1,2,. ..,n); 


^  y$>  (t  1)  =  /  1 1  2  dt  |  i 

y  '1)0 


s  bounded. 


We  have  already  seen  that  for  a  particular  rank  1  case  with 

eigenvalues  as  in  Figure  1  there  always  is  such  a  path  connecting  each 

point  z  in  jarg  zj  <  *7^  -  a  -  5  with  +co,  where  5  is  an  arbitrary 

3 

positive  number  less  than  ^  jt.  For  arbitrary  rank  r  +  1  it  follows 

that,  since  jarg  q,  ,(z)  -  arg  [(qc  -  q .  )zr+^/(r  +  l)]j  can  be  made 

ij  10  jo 

arbitrarily  small  by  taking  jzj  sufficiently  large,  there  is  a  path  P 
satisfying  the  above  conditions  which  connects  any  z  (Jzj  sufficiently 


large)  in 


arg  z  - 


2itk  .  1 

r  +  1  !  =  r  +  1 


'2* 

2 


-  a 


and  £  =  exp 


&tk  sT-l 

r  +  1 
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Let  us  set 


(3.3.8) 


Dj(z)  =  Q'(z)  -  qj'(z)i. 


If  a  vector  Wj(z)  satisfies  (2.3.3) »  then  by  (3.3.2)  the  error  vector 


(3.3.9) 


6jm(z>  "  (Wj(z)  -^jm(z)) 


-q  M) 


satisfies 


(3.3.IO)  e.  (z)  -  D'.(z)  e.  (z)  =  B(z)  €.  (z)  -  R.  (z). 

*  dz  jm'  '  j  jm'  2  '  '  jm  jm'  ' 

z 


Now  if  € .  .  (z)  is  a  solution  of 

jm;kv  ' 


(3.3.II)  € .  .  (z)  =  r 

jm;kv  '  j 


z  D  .(z)  -  D  .(t) 


r  '[t“2  B(t)  €.  ,  (t)  - 


L 


jm;k 


R.  (t) ]  dt 

jm 


where  the  path  of  integration  is  chosen  as  described  above,  then 
Gjnrk(Z)  simultaneously  also  satisfies  (3.3. 10). 

Let  us  now  turn  to  the  interior  eigenvalue  case.  We  recall  that 
we  assume  all  eigenvalues  of  to  lie  in  the  shaded  region  indicated 

in  Figure  2. 


Take  for  example  a  rank  one  case  for  which  the  coefficient 
matrix  C(z)  is  holomorphic  for  all  |z|  sufficiently  large.  If  the 
eigenvalues  of  Cq  (or  A^)  are  as  indicated  in  Figure  2  then  given  any 
number  S  in  0  <  5  <  it  -  a,  there  is  a  path  '(p  joining  -00,  any  point 


z  in  |arg  z-—  jtjfit-a-b  and  oo, 


and  consisting  of  at  most 


(8.?.?) 


»  t*' 


,9VOC  >  >dr.  -’.£5  :!  .00,(3  8i  3 » Jf  t  1-  •  J o-  ttrfj  - ^ :  rf 
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three  straight  lines,  none  of  which  contacts  the  origin  such  that  Re(A^t), 

Re(A  t)  increase  monotonically  while  Re(A  t)  decreases  mono tonic ally. 

2  3 

If  in  addition  a  domain  about  the  origin  is  to  be  avoided  by  ^  |z|  may 
need  to  be  taken  large  if  5  is  small.  Similarly  it  is  easy  to  see  that 


there  is  such  a  path  joining  -  oo,  any  point  z  in  |arg  z  +  —  jt  |  <  jt  -  a 
and  oo  .  Thus  these  two  families  of  paths,  one  of  which  passes  around 
the  origin  in  a  positive  sense  and  the  other  in  a  negative  as  we  traverse 
from  -  oo  to  oo,  cover  the  complete  vicinity  of  infinity. 


-  & 


In  general  assume  C(z)  is  holomorphic  in  a  domain  yQ-  which 
contains  (or  extends  to  infinity  in)  one  of  the  sectors 

|arg  z  -  (2k  t  — )  ^y|  <  (?t  -  a)/(r  +  l)  for  some  integer  k  in  0  <  k  <  r 
and  all  z  sufficiently  large.  Divide  the  integers  N  =  {l,2,***,n}  into 
two  disjoint  classes  and  such  that  if  A^  is  in  |arg  A|  <  a 

then  i  e  N^,  while  if  A^  is  in  |arg  A  -  jt  |  <  a  then  i  e  N^. 

Note  that  j  may  be  taken  either  in  or  in  N^. 


We  next  reorder  and  if  necessary  relabel  the  elements  of  the 
equation  (3-3*10)  30  that  we  may  write  them  in  the  form 


tz  €jml(z)-Dh(z)ejml(z)  =  z’2[a(z)ejml(z)  +  P(z)ejm2(z)]-Rjml(z) 
ejm2(z)'Dj2(z)ej»2(z)  =  Z’2[7(z>6jml(2)  +  6(z^  jm2(  z)  ]-Rjn»2(z) 


(3-3.12) 


where  the  top  line  contains  all  rows  of  (3. 3 *1.0)  such  that  i  e  N^  and 
the  bottom  line  contains  all  rows  of  (3.3. 10)  such  that  i  €  N^.  The 
diagonal  matrix  D'.(z)  in  (3.3. 10 )  is  accordingly  subdivided  into  two 
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diagonal  matrices  D^(z)  and  D'o(z);  the  matrix  B{z)  is  accordingly 

partitioned  into  four  blocks  a(z),,  3(z),  ?(z)  and  5(z),  and 

Ro  (z)  is  aplit  into  two  vectors  in  the  same  manner  as  €.  (z)  was  split 
jm  jm 


Let  C-  be  the  smallest  closed  disc  with  center  at  the  origin 
containing  the  zeros  of  q!^(t)  (*•  =  l>2»...»n;  i  j)  and  let 
C&  =  oQ  -  £  -  {°°}.  We  next  define  a  region  <?$(£,,  $  £01  )  to  be  the 
union  of  all  points  z  (z  ^  z  /  su°b  that  there  is  a  path 

connecting  z  and  ( in  that  order),  where  is  one  of  the 

points  oo  exp  [(2k  +  l)  jt  v- T/(r  +  l)],  and  ^  ^  =  oo  exp  [2kjt  \f-T/(r  +  l)  ] 
satisfying  the  following  interior  eigenvalue  conditions 


that 


1)  Except  for  £lk»  p  lies  entirely  in  qQ  ; 

2)  P  consists  of  a  finite  number  of  Jordan  arcs  t  =  t(s)  such 

Mil  , 


ds 


is  continuous  and  non-vanishing; 


3)  For  any  two  distinct  points  t^  and  t  in  the  order  £ 

m  IP  we  have  | exp  [q..(t0)  -  q..(t  )]|<  1  if  i  £  N 

■t  .J  ^  *•  —  -t 


lk: 


*2k 


|exp  [qij(t1)  -  qij(t2)l|  <  1  if  i  €  n2? 

U)  With  a,  3,  7  and  6  defined  by  equation  (3°3°3C>) 
|  [a  +  &  +  /(a  -  5)^  +  437)  Yp  (t’1)  <  1. 


It  follows  that  the  vector  of  analytic  functions  which  satisfies 


q;.  '  '  : 


l:\Oi rfainsv-con  l  -„a  yw«i Inoa  ai 


■ 
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(3-3.13) 


ejmljk(z) 


£jm2;k(z) 


f'Z  eDjl(z)'Djl(t) 


'lk 


r  eDji(z)-Dj2(t) 

JL 


(2k 


/ 


-2 


a(t)  p(t) 


7(t)  B(t) 


G jmljk^ t ^ 


ejm2;k(t^ 


V(t) 


R.  0(t) 
jm2 


dt 


where  the  integrals  are  taken  along  a  path  as  described  above,  is 
simultaneously  a  solution  of  (3.3. 12  ). 


3.3.3  Treatment  of  R..  ( z ) 
— - 1  jnr — u 


In  this  section  we  shall  accurately  estimate  the  integral 


(3.3.14) 


R*.  (z)  =  F  e 
J-J"1  J  5 


1J  J  R. .  (t)  dt 

ljm 


where  the  path  of  integration  satisfies  anyone  of  the  sets  of  conditions 
described  in  the  previous  section  and  £  is  a  point  appropriately  chosen 
at  infinity  such  that  as  t  traverses  from  £  to  z,  Re  ^ ( t ) 
decreases  monotonically .  R^^(t)  is  given  by  (3.3.7). 


The  chief  difficulty  surrounding  the  above  integral  occurs 
for  i  ji  j.  If  i  =  j  then  R^  (t)  =  0(t  m  *),  |t|  ->oo;  if  i  /  j 
then  R^jm(t)  =  0(tr  m  ^),  1 1  j  — ►  00 .  Thus  when  i  ^  j  we  integrate 


.(!  o  ■  ivlea  fi  ynis 


taion  a  3; 

♦  r  iK  Jlnoionoar  aseaai 

[\i  tl  ;c  <-  |J|  .(  ‘i)0  (•)  .  ,  OJrt5  t  .  t  11  t\i  70} 


70  - 


(3.3*14)  by  parts  to  obtain  a  good  bound  on  Rjjm(2)* 


For  i  ^  j,  the  equation  (3.3*7)  may  be  written 


(3.3.15) 


ljm 


q,*(2)  ^  -q..^2) 

(e 


/  \  'ij '  '  d  /  iiv  'TT  -mN  ,  /  /  \  /  \  r  mlT_  -m- 1 

R,.4„(z)  =  -e  dz  ^  J  )-z[q,  Jz)-(q,. „-q.,J2  +“JU,..,  z 


ljm 


ij 


10  jo 


z  ijm 


m+r+1  min(|i,r+l)  p-r-2 

F  F  ^qik'qjk^Uij,p-k  +  (M"r"1)Uij>p-r-l  +  l  ^BkUj,p.-r-2-k^i 

p=m+l  k=(i+l-m  k=o 


r-|i 


m- 1 

y  { fvk^^jk)}. 

k=o  1 


-m-2  . 


Thus  the  first  term  on  the  right  of  (3.3. 15)  bas  a  form  suitable  for 
bounding.  Consider  now  that  polynomial  part  P^m(z)  consisting  of  all 
of  the  right  hand  side  of  (3.3*15)  excluding  the  first  and  last  terms. 


It  is  clear  that  ^-y^2)  *-s  a  polynomial  with  dominating  term  of  order 
zr  m  1  as  Izl  — > 00 ,  On  substituting  P..  (t)  for  R..  (t)  in  (3.3.14) 

1  1  0  ijm  7  ijm 

and  denoting  the  resulting  integral  by  P*^(z)  we  have 


*  p2  q, ,(2)"q, ,(t) 

(3-3=16)  PMJz)=  /  e  J  J  q/.(t)  (  Pijm^t)  ^  dt 


/ 


e 


so  that 


(3.3.17) 


*  /  X 

P.  .  (z) 

ijm 


P.  .  (z)  pZ  q.  .  (z)-q.  ,(t) 

uf  7  +  r  e  iJ  lj 

qij(z) 


_d_ 

dt 


W,(0  ! 


dt 


x  1  \  fir-  jP"  v ) #  »P 

®uu  *  *  Vo/Hs)up  Vi°  ' 


' 


^isq  J  ilfrroi  ,-J  7  3*rf:J  w  >n  lab  t » .  roO  ,ani6ft:;od 
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Collecting  terms  we  obtain 


(3.3.18)  a*jm(*)=/? 


z  r  _  -m-l  \  d 

m  Uijm  +  («  -  1  !dl 


iim^K 

qS(t) 


-  e 


m-l 


k=o 


Vk^V^  ' 


-m-2 


dt  , 


q. .(z)-q. .(t)  j  ‘ \  m  o 

and  since  |e  1J  1J  |  <  1,  and  I77  (  ^TT'T'  )  |  <  0  (|t|  ) 


,  /  P .  .  ( t)  \ 

J 1 1  ->-00,  the  right  of  (3.3*18)  is  easily  bounded  and  can  be  utilized  for 


purposes  of  either  a  vector  or  a  norm  bound. 


On  the  other  hand,  for  i  =  j  we  have  by  (3.3.5)  that 


(3*3-19)  R*.  (2)  -  r 

J  jmv 


Z  r 


m- 1 


m  U  .  .  t 
Jjm 


-m- 1 


-  T  Ivk^^jJ. c 

j 


k=o 


-m-2 


dt 


and  obtaining  a  bound  for  this  equation  is  even  simpler  than  obtaining 
one  for  (3*3.18). 


3.3.4  Error  Bounds  for  the  Extr erne  Eigenvalue  Case 

We  shall  first  obtain  a  vector  bound.  By  Section  B.2  of 
Appendix  B  the  solution  of  the  equation  (3.3. 11)  exists  and  is  a  unique 
vector  of  functions  holomorphic  in  a  Riemann  surface  which  contains 
the  region  ylQ  (z,£^)  defined  by  the  extreme  eigenvalue  conditions  in 


Section  3*3*2. 


(I),fp-(*h.p 


90nlB  Ins 


o  s  bob  cod  <  '  i  >  <*i  (  [,£.£)  O  3»,gli  orfJ 


.  ru/od  mion  a  no  io?o*v  a  ledila  !to  ssioqiuq 


•*t: 


l-H*?  (w*w> 


.(8I.’:.£)  toI  o  .> 


^  ■■  i  ,  r  i  r  ■  ■  .  <•.  .«■  i  r  I  .-(  , 

Urxs  {XI,  aoJtiMupa  srfd  Id  i  o*.  uIor  sd}  a  xlbnaq jA 
"•  ^  ;  '  >•  '  ti  trf  noj-  j(Vj1  So  io3 


•*•«•«  no13392 


. 
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With  t^  and  t  in  the  order  tp>  z  on  we  define 


(m+1)  7.  .  =  sup 


i  jm 


l,t:2  ** 


(3.3.20) 


iqij(t2)"qij(tl)r  m+2  d_  / 

1  dt  v  ^(0 ;  t=ti  + 

m-l  /  V 

Z  {v*  WJ  •  £  (%trr) « t  c2 

k=l  J  1 


i  =  1,2, . . . ,n;  i  /  j) 


(m  +  1)  7 


jjm 


sup 
t  €  P 


m-l 

Xivu'^J 

k=o  J 


We  also  denote  7,  to  be  the  vector  with  i’th  element  7  :  i„e( 

jm  ijm 

(7.3.  =  7. .  ,  and 
jm  1  ijm 


(3.3.21)  B  = 


sup 

VS  €iP 


q41(0-q, ,(t.)  r  -j 

*  1J  2  1J  1  iB(t.)l 

^  -  is 


(i,s  —  l,2,,.,,n). 


On  combining  equations  (3.3.20)  with  (3.3.21)  and  substituting 
in  equation  (3.3- ll)  we  obtain 


(3.3.22) 


'Wz)|  <S  D1'2'  B  l£jm;k(t)l  +  "'V'”"1!  + 


Jm 


dt 


We  apply  Lemmas  C.l  and  C.2  of  Appendix  C  to  (3.3*22)  to  obtain 


Theorem  3.3.1  If  corresponding  to  an  extreme  eigenvalue  of  A 

(or  C  )  we  can  define  a  region  oQ(z>  £,  )  as  described  by  the  extreme 
O  11  R  ’  '  " - 

eigenvalue  conditions  in  Section  3*3*2,  then  the  equation  (2.3*3)  possesses 


^  £  Wtv> 


•  i,w.4j  i ,  - 


,ttT  (X  +  .) 


,ctt' "■ ■■-  ;  ;  0  mv  ;.i  «tf  03  T  ^ joi  <>b  o«I»  sW 


q  .tS 


).  'u»  bn*  .{.(,)  '3fv  (•■:  .«  1  !ool:;i.upn  f.  IrJtdno)  *0 


£2  St  Hi  t  \  -f,mg-'iiP3  il 
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an  actual  solution  vector  W ,  ^(z)  depending  on  L  and  on  arbitrary 

J  5  t\. 

positive  integer  m  such  that 

m- 1 


(3.3.23) 


where 


W  .  ,  (z)  = 

jmjtc 


U .  zS+€.  .  (  z  ) 

js  jmjk' 


1  <iAz) 

^  -J 


s=o 


(3.3.24)  Njm;k(z)|  <  exp  (B^t-1))  n^(U.m  O  +T^(7.m  f”’1)) 


each  vector  U  ^  ( k  =  0,1 


e  e  • 

$  J 


whenever  z  is_  in  ^(9  (z>^)«  Moreover  i£  —  vector  °f  funct 

holomorphic  on  a  Riemann  surface  which  contains  ^  (z,?k).  In  (3.3.23) 

)  is  the  j ' t h  column  vector  of  the  matrix 
U,  defined  by  Theorem  2,3.2;  q  „  ( z  )  is  the  j ' th  diagonal  element  of  the 
diagonal  matrix  Q(z)  ( equation  (2.3.14)).  In  (2.3.24)  B  is  an  nxn 
matrix  of  non-negative  elements  defined  through  equations  (2.3«4)  (with 
B  (z)  =  B(z))  and  (3.3*21),  while  the  i “ t h  element  (i  =  l,2,*,#,n)  of 
the  vector  7^  jLs  defined  by  equation  (3.3.20). 


10ns 


Let  us  now  obtain  a  norm  bound  for  e.  ,  (z)  in  equation 

jm;kv  ' 

(3.3.23).  For  this  purpose  we  define  a  vector  V^(z,t)  by 

q1j(z)-q1jit)  im-2  d_  /fijmii'S 

-  (e  1)c  at  ^nzr)+ 


(3.3.25)  (V  (z,t)) 

J  1 


+  e 


qij(z)-qij(t) 


m-1 


I  tVkt^V3! 


k=o 


m-1 

=  Y  ^m-k^^jk^j 

k=o 


(i  =  l,2,”*,n;  i  /  j) 


(i  -  J) 


. 


-  7^  - 


Using  any  one  of  the  definitions  of  norms  in  Chapter  I,  and  with  t^ 
and  t^  defined  as  for  equation  (3.3. 20)  we  define 


B  =  sup  || e  J  J  B ( t )  |! 

t j  s  t2  €  ft 


(3.3.26) 


7  . 
Jm 


sup  ||V  (tBit1)||  . 

tlft  2elP  J 


Theorem  3.3.2 
3.3.I  is  given  by 


A  norm  bound  for  the  vector  Gjm<>k(z)  AE  Theorem 


(3.3.27)  ||e,m  t(*)||  5  exp  {BT^Ct'1))  {||U, 


jm 


jm 


7^(t 


where  B  and  7 .  are  non-negative  numbers  defined  by  equations  (3.3*25) 

and  (3.3.26). 


3,3,5  Error  Bounds  for  the  Interior  Eigenvalue  Case 

By  the  interior  eigenvalue  conditions  of  Section  3*3.2  the 
solution  to  equation  (3*3.13)  exists  and  is  a  unique  vector  of  functions 
holomorphic  in  a  domain  (actually  a  Riemann  surface)  which  contains 
the  region  JQ  ( ^2k)  *  *^e  Pro°^  °f  this  last  statement  is  given  in 
Section  B.2  of  Appendix  B. 

Using  any  one  of  the  definitions  of  norm  given  in  Chapter  I,  we 
shall  start  by  obtaining  a  norm  bound  for  the  solution  of  equation  (3*3*13). 


■ 


?A  *o  .Q  u  130^8 
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Let  Vj(z,t)  be  defined  as  in  (3.3.25).  We  split  the  vector 

Vj(z,t)  into  two  vectors  V^(z,t)  and  V^(z>t)  in  a  manner  corresponding 

to  the  way  in  which  R  (t)  was  split  for  equation  (3* 5* 13)*  We  analogously 

split  the  vector  U.  into  U.  .  and  U.  such  that  equation  (3.3.13) 

r  jm  jml  jm2 

may  be  written 


*»  * 

<  — 1 

PZ 

ejml;k^z) 

0 

^lk 

(3.3.28) 

— 

P  z 

- 

£jm2;k(z) 

- 1 

CM 

■*_n 

O 

_ 1 

Dji(z)'Dji(t) 


0 


0  e 


B.2(z)-DJ2(t) 


-2 


<x(t)  p(t) 


/(O  8(t) 


‘j-S.kW 


TT  -m-1  „  /  . 

"m  Jjmlt  +  Vjl(z>t)t: 


TT  -m-1  „  /  \  -m-2 

-m  U.  t  +  V.0(z,t)t 
jm2  J2V 


>  dt 


J 


We  furthermore  define 


(O  +  sup  ||  ^  V  (t  ,t  )||  (t-"1'1) 

lk,  z  lk,  z 


(3.3.29 ) 


=  hv  v^gk(t"m)  + 


sup 

tl>t2 


shvti*t2)iiy'..52k(t‘B'1) 


where  t^  and  t^  are  any  two  points  on  l(p  in  the  order  tp>  £pk 

and  the  variations  are  taken  along  ^  .  Similarly  we  put 
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D„(t2)-D  (tl)  D  (t2)-D  (tx) 

a  =  sup  || e  J  J  a(t  )||,  (3=  sup  ||e  3  3  0(O|| 

>  t-2  sty  ,  t2  €  ^ 


(3.3.30) 

D^p(tl)“Dip(t2)  D^pCtl )~V±n(t2) 

7  =  sup  ||e  3  32  7  ( t  )  || ,  &  =  sup  || e  j2  3c~ 

ti,t 2  e $  ti.ta  e  ft* 


6(t2)|| 


Applying  the  analysis  of  Section  B.2  to  the  above  equations  we 
obtain  the  following  theorem. 


Theorem  5.3.3  If  corresponding  to  an  interior  eigenvalue  A  oj! 
Aq  (or  Cq)  the  interior  eigenvalue  conditions  of  Section  3.3.2  are 
satisfied,  then  the  equation  (2.3.3)  possesses  an  actual  solution  vector 


(3.5.31) 


w.  -  .(z) 
jml ;k  7 

m-1 

y 

U 

jsl 

/ 

c _ J 

s=o 

%2;k(z) 

k- 

U.  0 

-  Js2J 

-s 

Z  + 


€  .  -  .  ( 2) 

jml ;kv 


,ejm2;k(z' 


qj(z) 


at  each  point  a  of  O^)  ’  w~rier& 

(3.3.32) 


< 

^ml^lk’2^ 

4* 

vf  z(t_i)  0  ‘ 

nk,z 

E  f(V 

Y'l)] 

^ml^lk’^k^ 

Jlejm2;k(2)H_ 

^n,2(2’C2k) 

0  Yz  t  (t_1) 
Z,q2k  J 

^m2^1k^2k^ 

*  0 

and  where 


. 


'irtQ*rtJ  aoi^oric^  drfi  nJtajdo 
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(3.3.33)  E  f (Ex)  = 


1  -  (oh-5)x  +  Ax' 


a  -  AX 


7  5  -  Ax 


;  A  =  a5  -  |37  . 


The  bound  (3.3*32)  is.  valid  for  all  z  such  that  —  [a  +  S  +  \l (a-5)2+4(37 )°K( t  ^)<1 

2  tP 

The  elements  of  the  solution  vector  (3.3.3I)  depend  on  an(^  an 

arbitrary  positive  integer  m;  each  element  in  this  vector  is  holomorphic 

in  a  Riemann  surface  which  contains  (^k^  *  The  function  q^(z)  f s 

the  j'th  diagonal  element  of  the  diagonal  matrix  Q(z)  ( equation  (2.3*1^)). 

The  vector  U .  ,  the  j ' th  column  vector  of  U  (s=0,l,2,***)  (defined 
J  s  s 

in  the  proof  of  Theorem  2.3.2)  .is  partitioned  in  a  manner  described  in 

Section  3.3.2  into  two  lower  dimensional  vectors  U.  „  and  U. 

-  >  > - —  jsl  -  js2 

The  functions  if/  (,)  (s  =  1,2)  are  defined  by  equation  (3.3.29) 
and  the  numbers  a,  3,  7  and  5  are  defined  by  equation  (3.3.30)* 


Next,  we  shall  obtain  a  vector  bound  for  the  solution  of 
equation  (3* 3*13).  To  achieve  this  goal  we  proceed  directly  to  bounding 
the  elements  of  vectors  and  matrices  in  equation  (3.3.28).  Thus,  with 

a(t)  =  [(a(t)].s(i  e8j.se  Nj)],  p(t)  =  [{p(t)]is(i  £  8,,  s  E  Ng) ] 

(3-3-54) 

7(0  =  [{7(t))is(i  €  N2.  s  €  Nj ) ] .  S(t)  =  [{6(t)).s(i  €  N2>  s  €  Ng)] 


and  t^  and  t^  as  defined  for  equation  (3.3.29)  we  define 


* '  •  -  A  J,  ■  Si  £  £  ’jd-iSKl.S  iU'?!  »•’*  2J5LSL' £U.  $  ■'■- 3 :  og  XI 122  '  <  :;t  ' 
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D  (ta)-D  (tx)  D  (ta)-D  (ti) 

a  =  [  sup  |  {e  J  J  a(t  ))  |]  ,3  =  [  sup  |  [e  3  J 

ti,t 2  e  tx,t2  e  & 

(3.3.56) 


7  =  [  sup  |(e 

ti,t2  €^3 


Dj2(tl)-Djl(t2) 


7(t  )]  | ],  8  =  [  sup  | {e 

2  t1,t2e^ 


D.0(ti)-DJ0(t2) 


J2 


J  s(t2)].s|] 


and 


-m-1 


‘  m!Ujmll  I*  I  +  s“p  I  W^l  I1 

tnt2  €  (r 


-m-2 


(3.3.57) 


<Wt) 


=  m|U 


-m-1 


jm2iit  1+  7  jyvVii1 

cu  t2  €cr 


-m-2 


where  t^  and  t  are  points  on  IP  located  as  described  for  equation 
(3.3.29).  In  addition,  each  element  of  the  diagonal  matrices 


D1s(z)'Dis(t) 

eJ  J  (s  =  1,2)  has  an  upper  bound  of  1.  On  combining  the 


above  results  we  obtain  from  (3.3.28) 


(3.3.38) 


€jml,k^Z^ 


€jm2,k(Z) I 


< 


r 


j 


0 


’lk 


0  r 


2k 


f. 


{ 


a  3 


7  8 


€jml,k(t'l 


€jm2,k^t^ 


1 

lt~2l  + 

Vc> 

•  *** 

J 

I 


dt 


The  inequality  (3.3.38)  is  solved  in  Section  C.3  of  Appendix  C 
Collecting  the  results  pertaining  to  the  solution  of  this  inequality  we 


have 


io,  ,vp 9  toi  6.0  oe  b  hie  ^  ^  no  ainxoq  #ia  j  bis  3 
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Theorem  5.5.4  If  all  eigenvalues  of  the  matrix  ^ see 

equation  ( 3 . 3 • ^0 ) )  are  less  than  one  in  magnitude  then  the  vector  e  .  ,  ( z ) 

jnil  j  k 

( equation  (3*  3*  31 ) )  is  also  bounded  as 


(3.3.39) 


l£jml,k(z>l  5  exp  fa^lk>Z  ^ml(?lk’z) 


where 


+  Cl(Slk'*)[expl8T^«'1)  ^m2  (?lk’W  +  B]> 


f  (u,v) 
ms 


.  r 


u 


C  (t) 
ms 


dt 


(s  =  1,2) 


(3- 3. to) 


B  -  [I2  '  C2  Cl^lk>^2k)]  V 

•[expfa'/tr1))^^?^,^)  +  C1(eiU,£2k)eXp(s7^(?-1)3^m2(eilt,e2U)] 


cl(Clk>z)  =  tt’1[exp(aV^kjZ(?"1))-I1  ]p;  Cg  “  S‘1[exp(B'l^?  1 ) ) -Ig] 


and  where  1^  designates  a  k  x  k  unit  matrix,  and  I  designates  an 
(n-Hc)x(n-K)  unit  matrix;  k  being  the  number  of  elements  In 
(Section  3*3*2).  Moreover ,  an  exactly  similar  result  holds  for  €.  ,  (z) 

j  in^f  y  k 


A  norm  bound  which  is  sharper  in  the  neighborhood  of  infinity 
than  that  given  in  Theorem  3*^*3  can  be  obtained  by  a  procedure  similar 
to  the  one  used  to  obtain  the  vector  bound  in  Theorem  3.1+.4  (see  Appendix 

c). 


. 


•  !j  '  f  3  C  ix  i\  amvix  :•  fj  ls  ^3 
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Theorem  5.5.3  With  the  notation  of  Theorem  3.3.3  let  Q  =  max  (a, 5). 


(3*3*41)  c  =  exp  +  •^^)7^(^"1))sinh  <  1 


ip 


2  '  V 


then  the  norm  of  the  solution  e  ^  jc(z)  ( equation  (3.3.3I))  is.  bounded  as 

(3. 3. >*2)  Hejmi,k(z)H  5  exP  falflk,z  (?"1)}  ^mi(hk’z) 

+  Y(Z)t^hk>W+7(1-C)'1[D^„a(hk>W+E!('m2(hk>^c)1? 


where 


Y(z)  =  tf&  exP((^  +  ^^fl)}sinh  f2— Y{ T  7(z~1)} 


2  ^lk’Z 


(3.3.43)  D(x)  =  “  [exp  (xT^CiT1))  -  1] 

E  =  ~T  ^(aA-b)  ~  2  r)  +  D(a)  exp  (6  T^"1)}]* 


It  is  noteworthy  that  the  conditions  of  Theorem  3*3*4  can  be 
tested  using  the  Gersgoren  theorems  ([12]p.  227)  an(*  that  they  (as  well 
as  the  conditions  of  Theorems  3*3*3  anc*  3*3*5)  can  be  satisfied  by 
restricting  ourselves  to  contours  sufficiently  far  from  the  origin. 


If 


■ 
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3.4  The  General  Case 

5.4.1  The  Differential  Equation  for  an  Approximation 

We  shall  obtain  an  actual  solution  vector  of  the  system  (2.4.4)- 
(2.4.5)  corresponding  to  the  v'th  formal  solution  vector  (as  defined  in  equation 
(2.4.22))  in  the  j’th  column  of  blocks  of  formal  solutions  as  defined  in  the 
proof  of  Theorem  2.4.2. 


Starting  with  the  partial  sum 


(3.4.1) 


m-1 

V(z)  ■  (  I  V 

k=o 


e  J 


Q.(z) 


where  is  the  j '  th  column  of  blocks  of  the  matrix  (k  =  0,1,*  ••) 

as  defined  in  the  proof  of  Theorem  2.4.2,  and 


(3.4.2) 


Q.(z)  =  q.(z)  I  +  J  log  z, 

J  J  M-j 


q.(z),  I  and  J  also  being  defined  in  Theorem  2.4.2,  we  define  R.  (z) 
J  M-j  Jm 

by  the  differential  equation 

,  -2  Mz) 

(3.4.3)  ~  cpjm(z)  -  (Q '  ( z )  +  z  B( z)  ]  cpjm(z)  =  Rjm(z)  e  J  . 


On  expanding  (3*4.3)  and  using  (2.4.10)  we  obtain 


(3.4.4) 


R.  (z)  = 

jnr 


r+m+1  min(jj.,r+l) 

V  I  V  (Q  ut 

/j  1  Ls  s  j»h-s 

p=o  s=o 


u+  Q  ) 

j.H-s  ^js' 


p-r-2 


m-1 


B  U 

s  j,|i-r-2-s 


-  y  (b  (z)u  ^  z'm": 
/_j  \  m-s  jk  J 


s=o 


s=o 


4.  )  ."is  *rf:i  lo  .199  non.  Ig  Inwlo:,  ne  nli;3do  J  U  rfa  »W 

jU  '  fl  iir  *!>  V  *  *?  V  r:  i  uioa  a/b vo'i  ,!..}'^  srfj  o3  gftJh  ioqaaiioo  (?  4  9) 

’  ■  r  r  •  ^  -o  •  -v'-  >  ,  ; 1  {  •  11  ( {$<».  4.  ) 

.3.4.5.  nrtiosrf?  lo  'fcooiq 


j«  i  vt.iT  q  arfl  rfJiw  8fti  JT8l2  ^ 
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where  U*  =  [U.  if  0<s<m-l;  0  otherwise},  Q  (:>= 0,  l^'^r+l) 

j  s  j  s  —  —  o 

is  defined  in  (2.4.7),  Q.  is  the  j  ’  th  block  of  Q  ,  and  B  (z) 

J  s  s  s 

(s  =  1,2,***)  is  defined  in  (2.4.8). 

It  is  necessary  to  examine  the  (i,j)'th  block  R..  (z)  of 

J  ^ 

equation  (3*4.4).  On  making  use  of  (2.4.21)  in  the  (i,j)'th  block  of 
(3,4.4)  we  obtain 


(3.4.5)  R..  (z)  =  q'  (z)U..  z  m  -  [q/.(z)-(q.  -q.  )zr  P]U..  z  m 

ljm  ij  7  ljm  lj  ip  Jp  ljm 


m+r+1 

V" 

L 

(i=m+p+l 


min(|i,r+l) 

^qis  q  js  ^Ui  j  ,  (i-s  +  JpiUij  1  (i-r-1  Ui  j  ,|i-r-lJ|i^ 


s=|i+l-m 


(i-r-2 


+  (n-r-l)U.  4  ,  +  V  (b  U.  _  \ 

ij, H-r-1  [  s  j,|i-r-2-sfi 

s=o 


m-l 

-  y  -Tb  (z)  u.  \ 

L>  {  rn-sv  JS 


-m-2 


s=o 


if  i  ^  j  and  for  some  p  in  0  <  p  <  r  ,  q^-q  /  ot^er  ^an<^ 


if  q.p  =  qjp  (0  <  p  <  r),  then 


(3*4.6)  R..  (z)  =  [(q.  .  -  q.  .  +  m)U  .  -  U . ,  J  +  J  U 

7  ijmv  '  i,r+l  j  ,r+l  7  ijm  ijm 


ijm 


]  z 


-m-l 


m-l 


-  y  {B  ( z)U .  } 

/  ,  m-s  7  js  i 


-m-2 


s=o 


where  (3.4.6)  also  includes  the  case  i  =  j. 


' 


/  'V  cil  b*  **9f>  al 


lo  Mooli  ril’(t.i)  Jri.  nf  (iS  d.S)  lo  isu  jctte  aQ  .(4.4. j) 

oi*Jdo  (4.4. {) 


ffl-i  >Wl4®« 


i  :  Uq  orfi  riO 


9860  df(.7  29bula«i  0«X«L  /3,4,f)  dtSffw 
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Let  a  block  W^(z)  of  (a  vectors  of  holomorphic  functions  satisfy 
the  differential  equation  (2.4,4)-(2.4.5).  The  error  block 


(3.t.7) 


Tl.  (z)  =W.(z)  -cp.  (z) 
jm  J 


satisfies  the  equation 


(3.4.8)  ~  4jm(z)  -  CQ(Z)  +  z"P  B(z)]t^m(z)  =  -  Rjm(z)  e  J 


Q  .(z) 


where  R^(z)  is  given  by  (3.4.5)  or  (3.4,6)  above. 


For  later  convenience  we  partition  the  equation  (3.4.8)  into  two 
simultaneous  blocks 


(3-M) 


dz 


* 

V(z) 

V,(z )+w.y>  o' 

-2 

a(z) 

0(z) 

’hmd2) 

+  z 

.’Wzl 

0  D^zJ+I^z) 

7(z) 

5(z) 

11 jm2 ^ 

V(z) 


R .  (z) 

jm2v  ' 


Vz) 


where  n.  .(z'l  contains  all  those  i'th  blocks  for  which  i  e  N,  ;  q.  „(z) 
jml  '  1  jm2 

contains  all  those  i'th  blocks  for  which  i  €  N  where  N  and.  N  are 

C.  L  C, 

two  disjoint  sets  with  U  =  [l, 2, •••,&).  Precise  definitions  of 
and  Np  will  be  given  in  the  following  section.  The  matrix  B(z)  is 


doiitu.  to:  8?io  :  '  <  li*  at-  i  win  .  (*),  n  «,i»dw 
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accordingly  partitioned  into  blocks  <x(z),  3(z),  7(z)  and  &(z),  and 

R^m(z)  is  similarly  partitioned  into  R^^(z)  and  The  matrices 

D^(z)  and  D^(z)  are  partitions  of  the  matrix  D^(z)  =  Q(z)  ~  Iq^(z). 

5.4.2  The  Different  Possible  Cases 

Some  of  the  possible  cases  that  can  arise  in  the  most  general 
case  have  already  been  covered;  namely  the  case  of  a  regular  point,  the  case 
of  a  regular  singular  point  and  the  case  of  an  irregular  singular  point  when 
the  eigenvalues  of  the  lead  coefficient  matrix  are  distinct.  In  addition  we 
have  assumed  in  all  of  the  above  cases  that  the  end  points  of  integration  in 
the  integral  equations  defining  solutions  can  always  be  chosen  at  infinity. 

By  this  last  assumption  we  have  merely  assumed  that  the  coefficient  matrix  is 
holomorphic  in  certain  domains  extending  to  infinity.  However,  when  two  or 
more  eigenvalues  of  the  lead  coefficient  matrix  are  the  same,  or  more 
generally  for  C(z)  as  described  in  Theorem  2.4.1  it  will  not  always  be 
possible  to  achieve  end-points  of  paths  at  infinity,  even  when  C(z)  is 
holomorphic  for  all  |z|  sufficiently  large. 


In  our  treatment  of  the  general  case  we  shall  cover  all  the  possible 
remaining  cases,  including  those  already  covered  in  previous  sections. 


It  is  demonstrated  in  [8]  that  for  each  fixed  j  the  complete 


vicinity  of  infinity  can  be  divided  into  a  finite  number  of  sectors  S 


k,j 


(k  =  1,2,*,#,N)  such  that  for  all  i  (i  =  1,2,  •••,<£)  for  which  there 


is  an  integer  p  in  0  <  p  <  r  such  that  q  -  q  ^  0  we  have  either 
Re  q^(z)  bounded  above  or  else  Re  q^(z)  tending  to  +  oo  as  |z|  ->oo 


,JDt  »  i'>  9-t6  XjU  ,f!l  Jft  :,':9m»03  bB9l  SrfJ  Jo  a®l/rBVU3<jl9  >rfU 


■ 

ro  ov  I  I.  -  iw  ,i  >3*/oH  j  ni-'ni  oJ  rrtbf  .  -  arrlainol  nXu J  90  nJt  ,>M<pro~io  £orf 
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in  a  particular  sector  S,  . .  This  fact  enables  us  to  classify  all  the 

kj 

remaining  cases  under  two  different  headings:  (a)  those  cases  for  which  it 
is  possible  to  define  an  actual  solution  vector  by  a  single  Volterra 
integral  equation;  and  (b)  those  cases  for  which  we  require  a  simultaneous 
pair  of  Volterra  integral  equations  to  define  a  solution  vector. 


In  order  to  characterize  cases  (a)  and  (b)  above,  we  first  establish 

our  notation.  {£}  e  c&  (  cB  is  the  domain  in  which  c(z)  (equation  (2.b,h)) 

is  holomorphic)  will  denote  the  fixed  limits  of  integration:  in  case  (a), 

{£]  =  £;  in  case  (b),  {£}  =  £  ,  £  .  <5  will  denote  a  closed  disc  with 

center  at  the  origin  and  radius  a.  If  for  all  i  (i  =  1,2, q.  = 

■*“P 

q^p,  p  =  O,l,o..,r,  then  a  =  0  if  every  member  of  {£)  is  unbounded,  while 

if  min  tUI)  <  oo  we  set  a  =  min  {|£|}.  If  in  case  (b)  all  members  of  (£} 

are  bounded  then  we  assume  £  ^  £^  an(*  a  =  |£^|  =  j£,p|.  If  for  at  least 

one  i  e  (l,2,...,i)  we  have  q  ^  q  for  some  p  €  (0,l,...,r),  let 

a  denote  the  maximum  modulus  of  the  set  of  zeros  of  q'  (z)  for  all  such  i, 
o  ij 

set  a  =  a^  if  min  {|£|}  =  00  while  if  min  {|£J}  <oo  we  choose  {£} 

such  that  min  {|£|}  >  a^  and  no  member  of  {£}  coincides  with  a  zero 

of  q!  ,(z) .  Let  e$*  -o9-L  -  S  where  S  is  empty  or  [oo  e1^!  -  oo  <  Q  <  oo] 
ij  a 

according  to  whether  or  not  r  <  -2. 


Case  (a)  above  is  related  to  the  extreme  eigenvalue  case  of  the 
previous  section.  It  jls  characterized  as  follows .  We  define  a  region 
£)  to  be  the  union  of  all  points  z  £)  such  that  there  is  a 
path  connecting  z  e  cB{  z,  £)  and  £  which  satisfies  the  following 
conditions ; 

l)  Except  for  £,  if  £  €  <$*  -  0$  y  lies  entirely  in 


rfiirfv  30*  asa/io  aeorfi  (d)  bo*  ;no23aopa  19359301 


ic  .  v  noi  i'll  3  £>a.ti9.  3  i^au;  >  J>J3g93nl  ai-  >iXoV  w;n 


«nol3&3oa  loo 


rfJiw  ssxb  twRoIo  6  JOTSb  mw  i)  .,J.,3«(3)  ,(d)  nu  ol 
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that 


2)  F  consists  of  a  finite  number  of  Jordan  arcs  t  =  t(s)  such 

dt-lsl  , 


ds 


is  continuous  and  non-vanishing; 


3)  All  elements  of  the  matrices  exp  [Q  (t)  -  Q.(t)  -  I  [<1.(t)  - 

l  J 

qj(t)])  (i  =  1,2,. ...t;  i  /  j)  are  bounded  for  all  points  t,  r  in  the 
order  £,  ,  t,  t,  z  on  (P  ,  If  J  /  0,  or  if  for  some  i  (i  =  1,2,..., .0; 
i  jt  j)  we  have  q^  =  q  for  p  =  0,1,..., r  (or  both)  we  also  require 
that  T/t  be  bounded; 

(t^  ^)  is  bounded  for  each  fixed  q  where  q  may  be 
taken  zero  if  =  0  but  otherwise  q  is  an  arbitrary  positive  number 

less  than  1. 


It  is  shown  in  Appendix  B  that  if  0&(z,  £)  contains  a  point 
zq  (zq  bounded,  £  £)  then  z,  £)  is  in  fact  a  non-empty  open 

connected  domain. 

We  observe  here  that  the  region  ^$(z,  £)  defined  by  the  above 
conditions  includes  the  corresponding  region  defined  in  Section  3„1  and  3-2.2 
in  the  case  when  the  general  system  reduces  to  that  at  a  regular  point,  or 
a  regular  singular  point,  and  it  similarly  includes  the  region  defined  by 
the  extreme  eigenvalue  condition  in  Section  3.3.2  in  the  extreme  eigenvalue 
case  when  the  lead  coefficient  matrix  of  the  general  system  has  all  its 
eigenvalues  distinct. 


Case  (b)  above  jjs  related  to  the  interior  eigenvalue  case  of  the 
previous  section.  I_t  ij3  characterized  as  follows .  We  define  a  region 
Tl)  such  that  there  is  a  path  S'3  connecting  £^,  z  and 
(in  that  order)  satisfying  the  following  conditions: 

l)  Except  for  and  £  which  may  be  boundary  points  of  , 

lies  entirely  in  o0  ; 
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such  that 


2)  consists  of  a  finite  number  of  Jordan  arcs  t  =  t(s) 

dt(s) 


ds 


is  continuous  and  non-vanishingj 


3)  It  is  possible  to  subdivide  the  integers  i  =  1,2, ...,i  into 


two  non-empty  disjoint  classes  and  such  that  for  all  points  t 

and  r  in  the  order  t,  r,  on  >  all  elements  of  the  matrix  Q^(t, 

exp  [Q .  (t)  -  Q  (t)  -  I  [q . (t)  -  q.(t)]}  are  bounded  if  i  e  N  (i  /  j) 

1  1  J  J 

while  if  i  e  (i  ^  j)  then  all  elements  of  the  matrix  Q^(t,  T)  are 
bounded.  If  J  ^  0,  or  if  for  some  i  (i=l,2,...,i;i^j)  we  have 

"j 

q.  =  q.  for  p  =  0.1,.... r  or  both  then  we  assume  without  loss  of 
HiP  JP 

generality  that  both  i  and  j  are  in  N^,  that  is  bounded  and  that 

t/r  remains  bounded  for  all  points  t,  r  in  the  order  t,  t,  ^  on  ; 

4)  ~  [a  +  S  +  \T(a  -  &)2  +  407  }l^(t^  ■*■)  <  1,  where  q  may  be 
taken  to  be  zero  if  J  =0,  but  otherwise  q  is  an  arbitrary  positive 
number  less  than  1,  and  a,  3,  7  and  5  are  defined  by  (3.4.31). 


It  is  clear  that  in  the  case  when  the  eigenvalues  of  the  lead 

coefficient  matrix  of  the  general  system  are  distinct  then  the  domain 

,  0)  defined  by  the  above  conditions  includes  the  domain  (?$(£  ,  £  .  ) 

JL  c.  lK 

defined  by  the  interior  eigenvalue  conditions  in  Section  3- 3*2. 

We  have  shown  in  Appendix  B  that  if  q)  contains  a 

point  z  (|z  |  bounded)  distinct  from  £  and  t  then  £$(£  ,  £  ,  q) 

o  o  12  12’ 

is  in  fact  an  open  connected  domain. 

If  r  >  0  (equation  (2.4.4))  we  can  again  represent  the  eigen¬ 
values  of  Cq  as  points  in  the  complex  plantand  there  enclose  them  by  a 
smallest  strictly  convex  closed  polygon.  As  opposed  to  the  distinct  eigen- 
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value  case  it  is  now  possible  to  have,  several  formal  solution  vectors 

corresponding  to  a  particular  eigenvalue  q  ,  of  C  (see  e.g.  equation 

o  j  o 

(2.4.22)),  Nevertheless,  if  q  is  an  interior  point,  and  not  a  boundary 
point  of  the  polygon  it  will  be  necessary  (if  at  least  one  end-point  of 
integration  is  taken  at  infinity)  to  express  the  error  vector  by  a  simultaneous 
pair  of  Volterra  vector  integral  equations.  If  q  is  an  extreme  point 
of  the  polygon  and  only  one  formal  vector  solution  corresponds  to  q  it 
is  again  possible  as  in  the  previous  section,  to  express  the  error  vector 
by  a  single  Volterra  integral  equation.  If  q^  is  an  extreme  point  of  the 
polygon  and  more  than  one  formal  vector  solution  corresponds  to  q^,  or 
else  if  q^  is  on  an  edge  of  the  polygon  (and  not  necessarily  an  extreme 
point)  we  can  always  determine  whether  we  can  express  the  error  vector  by  a 
single  or  a  simultaneous  pair  of  Volterra  integral  equations  by  testing 
whether  or  not  the  above  generalized  extreme  eigenvalue  conditions  can  be 


satisfied, 


5.4.5  General  Treatment  of  (z) 


We  shall  accurately  estimate  the  integral 


(5-4.10) 


pz  Q  (z)-Q.(t)  Q  (t) 

R*.  (z)  =  -  /  e  1  1  R..  (t)  e  J  dt 

ljm  ,/ ^  ijml 


where  (t)  is  given  by  equation  (5.4.5)  or  (5.4.6),  and  the  contour 

integral  is  taken  along  any  one  of  the  paths  specified  in  Section  5,4.2. 
In  order  to  conveniently  evaluate  the  dominating  term  in  this  integral  we 
write  (5.4.5)  in  the  form 
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(3.4.11)  R. .  (t)  =  [Q'.(t)U  +;U..  -  U.,  Q.(t)J  t 

'  7  ijm  l  iim  t  nrn  ijm  i'  ' 


-m 


ijm  t  ijm  ijm  j 


m-1 


-m-2 


k=o 


where 


(3.4. !2)  PljB(t)  =  -  ([S,p+r‘iJ>p+1)tr'p'1+  •••  +(q.  +  “t'1 


i  jr+1  j  ,r+l 


+  J  t"1]  U.,_  -  U, J  t~1^)  t~m 


ijm  ijm  ji 


m+r+1  min(r+l ,u) 


X  A  ^jk^ij  ,n*k+J|i.Uij  ,|i-r-l  Ui  j  s pi-r-1  Jp. 


|a=m+p+l  k=p.+  l-m 


p-r-2 


+  (n-r-l)U..  .+  )  (B.D,  _  .}. 1 

'  ij  k  j  ,(i-r-2*k  i  ’ 


,r-n 


k=o 


The  first  term  in  (3.4.11)  is  already  in  a  form  such  that  it  can 

be  easily  bounded.  Assuming  that  for  some  integer  p  in  0  5  P  <  r  we 

have  q.  /  q.  ,  let 
IP  '  JP 


(3.4.13) 


».  r  2  Q.,(z)-Q.(t:)  Q  ■  (t) 


dt 


ijm 


where  P  ^(t)  given  in  (3.4.12).  It  is  convenient  to  define  S^(z,t)  by 


q^CzH^Ct)  Q1(z)"Qi(t)  Qj(t)~q,(z) 

(3.4.14)  S  (*, t)e  iJ  1J  =  e  1  1  P  (t)  e  J  J 
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Clearly,  knowing  the  form  of  Q^(t)  the  matrix  S„(z,t)  can  be  explicitly 
expressed.  In  any  case  S^(z,t)  is  nothing  more  than  P^m(t)  multiplied 


by  a  matrix  polynomial  in  log  (  — j  . 


Thus 


(3.1*.15)  PtJ  («)*  -/£  s^z.Oe  J  J  (-qij(t)  dt 


qij(z)-qij(t) 


"♦  r  j.  f-H 

■  6  Ji  dt  Wj 


_  ](z,tK  <^(0^(0 

Jr  dt  \  q'  (t)  )e 


so  that,  if  for  some  integer  p  in  0  <  p  <  r  we  have  -  q  0,  then 

t=z 

t-£ 


* 


(3.4.16)  Rljm(z)  =  e 


Qi(2)-Qi(t)„  .m  Qj(t) 


0,,  t  e 
ijm 


q1(z)-qij(t)  -,t=z 

-  S.,(z,t) 


e 

+  - T7 


(t)  -ij 


t=5 


r 


2  r 

5T 


’ll 


(z.t) 


.m-1 


IP 


k=o 


-  e 


,Q1(t)] 


dt 


On  the  other  hand  it  follows  readily  from  (3.^.6),  that  if  i 
or  if  i  /  j  but  qip  =  q^  (0  <  p  <  r)  then 

*  Qt(z)-Q,(t)  Q  .(t)-itssz 

(3.4.17)  R  (z)  =  e  1  U  t  eJ 

ijm  ijm  t_y 


=  j, 


pz  r-*1 


-m-2  Qj^  . 
e  J  dt 


k=o 


> 

I  xjfc  ..law  sxij  (^).p  I©  ir*io3  mil  iwonrf  t\l  ;a»IO 

ttl  lai  >ny [ oq  xl«»m  a  * 


M.p- 


•-  ‘ l"-  ( m  f-'  5 


■fir*  p 


f  a  5  -  ■)  <  i  ^in  »  »  :o'l  i  ,  •  »:.<i J  or 

H  .-Jr...  - 


<*),'  (J)  9-1*)  P 

‘  *  V.'  •  - <■>. 


«  (dl.d.vM 


'*  9  T3  o  .55 


*  9  *  +  mmr 
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The  right  hand  sides  of  both  equations  (3.4.16)  and  (3.4.17)  are 
now  in  a  form  easy  to  bound,  under  the  conditions  of  the  paths  specified 
in  the  previous  section. 

In  the  next  two  sections  we  shall  obtain  the  error  bounds.  The 
analysis  is  carried  out  under  the  assumption  that  end  points  are  fixed 
bounded  or  unbounded  points. 


Considerations  which  motivate  our  procedure  of  bounding  the 

*  -q.(z) 

elements  R. .  (z)  e  J  (equations  (3.4.16)  and  (3.4.17))  are  the  following. 

ljm  *  -qi(z) 

We  observe  that  each  element  of  R,  .  (z)  e  J  is  of  the  form 

ijm 

(3.4.18)  P(z,£)  =  f(z)  -  f(£)  +  g(z)  -  g(£)  +  f  h(z,t)  t  m  2  dt 

where  f(z)  <  0(|z|  m| log  z | ^  V)  (see  equation  (2.4.22)),  g(z)  <  0(|z|  m  1 
|  log  z|^  v)  and  |h(z,t)|  <  O(|log  z|^  |  log  t|^)  |z|,  |t|  -»-oo,  where  k 
and  2  are  non-negative  integers.  The  reason  for  introducing  terms  like 
t  \  t  T^1  where  rj  and  are  non-negative  numbers  less  than  one  is  to 

conveniently  eliminate  logarithmic  terms  while  at  the  same  time  obtain 
good  bounds.  Hence 


(5.4.19)  |p(..oi  <  -  ?ni’ml 


,s(z)..-  s(S) 


f  h(z,t)t  T^1t^1  m  2dt 


+ 


abrr/jod  boos 
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|p(z,£)|  <  |— ~ — 1  |(ti  -  m)  r  t^1  m  *  dt  | 
1  1  =  '  1  '  J  ^  1 


I  r  feLCS-.^gt^-stO)  +  h(z,t)t^4  tTll"m'2  dtl 

■Jj  J 


sup 


f(T)  -  f(0 


i  |T  «>  V’li-.jm-'J  'y 


-m, 


f  Sup  Igfr)  -  g(C)  h(r,t)t  T'1.  / 

It.T  e  y  'Tni->n-X  _  ni  .  m  .  i1  I  V 


lli-m-1 


) 


where  in  the  last  limit  superior  t  and  r  are  points  on  in  the 
order  £,  t,  t,  z. 


Thus  let  be  a  fixed  positive  number  which  we  take  to  be 

zero  if  J  =  0  (i  =  1,2 ,***,i)  and  in  0  <  t)i  <  1  if  at  least  one 
^i 

J  4  0.  It  is  convenient  to  define  vectors  U^*(t),  vY^(t)  and  V^(r,t) 
"  '  jm  7  jm  jm  7 


^i 


by 


UCJV*(T)  _  I  i-iJ|T  6 _  ~6 _ 

n.  -m  n.-m 


_m  Qj(0-qj(T>ia»V 


C(t)  ■ 


Jll-m  _  ,^1-™ 

*  (i  =  1>2, • • • ,i) 

r  q.,(T)-q..(0  'iw 

jSij(T>T)/q£j(T)"e  J  J  Si^(r,0/qij(0  p 

(i/j)* 


ri1  -m-1  /,rj1-m-l 

t  -  b 


(3.4.20) 


=  0  otherwise; 


ij  '  4ij  '  hh-2 

J  J  t  + 


V  Q,(t)-q,(T)  -|0)v 

j  h  C*/J) 


C(T-C)  - 


k=l 


km-l 


-qi  r  Q.(t)-q.(T)  ->,a)v 


m+l-q-L 


L  lB-k(t)Vie  J 

k=o 


y 


otherwise. 


s-.-xr,  +j,  . 


. 


■ 
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In  (3.4.20)  (A.}WV  =  indicates  the  (u),v)'th  element  in  the 

J  ^  J 

(i,j)'th  block  of  the  matrix  A,  and  the  notation  (i  /  j)  to  the 

right  of  a  particular  equation  indicates  that  the  particular  equation 

contains  all  those  i'th  elements  for  which  there  is  an  integer  p  in 

0  <  p  <  r  such  that  q.  -  q .  4  0. 

=  -  Hip  HJP  r 

5,4.4  Error  Bounds  for  Extreme  Eigenvalue  and  Related  Cases 


If  a  block  T|^(z)  of  vectors  of  holomorphic  functions 


satisfies 


(5.4.21)  r]  (z)  =  f  exp[Q(z)-Q(t)]{t“2B(t)  r\  (t) 

J  £  Jm 


Q,(t) 

R.  (t)e  J  }dt 
jm  J 


where  the  path  Ip  of  integration  satisfies  the  extreme  eigenvalue  conditions 
in  Section  3»^-»2>  then  it  follows  that  ri jm( z )  also  satisfies  (5.4.8). 

q.(z) 

Let  us  put  n .  (z)=€.  (z)eJ  ,  and  let  us  consider  the  v'th 
'jm'  '  jm  ' 

column  vector  (z)  of  e.  (z).  With  reference  to  the  discussion 

jm  '  jm  ' 

surrounding  the  equation  (2.4.22),  the  integral  equation  for  €^m(z)  takes 


the  form 


(3.4.22)  e\  (z)  =  r  e 


z  D4(z)-D.(t) 


jm 


J  5 


S=V 


s 

where  R.  (t)  is  the  s'th  column  vector  (s  =  l,2,***,p.)  in  R.  (t), 
jmv  '  x  *  *  'j'  jmv  ' 

and  D  (z)  =  Q(z)  -  I  q  (z). 

J  J 


Thus,  using  (3.4.20)  and  taking  t  and  r  in  the  order 
£ ,  t ,  r ,  z  on  $>  we  define 


’  1  .  A •*.  ’  )  n  I 


ii  .iiionoa  «ttiewo»s*3  «c3  M  [3  e  .1*.  1,  „o,  T;  *  o  riJ^q  <d3 
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B  = 


M 


sup 

t,T  e  P 


D.(T)-D(t) 


B( 


t)t-^ 


ks 


(3.4.23) 


u 


U)V* 

ijm 


sup 

r  €  p 


U 


GOV* 

ijm 


ij 


sup 

t,T  €  P 


where  the  notation  {A}^  indicates  the  (k,s)'th  element  of  the  matrix  A, 
and  the  fixed  number  r>  is  zero  if  J  is  zero  but  is  otherwise  an 
arbitrary  positive  number  less  than  1. 


We  are  now  in  position  to  envoke  the  analysis  of  Section  B.l 
of  Appendix  B  to  obtain 


Theorem  5.4.1  If  corresponding  to  a  formal  vector  solution  W^(z)  °f 
equation  (2.4.4)  we  can  define  a  region  £$ ( z ,  £  )  ajs  described  by  the 
extreme  eigenvalue  conditions  in  Section  3.4.2,  then  the  equation  (2.4.4) 
possesses  an  actual  solution  vector 


(3.4.24) 


m-1 


s=v 


us.  (lo-g  z)S'V ) 
Jk  (s  -  v):  J 


z) 


of 


functions  holomorphic 


in 


such  that 


(3.4.25)  |«J.(.)|  <  exp  {»7^(t^,}{«;^-)  +  7]m%  (t^-""1)} 


for  all 


vectors 


z  in  0$  (z,£)>  where  the  matrix  B  and  the  elements  of  the 

v*  v  /  1  \ 

and  7  are  defined  by  equation  (3.4.23).  The  solution 
depends  on  £  and  an  arbitrary  positive  integer  m.  The  function 


(*)j»  (  )4a 


l  j  {  a;  p  l»<  u/i  on.  b  e 

4 


M 1B1  Y.  1  0.-  .  ■;! Ck.iv.  ■  ;  : .■ 


. 


< 


-  95  - 


q^(z)  i§.  defined  in  Theorem  2.4.1.  The  integers  |i  and  v  and  the 

g 

vectors  U.,  (k=  0,1,2,  ••• ;  s  =  1,2,  •  •  •  ,|i.)  are  defined  as  for  equation 


jk 

(2.4.22).  If  each  J  =  0  (i  =  1,2, 


j 


then  the  numbers  r)  and 


above  can  be  taken  to  be  zero:  if  J  =0  then  q  can  be  taken  to  be 


j 

zero:  otherwise  r\  and  are  arbitrary  positive  numb er s  less  than  1. 


V  /  \ 

Let  us  now  obtain  a  norm  bound  for  the  vector  €.  (z)  in 

-  / 

Theorem  3*4.1.  To  achieve  this  goal  we  assume  the  definitions  (3.4.20) 

for  the  vectors  (z),  V^(z),  V^(z,t),  and  define 

jm'  '  *  jm  jm' 


»AT)-DAt) 


B  = 


i  -  n 


sup  || e 
t,T6(P 


J  B(t)t_T1|| 


(3. it. 26)  «£  =  T«»P  B^(t)|| 


7V.  =  sup  ||vY1(t)  +  v!2(r,t)|| 
jm  t.reiP  Jm  Jm 


where  t  and  t  are  on 


IP 


in  the  order  t,  t,  z. 


Theorem  3.4.2  A  norm  bound  for  the  vector  €.  (z)  in  Theorem  3.4.1 

=  — - jm - 


is  given  by 


( 5.4.27)  ||^m(z)t|  <  exp  {Br^(tT'-1)l{u;:^(t^-m)  + 

where  the  numbers  B,  U^*  and  are  defined  by  equation  (3.4.26). 

J  ^  J 

The  numbers  T|  and  are  defined  as  in  Theorem  3*4.1. 


3.4.5  Error  Bounds  for  Interior  Eigenvalue  and  Related  Cases 

q,(z) 

Substituting  ^mk(z)  =  ^mk(z)  e  J  (k  =  1>2)  in  (3*^*9)  and 
integrating  as  in  (3.3. 13)  and  (3.3.28)  we  obtain 


ni  baalUb  «1  («)  p 


' 


70:  s>v  ‘>rfi  7oi  jbric;  ^  jr.j  *  £  UJdo  won  %u 
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i-r  oi 

,£jm2(z) 

i 

o 

l/T  ~  ^ 

D  N 

I _ 

„  «v>v> 

(3.4.28) 


’a(t) 

p(0 

dt  + 

r_  *  .  o 

v(z) 

_/(t) 

8(t) 

where  in  the  notation  of  (3.4.9) 


(3.4.29) 


*7v  m  % 

v(z)  = 


-  /  e  jk  Jk  R.  .  (t)  e  J 

J  r  Jmk 


dt  (k  =  1,2). 


For  purposes  of  bounding  the  solution  of  (3.4.28)  it  is  again 

convenient  to  use  the  definitions  (3.4.20).  We  thus  partition  each  of 

v*,  \  vl /  \  v2/  \ 

the  vectors  U.  (r),  V.  (r),  V.  (T,t)  into  two  lower  dimensional 

jm  jm  '*  jm'  ' 

vectors  U  Vjmk<T)’  Vji<T’t)  <k  -  X*2>  in  a  manner  analogous 

*  /  \ 

to  the  way  in  which  Rjm(z)  was  split  into  lower  dimensional  vectors 
Rjmk(z)  (k  =  1,2).  It  is  of  course  presumed  that  the  point  £  in  (3.4.20) 
is  appropriatly  changed  to  if  i  e  (k  =  1,2;  see  the  interior 

eigenvalue  conditions  of  Section  3.4.2). 


We  shall  first  obtain  a  norm  bound  for  the  solution  of  (3.4.28). 
It  is  thus  convenient  to  define  scalar  functions  ^/mjc(  »  )  by 


9 


(i)fl  (J)* 

.(s.x  -  i) 


( ^ 3  ►  :jo!3£jO r  9;  3  nl  979riw 


9^'  «i  (i  ^3,  « ?  J  3  noiTiiiot  j  jrri  i/._  I  »  p*>eoq7iKj  to *2 


»»j  w)  (j.  ) ;  ;  ,(  )Jv  .(,)>  ,<5 


.uogoi  „i  (S.X  .  *)  ,(.*),  ,(r)  ,% 


»io3o»y  Xsitolsnsmlh  1  .woX  oJrcl  JXXqe  *+.*  (*)  M  diiriv  nt  T«w  arfj  oi 

.(S.X  »  it)  (*)  *1' 


roXMJW  »rfi  »5  i2.X  .  *)  j*  »  J  »  ,2  oj  isgn.rfo  tX3*Xiqoiqqf,  „i  • 

*  •  *  0  nol3w.  D8  9ff  Toi  bnjicd  •r*c.  ■•-,  -.  sjro  JeiiS  Ila/fa  »W 
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jJC^rvc + 

c2  e  r  1 


(3.4.30) 


*  tl,;r.  pi 1  "S'1**’  * 


v-y  •  t  («’*’)  • 

ci  e  (r  <= 


+  tir.pf|Vi(tl)  +  V^(tl’t2)ll^^2  (t’1l’m’1) 


where  t  and  t  are  located  on  in  the  order  £  ,  t  ,  t  , 

L  d  JL  1  d  t 

the  fixed  number  is  chosen^  as  for  ( 3.4.20). 


and 


Similarly,  with  t^  and  t^  chosen  as  for  (3.4. 30)  and  r\ 
zero  if  J  is  zero  but  otherwise  t\  an  arbitrary  positive  number  less 


j 

than  1,  we  define  non-negative  numbers  a,  3,  7  and  &  by 


-7 


a  =  sud 


Du(t2)’Du(ti) 

»  J  J  o(t  )||, 


fcl»  t2  e  P  1  " 

..Q  D.^t^-D.^tJ 


3  =  sup 


(3.U.31) 


ti»t2  elP 


7  =  sup 

ti,t2  eP 


6  =  sup 


1*7 


t  -11  »,p(ti)-D1p(ta) 


PO^II. 


1  -  T, 


e  J27 


j2  7(t0)||. 


t£l_  D.2(tl)-D.2(t2) 


tl  » t2  €  IP  ^ 


B(t2)||. 


A  variation  is  of  course  possible,  i.e.  we  could  choose  one  particular  rj 
for  and  a  different  r\  for  For  simplicity  of  expression  we 

have  used  one  rj  in  (3.4.30)  as  well  as  in  (3.4.31). 


. 


•$  3  S3t 


X  «S)  «*>  8 19 dr.  i  ml*sb  av  tI  and* 


(i3)c>a-(si),ta 
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On  combining  the  above  definitions  with  Theorem  B.2  of 
Appendix  B  we  obtain 


Theorem  3.4.3  If  corresponding  to  a  formal  solution  vector 

W.  (z)  (c.f .  equation  2.4.22))  we  can  define  a  region  ,  £  ,  r| ) 

1  2 

as  described  by  the  interior  eigenvalue  conditions  in  Section  3.4.2, 
then  the  equation  (2.4.4)  possesses  an  actual  solution  vector 


'V 


(3.4.32) 


jml 


(z)\ 


_( z)y 

v  jm2v 


(f°g  z) 

(s-v): 


s-v 


+ 


qj(z) 


of  functions  holomorphic  in  (JD{ £  ,  £  ,  r) )  such  that 


^ml(h’z) 

Yi  0] 

3^ 

I—* 

tor 

ro _ 

(3.4.33)  < 

►  << 

>  +< 

JL 

,  E  f[E7^(t ,)'1| 

t  ► 

Met  „(z)|| 
U1  jm2v 

K2(z’r>s\ 

0  (t^'h 

^  2  J 

L  1  2 

v  j 

where  E  f(E^)  ijj  defined  by  equation  (3*3*33)«  The  elements  in  the 
solution  vector  ^jm^(z)  depend  on  (k  =  1,  2)  and  an  arbitrary 

positive  integer  m.  The  function  q^(z)  liL  defined  as  in  Theorem  2.4.1. 
The  vector  U^k  is.  the  s^  column  vector  in  the  jt^  column  of  b locks 
of  Uk  (defined  in  the  proof  of  Theorem  2.4.2);  Tt  jLs  partitioned  into 
lower  dimensional  vectors  as  described  in  Section  3 • 4 . 2 .  The  functions 
^jmk(  y  )  (k  =  1,  2)  are  defined  by  (3,4.30)$  q  is.  an  arbitrary 
positive  number  less  than  one,  or  zero ,  depending  upon  whether  or  not 


0 


(«<* .?) 


2  rf  t  Jil  J5i  ?jt;  -l  .  •  r.  t  ,  U  1030$-'  arfl 
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there  is  at  least  one 
numbers  a,  0,  7  and 


^  0  (i  =  1,  2,  ...  ,  i).  The  non-negative 

are  defined  by  equation  (3.4.31). 


In  order  to  obtain  a  vector  bound  for  the  vectors  eV.  „(z)  and 

-  jml  ' 

€jm2(Z)  in  e<*uation  (3.4.28)  we  put  a(t)  =  [{a(t)]gh  (s,h  =  1,2, ...,k)]; 
3(0  =  [{3(0}  (s  =  1>2,...,k;  h  =  k  +  1,  k  +  2,  ...  ,  n)]; 

7( t)  =  [[7(t))sh  (s  =  k+1,  k+2,  ...  ,  n;  h  =  1,2, ...,k)]; 

6(t)  =  [{8(t))  ^  (s,  h  =  k+1,  k+2,  ...  ,  n)]  where  by  our  partitioning 
of  the  matrix  B(t)  into  blocks  3(t),  |3(t),  7(t)  and  5(t),  the  matrix 
a(t)  is  k  x  k  ,  and  similarly  for  3(t),  7(t)  and  6(t)  as  indicated 
here.  With  r)  defined  as  in  the  above  theorem  we  define  matrices 
a,  3,  7  and  5  of  non-negative  elements  by 


a  = 


1  -  n 


7  = 


1  “  I 


sup  .f  -n  .  .. 

e  a(Cl)] 

sud  -n  ^  \  1  ( ^2)  j  -i  ( ) 


sh 


sh 


(3. i+. 31*) 


p  = 


1  -  n 


'  sud  -r,  D,p(tl)-Di?(t2) 

*  /(t2)4h 


5  = 


1  -  Tj 


SUD  -Tl  D4p(tl)"Dip(t2) 
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where  ti  and  t2  are  defined  as  for  equation  (3.4.30*  [(tt^i))^!  = 

a^i)  =  [{a(t1)gh  (s,h  =  1,2,...  ,  k)]  and  similarly  for  3 ( t x ) , 

7( t2)  and  &(t2)  as  above. 

Next,  with  the  vectors  ^*^(0*  Vjmk^  and  Vjrak^T,t^  (k=1»2) 
defined  as  for  equation  (3.4. 30)  we  define  vectors  ^m^(  >  )  (k  =  1,2) 
of  non-negative  elements  by 
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(3.4.35) 


:1>^PC!P{|V>l(t2)  +  VM(Vtl)lK1.«(t1r,'1) 


Vz>y  = 


sup 

?  €  P 


u*v0(t0)|A^  ,  (t’1l'm) 

jm2v  2n  '  z,£0v  7 


t  7  .>£<■.>  * 

ti,c2  ey  d 


where  the  limit  superiors  are  taken  over  each  individual  element  of  the 
vectors.  The  number  T):.  and  the  points  t^  and  t  are  as  described 
for  equation  (3.4.30). 


The  remaining  details  concerning  the  proof  of  the  following 
theorem  are  carried  out  in  Appendix  C. 


Theorem  5 .4 .4  If  all  eigenvalues  of  the  matrix  C 

equation  (3.4.37))  are  less  than  one  in  magnitude  then  the  vector 
(equation  (3.4.32)  is,  bounded  as 


(see 


(3.4. 36)  l£Jmi(z)l  5  exP  (t n*1)}^ml  (£rz> 

+  C1(^1,z)  exp  (t71'1)!  +  8 


where 
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B  =  [I2  -  C2Cl(?1,?2)]  C2  • 

(3.U.37)  [exp  +  Cl«l*yexp  (^'1)^a2(^)J 

Cl(^l»z)  =  a  1[exp{aT^  jZ(tn  1)}  -  1^3;  =  5  1  [exp(5j^(tr]  1)}-I2] 


and  where  1^  designates  a  k  x  k  unit  matrix,  and  designates  an 

(n  -  k)  x  (n  -  k)  unit  matrix,  k2  being  the  number  of  elements  in 
<x(t)  (equation  (3.4.9)).  The  vectors  $  ,  (  ,  )  (k  =  1,2)  are  defined 
by  equation  (3*4.35)*  while  the  matrices  a,  3,  7  and  5  are  defined 
by  (3.4.34).  The  numbers  T)  and  are  defined  as  in  Theorem  3.4.3* 

A  norm  bound  which  is  sharper  in  the  neighborhood  of  infinity 
than  that  given  in  Theorem  3.4.3  can  be  obtained  by  a  procedure  exactly 
similar  to  that  for  the  above  theorem.  The  proof  is  given  in  Section  C.3 
of  Appendix  C. 


Theorem  3*4.3  With  the  notation  of  Theorem  3*4.3  let  0  =  max  (a, &) 


If 


{(»  *  “"h 


<  1 


then  the  norm  of  the  vector  e^^(z)  (equation  (3*4*32))  is.  bounded  as 
(3.it.39)  IWjml(z)ll  <  exp  (^-z) 
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D(a  +  6)  -  2D(o)  +  D(a)  exp 
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CHAPTER  IV 


CONCLUSION 


In  this  chapter  we  summarize  the  contents  of  the  thesis,  and  state 
some  unsolved  problems. 

In  this  thesis  we  have  found  error  bounds  for  partial  sum  approxi¬ 
mations  of  formal  solutions  to  actual  solutions  of  the  system  of  differential 
equations 


(U.l) 


dX  h  .  ,  %  „  h 
—  =  t  A(t)  X  =  t 


s-1 

T 

k=o 


-k 


Ak  t'“  +  Ag(t)  t  S  X  (s  =  1,2,  *  •  •  ) 


2 

where  the  n  elements  of  the  n  x  n  matrix  A  (t)  (s  =  0,1,2,***)  are 

s 

holomorphic  and  bounded  in  a  domain  (which  may  be  a  Riemann  surface) 
extending  to  infinity  and  where  h  is  an  integer. 


The  first  step  in  solving  (4.1)  is  to  transform  this  equation  to 
canonical  form.  In  this  connection,  the  existence  of  an  n  x  n  matrix 
polynomial 

N 

(4.2)  T(t)  =  £  Tk  t'k/p 

k=o 

where  p  and  N  are  positive  integers  and  the  T^'s  (k  =  0,1,#**,N)  are 
n  x  n  matrices  such  that  when  the  transformation  X  =  T(t)W  is  made  in 
(4.1)  the  resulting  system  has  the  desired  canonical  form  (see  Section  2.4) 
is  established  in  [32], 


In  the  majority  of  the  cases  when  h  >  0  fractional  powers  will 
not  be  required  in  the  transformation  to  canonical  form;  instead  T(t) 


.  aural  dcroq  bsvloanu  amou 


i  it  lo  o  :  Toe  t/w< j  ,  oj  anoltuloi  Ifu  :o\  lo 


arxut  on*.;!  tS  a  x>  oi.rfaiiJw)  >  .  alar  ; .  b*;*  ro<  bos  slrfqi.  -cojo.l 


c  1.  .-j  ,  i  crao-  jJ  [.«!)  tr  tw:.s  ni  ,  *.*«  j.-.  ariT 


I  til  "in-' »  svDKoq  9i%  Vf  btr  q  b  fy 


-  104  - 


will  be  of  the  form 


r+1 


(U.3) 


k=o 


with  (k  =  0,l,*#*,r+l)  an  n  x  n  constant  matrix.  The  resulting 

transformed  system  has  distinct  eigenvalues  in  the  lead  coefficient  matrix, 
and  has  an  irregular  singularity  of  rank  r  +  1  =  h  +  1  at  infinity.  Other 
possibilities  can  however  also  arise;  well-known  ones  among  these  are  the 
case  when  the  coefficient  matrix  has  a  regular  singular  point  at  infinity 
(r  =  -l),  or  the  case  when  the  coefficient  matrix  has  a  regular  point  at 
infinity  (r  <  -2).  In  our  transformations  to  canonical  form  in  Chapter  II 
of  the  thesis  we  have  given  special  consideration  to  each  of  the  above  cases; 
we  have  explicitly  constructed  transformations  which  transform  the  system  (4.1) 
to  canonical  form  in  the  case  when  h  =  -1,  and  in  the  case  when  h  >  0  and 
the  lead  coefficient  matrix  (A  )  has  distinct  eigenvalues.  For  the  very 
complex  general  case  we  have  merely  stated  the  results  in  [32].  We  have 
however  proved  modifications  of  the  usual  canonical  forms  for  the  general 
case  and  for  the  case  when  h  =  -1,  which  are  as  easy  to  achieve  as  the  usual 
forms.  In  addition,  these  modifications  are  an  aid  to  obtaining  error  bounds 
in  Chapter  III.  For  present  purposes  it  is  convenient  to  denote  the  general 
transformed  system  by 


s  - 1 


(4.4) 


dW  r 
—  =  z 


C(z)W  =  z 


C.  z'k  +  C  (z)  z"S  W  (s  =  1,2 
k  s'  ' 


•  •  • 


) 


k=o 


where  r  is  an  integer  and  C(z)  has  the  properties  described  for  A(z) 
(see  Section  1.5). 
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Following  transformation  to  canonical  form  we  have  explicitly 
constructed  formal  independent  series  solutions  for  the  transformed  system 
in  all  cases.  These  formal  series  do  not  in  general  converge;  in  the  cases 
when  they  do  converge  the  formal  series  solutions  are  actual  solutions  of  the 
transformed  system.  In  particular,  if  in  (4.4)  r  <  -1  and  the  coefficient 
matrix  C(z)  has  a  power  series  expansion  which  converges  in  some  neighbor¬ 
hood  of  infinity  then  the  formal  series  solutions  obtained  in  Chapter  II 
converge  in  this  same  neighborhood  of  infinity  and  are  actual  solutions  of 

(4.4). 


Although  the  formal  series  solutions  generally  diverge  they  are 

far  from  being  useless.  It  is  known  for  example  [8]  that  in  all  cases  the 
formal  solutions  of  (4.4)  are  asymptotic  expansions  of  actual  solutions  of 
this  equation  in  certain  sectors  of  infinity.  Moreover  a  partial  sum  of  a 
formal  solution  is  easy  to  compute.  In  Chapter  III  of  the  thesis  we  obtain 
error  bounds  for  partial  sum  approximations  of  formal  solutions  to  actual 
solutions,  in  all  cases.  We  have  again  first  given  special  consideration  to 
the  cases  r  =  -2,  r  =  -1  and  to  arbitrary  non-negative  r  when  the  eigen¬ 
values  of  Cq  are  distinct.  We  then  concluded  Chapter  III  by  finding  error 
bounds  for  the  general  system  (4.4)  in  canonical  form;  the  results  of  this 
section  thus  include  the  previous  special  cases  together  with  all  other  remaining 
cases.  The  bounds  are  in  all  cases  given  in  two  forms;  as  vector  bounds  and 
as  norm  bounds , 

For  purposes  of  easing  the  difficulties  in  the  analysis,  the  simplest 
case  of  a  regular  point  (r  =  -2)  at  infinity  is  considered  first.  The 
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case  of  a  regular  singular  point  (r  =  -l)  is  somewhat  more  complex; 
it  was  thus  necessary  to  further  restrict  the  paths  of  integration  in  the 
Volterra  integral  equations  defining  the  solutions.  In  the  case  of  an  irregular 
singular  point  (r  >  0)  when  the  lead  coefficient  matrix  has  distinct  eigen¬ 
values  two  cases  distinguished  themselves:  the  case  when  it  was  possible  to 
define  the  error  vector  by  a  single  Volterra  vector  integral  equation; 
and  the  case  when  we  required  two  simultaneous  vector  integral  equations  to 
define  the  error.  The  error  bounds  are  considerably  sharper  in  the  former 
case  than  in  the  latter.  We  have  therefore  given  a  geometric  criterion  to 
distinguish  these  two  cases. 

By  our  method  of  constructing  formal  solutions  there  is  a  one-to- 
one  correspondence  between  the  formal  solution  vectors  and  the  eigenvalues 
of  the  lead  coefficient  matrix.  Thus  if  we  consider  the  eigenvalues  of  Cq 
to  be  points  in  the  complex  plane  and  to  be  enclosed  there  by  a  strictly  convex 

■3f 

closed  polygon,  then  corresponding  to  an  eigenvalue  which  is  an  extreme  point 
of  the  polygon  it  is  possible  (but  not  necessary)  to  express  the  error  vector 
by  a  single  Volterra  integral  equation;  corresponding  to  an  eigenvalue  which 
is  an  interior  point  and  not  a  boundary  point  it  is  necessary  to  use  a 
simultaneous  pair  of  Volterra  vector  integral  equations.  Corresponding  to  an 
eigenvalue  which  is  on  the  edge  of  the  polygon  but  not  an  extreme  point  it 
may  or  may  not  be  possible  to  express  the  error  by  a  single  Volterra  vector 
integral  equation. 

* 

We  have  assumed  here  that  end-points  of  integration  in  the  Volterra  integral 

equations  be  taken  at  infinity. 
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In  all  of  the  above  special  cases  it  is  possible,  under  the  assumption 
that  the  coefficient  matrix  is  holomorphic  in  a  sufficiently  large  domain,  to 
choose  end-points  of  integration  in  the  integral  equations  defining  the  solutions 
at  infinity.  This  is  no  longer  always  possible  in  the  general  case,  even 
when  the  matrix  C(z)  in  (H.4)  is  holomorphic  for  all  |z|  sufficiently 
large.  Nevertheless,  although  it  may  be  necessary  to  choose  one  finite  end 
point  of  integration  it  is  always  possible,  even  under  the  assumption  that  we 
are  constructing  actual  solutions  for  which  the  formal  solutions  are  asymptotic 
expansions  of  actual  solutions  to  express  the  error  vector  by  either  a  single, 
or  at  most  a  simultaneous  pair  of  Volterra  integral  equations. 

The  error  bounds  are  in  all  cases  carefully  estimated.  The  problem 
yet  to  be  solved  is  to  find,  among  all  the  possible  paths  those  that  minimize 
the  bounds;  in  this  connection  a  method  for  finding  the  paths  which  minimize 
the  total  variations  is  given  in  [6S 35 ]• 

For  sake  of  illustration  we  have  completely  solved  the  n  =  2  case 
in  Appendix  A  of  the  thesis.  This  is  accomplished  by  first  explicitly 
illustrating  with  a  flow  chart  a  method  of  transforming  the  general  n  =  2 
case  to  canonical  form.  Formal  solutions  are  then  constructed  for  the  trans¬ 
formed  system,  and  error  bounds  are  obtained  for  formal  partial  sum  approxi¬ 
mations  to  actual  solutions  of  the  transformed  system.  Assuming  the  domain 
in  which  the  coefficient  matrix  is  holomorphic  to  be  sufficiently  large,  it 
is  in  this  case  always  possible  to  obtain  actual  solutions  for  which  the 
formal  solutions  are  asymptotic  expansions  by  use  of  a  single  Volterra  vector 
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Appendix  B  contains  the  proof  of  existence  and  uniqueness  of 
solutions  of  integral  equations  in  the  thesis. 

In  Appendix  C  of  the  thesis  we  solve  some  real  variable  integral 
inequalities  which  are  useful  for  obtaining  error  bounds. 

A  difficulty  can  arise  in  practical  applications  if  solutions  are 
desired  corresponding  to  eigenvalues  in  the  lead  coefficient  matrix  of  the 
general  system  which  are  equal,  or  nearly  equal.  The  main  reason  for  this 
difficulty  is  the  problem  of  transforming  a  matrix  to  Jordan  canonical  form; 
this  problem  cannot  be  satisfactorily  solved  on  a  computer  in  the  case  when 
two  or  more  of  the  eigenvalues  of  the  lead  coefficient  matrix  are  equal,  or 
nearly  equal,  due  to  the  fact  that  only  a  finite  number  of  decimal  places  can 
be  carried. 

Now  it  is  shown  in  [5]  (Vol.  1,  Memoire  Troisieme)  that  in  the  case 

when  the  coefficient  matrix  of  the  system  (4.1)  has  a  convergent  power  series 

expansion  about  t  =  00,  the  actual  solution  matrix  of  this  system  which 

reduces  to  the  unit  matrix  at  an  arbitrary  point  t  in  the  domain  in  which 

the  coefficient  matrix  is  holomorphic  is  an  entire  function  of  t  ^  as  well 

as  of  the  matrices  A^,  A^,  A^,  Thus  a  small  change  in  the  eigenvalues 

of  Aq  produces  only  a  correspondingly  small  change  in  the  actual  solution 

matrix.  What  is  the  corresponding  behaviour  of  a  formal  solution  matrix  when 

two  or  more  eigenvalues  of  Aq  approach  each  other?  Perhaps  by  choosing 

appropriate  functions  of  the  eigenvalues  as  coefficients,  it  is  possible,  as 

* 

in  the  case  of  certain  second  order  differential  equations  to  take  linear 

* 


This  has  been  illustrated  to  me  by  Professor  E.L,  Whitney. 
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combinations  of  the  formal  vector  solutions  which  correspond  to  the  eigenvalues 
such  that  a  continuous  change  in  the  eigenvalues  produces  a  continuous  change 
in  the  formal  vector  solutions,  A  quantitative  investigation  of  this  problem 
would  be  desirable;  the  solution  of  which  might  very  well  throw  light  upon 
the  connecting  problem  ([47]  pages  7,  38)  as  well  as  solve  the  difficulty 
arising  as  a  result  of  the  fact  that  it  is  not  always  possible  in  practice 
to  transform  a  given  matrix  to  its  Jordan  canonical  form. 
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APPENDIX  A 


THE  n  =  2  CASE 


In  this  appendix  we  shall  completely  solve  the  general  n  =  2  case 
At  the  outset  we  explicitly  illustrate  by  flow  chart  a  method  of  transforming 
the  general  singular  n  =  2  system  (equation  (A.l))  to  canonical  form.  By 
applying  the  results  of  Chapters  II  and  III  of  the  thesis  we  then  obtain 
formal  solutions  for  the  transformed  system  and  error  bounds  for  partial 
sum  approximations  of  formal  solutions  to  actual  solutions. 


Let  us  assume  that  we  are  given  the  particular  system  X' 
(h  an  integer)  or 


z^A( z )X 


r  / 

xn 

/  'i 
X12 

V, 

'au(z) 

( - 

N 

V~C\J 

(A.l) 

X21 

x' 

22 

n 

=  z 

a21(z> 
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where  for  each  arbitrary  positive  integer  s 


s-1 


(A.2)  a^z)  =  y  aijk  z'k  +  aiJs(z)  z 


00 
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'  L  aijk  z 

k=o 


(i.J  =  1.2) 


and  where  we  assume 
and  bounded  in  some 


that  the  functions  a^(z)  (i,j  = 
domain  oB  extending  to  infinity. 


1.2) 


are  holomorphic 


We  begin  by  transforming  the  system  (A.l)  to  canonical  form.  The 
flow  chart  on  the  following  pages  indicates  the  steps  involved  in  trans¬ 
forming  the  most  general  singular  n  =  2  case  to  canonical  form. 


The  chart  is  constructed  along  the  lines  of  the  proof  of  Turrittin' 
theorem  [32],  and  therefore  contains  all  the  transformations  which  Turrittin 
used  to  prove  his  theorem  for  arbitrary  n.  It  follows  by  inspection  that 
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FLOW  CHART  ILLUSTRATING  TRANSFORMATION  OF  n  =  2  SYSTEM 

TO  CANONICAL  FORM 
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Make  transforma¬ 
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the  process  eventually  terminates  in  all  cases. 

Box  #2  in  the  chart  contains  the  normalizing  transformation,  box  #8 
the  exponential  transformation,  box  #19  the  shearing  transformation,  box 
the  root  equalizing  transformation,  and  box  #27  the  zero  inducing  transforma¬ 
tion.  Box  #7  also  contains  a  modified  form  of  the  zero  inducing  transformation. 

Although  it  is  evident  by  studying  the  above  flow  chart  that  the 
desired  canonical  form  can  be  achieved  without  use  of  the  exponential  trans¬ 
formation,  we  have  included  this  transformation  for  two  reasons:  firstly, 
it  makes  the  above  procedure  easier  to  understand,  and  secondly,  we  wanted 
the  chart  to  resemble  the  proof  of  Turrittin's  [32]  theorem  as  much  as 
possible.  The  resulting  system  thus  takes  the  form 

(A. 3)  ^  =  zr  C(z)W  =  [Q’(z)  +  z"2  B(z) ]  W 

where 
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where  Q'(z)  in  (A. 3 )  is  in  square  brackets  in  (A.U),  B(z)  =  [b^.(z)]  and 
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If  r  <  -2  then  Q'(z)  =  0.  If  for  some  integer  p  in  0  <  p  f  r+l 

ji  0  then  £,  =  0.  Otherwise  £,  may  be  either  zero  or  one.  The 
functions  b^(z)  ( i, j  =  1,2)  have  the  expansion 
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where  for  each  positive  integer  s  the  elements  b  (z)  (i,j  =  1,2) 

ijs 


are  holomorphic  and  bounded  in  a  domain  j)$l  extending  to  infinity. 


According  to  the  analysis  of  Chapter  II  the  system  (A. 3)  has  a 
formal  solution  matrix  of  the  form 
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Suppose  that  there  is  a  smallest  integer  p  in  0  <  p  <  r+l 

such  that  q.  -  q_  /  0.  Then  the  U.  .. 
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=  0  if  i  ^  j  and  0  <  k  <  r+l-p  (i,j  =  1,2) 
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in  this  case  each  column  vector  can  be  computed  independently  of 


If  on  the  other  hand  q^  -  q  =0,  p  =  0,l,***,r+l,  then 
»  Oj  Uq  =  i)  is  computed  from 


tv 
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£  =  0  or  1. 


r+1 

E  {  )  is  to  be  replaced  by  zero  if  p  =  r+1. 

s=p+l 
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In  either  case,  whether  £  =  0  or  1,  with  U 


o 


I  equation  (A.10) 


uniquely  determines  each  U^,  k  >  1.  If  6  =  0  each  column  vector  can  be 


computed  independently  of  the  other.  If  6,  =  1  the  second  column  vector 
in  (k  =  0,1,°" •)  can  be  computed  without  computing  the  first.  On 

the  other  hand  if  6=1  and  the  first  column  vector  in  (k  =  0,1,2, •*•) 

is  desired  the  second  column  vector  must  also  be  computed,  and  the  computa¬ 
tions  fcr  the  second  column  vector  must  precede  the  computations  for  the 
first. 


We  shall  obtain  an  actual  vector  solution  corresponding  to  the 


first  column  vector  of  formal  solutions,  and  we  shall  obtain  a  bound  on  the 
difference  between  this  actual  solution  and  a  partial  sum  of  a  formal  solution. 
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is  the  smallest  integer  such  that  q  -  q^  /  0. 
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be  an  actual  solution  vector  which  satisfies  the  equation  (A. 3),  Then  the 
error  vector  T-j^2)  =  w|(z)  ’  cp^m(z)  satisfies  the  differential  equation 

_p  q,(z) 

CA*  x7 )  T'lm(z)  '  =  Z  B^Zblm(z)  '  Rlm(z)  e 


Let  us  choose  a  non-negative  number  a  such  that  if  q^  =  q  , 


p  =  0,1,...,  r  then  a  =  0  if  |  £  |  =  oo  and  a  =  |£|  if  |£|  < 


00, 


If 


for  some  p  e  (0,l,...,r)  q,  ^  q_  then  we  set  a  =  a  where  a  = 
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max  (|u|  q^Cu)  =  if  1^1  =  oo,  while  if  |  £  |  <oo  we  assume 

£  e  oO  ,  £  does  not  coincide  with  a  zero  of  q^( z)  and  |£j  >  a^.  Let 
00*  =  0®  -  £  a  "  s  where  £  is  the  closed  disc  with  center  at  the  origin 

cL  cl 

and  radius  a,  S  is  empty  or  {oo}  according  to  whether  or  not  Q!(z)  =  0, 
We  now  define  a  region  z,  £)  to  be  the  union  of  all  points  z  (J  £) 

such  that  there  is  a  path  tP  connecting  z  and  £  which  satisfies  the 


following  conditions; 


l)  Except  for 


if  £  €  erf?  -  J,  f  lies  entirely 


in  o&; 


2)  if3  consists  of  a  finite  number  of  Jordan  arcs  t  =  t(s)  such 

that  ^  is  continuous  and  non-vanishing; 

3)  exp  [q2l(-r)  -  qgl(t)]  is  bounded  for  all  points  t,  t  in 

the  order  £,  t,  t,  z  on  ft  .  If  q  =  q^  for  p  =  0,1,..., r  and 
either  q^  q^  or  i*  ^  0  we  also  require  that  r/t  remain  bounded; 

b)  Yp  (t^  is  bounded,  where  q  may  be  taken  to  be  zero  if 

£  =  0  but  otherwise  q  is  an  arbitrary  positive  number  less  than  1. 


For  purposes  of  obtaining  actual  solutions  for  which  the  formal 


solutions  are  asymptotic  expansions,  and  in  order  to  obtain  sharp  error 
bounds  it  is  again  desirable  (whenever  possible)  to  choose  |£j  =  oo. 


OS’J  '?  li  •  )  t  tv  •  J-  y  i  JuXo  r  V'iD'.,  ,7B  j 
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If  it  is  possible  to  choose  £  such  that  |£|  =  oo  ,  then  by  making  an 
appropriate  rotary  transformation  z  =  co|  with  to  a  constant  and  jto|  =  1 
it  is  possible  to  make  £  =  oo  an  end-point  of  integration.  Moreover 
if  the  case  when  coefficient  matrix  C(z)  is  holomorphic  for  all  z 
sufficiently  large,  the  formal  solutions  of  the  system  (A. 3)  will  then  be 
asymptotic  expansions  of  actual  solutions  ^jm.k(2)  (which  we  shall  obtain 
below)  in  sectors  |arg  z  -  — |  <  (^|  “  <0  (k  =  0,l,*“,r-p) 

where  p  is  the  smallest  non-negative  integer  in  0  <  p  <  r  such  that 
f  q^p,  and  e  is  an  arbitrary  positive  number  less  than  j 

in  this  last  case  the  corresponding  end-points  of  integration  are 
=  00  exp  [  k  3  (k  =  0,1, •  •  •  ,r-p). 


In  Appendix  B  it  is  shown  that  if  p$(z,£)  contains  a  point 
distinct  from  £,  then  0$(z  , £)  is  an  open  connected  domain. 


It  follows  that  if  rl^Tn(z)  satisfies  the  equation 
(A. 1.8)  Blm(z)  =  f  exP  [Q(z)-Q(t)3  jt"2  B(t)  -  Rlm(t)  e  1 

where  the  integral  traverses  a  path  as  described  above  then  q  (z) 


also  satisfies  (A.I7).  On  substituting  q^(z)  =  €^m(z)  e 
this  equation  may  be  written 


dt 


in  (A. 18) 


(A. 19)  €lm(z)  =  ^  exp  [ (Q(z)  -  Iq1(z))-(Q(t)-Iq1(t))]  -jt  2  B( t )e lm( t ) -Rlm( t )|dt 


Let  us  now  estimate  the  integral 


»  *9  iitv«U  Is?  ^inl  *if 3  91,-rfv 


***  V  adw  nO  .(TI.a)  'i*ilal3ue  oel* 
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(A. 20)  R*m( z )  =  -  P  exp  [(Q(z)-Iq1(z))-(Q(t)-Iq1(t))]  R^t)  dt 


If  for  some  integer  p  in  0  <  p  <  r+1  we  have  q^p  /  q^p 


it  is  again  convenient  to  define 


* 


m- 1 


(A. 21)  p21(0  =  R21m(t)-q2'1(0u21mt-m+  J  [bu  (t)»u<  +  b12jm.s(t)u21s]t 


50-2 


s=o 


where  Rpim(t)  Siven  in  (A. 14).  Then,  proceeding  as  in  Section  3.4 
we  find  that 


m- 1 


*  /  V 

RIlm  z> 


”n 7 


[b  (t)u11  +b.  (t)u  .  ]t  m  ^  dt 

11 ,m-s v  '  11s  12,m-s v  '  2lsJ 


s=o 


(A. 22) 


R  (z)  =  u  (z  -e 
21m v  '  2lmv 


qo-.(z)-qoi(0 


*  .  v  ,  -m  n21x“'  'l2r3y-«>‘  * 


«  )  +  Wz> 


rz  q2i(z)-q2i(t) 

+  /  e 

'H 


d_ 

dt 


^21(t)\ 

vq2l(t)y+ 


m-1 


[b21>m-s(t)uilS  + 


s=o 


where 


b  ( t ) u  ]  t'm'2'l  dt 

22,m-s v  '  21s  ! 


(A. 23) 


*  ,  .  ppi ( z)  q21(z)-q21(0  p21(0 


72lm^Z^  qoi ( z)  ’  e 


l21 


q|i(J) 


On  the  other  hand,  if  q^(z)  =  q^z),  then 


Note  that  p 2^(t)  =0  if  the  smallest  integer  p  in  0  <  p  <  r+1  such 
that  q2p  -  q  /  0  is  r+1. 


■ 


0)  a  5id.*/v 

Jr.rfl  bnii  *w 


» 


j  '  -J  0  (Ti  q  >1  911  15  1  tl  0  -  (•,  i  93oW 
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(A. 24) 


Rllm(zK  =  /  1 

R21m(z)'  ^  l°4 


1 

log  ^ 


s=o 


11 ,m-s 


21 ,m-s 


(0 

(t) 


b12,m-s^t^ 


Us 


b22>m-s(t)y  Sis 


u 


12s 


u 


22s 


t"m-2  dt 


v£  log  t' 


We  thus  define  non-negative  numbers  u_  (i  =  1,2)  by 

llm 


Ullm=  SUP 


T  €  jf' 


(U11m  +  £  UTOm  log  t)t  m  -  (u,,  +  £  u  -  log  £)£ 


-m 


llm  s  12m  &  '  v  llm  s  21m 


Si-m  _  j»Tl1“,m 


(A. 25) 


u 


21m 


sup 


,  ,  .  ,  -m  q21^‘q21^ 

(u21m+?U22ml0ST)T  e 


(«2lm+5[l°g  f(u1im+u2lralosO+u22mlog?])5 


-m 


«  P 


T”ni“m  _ 


-m 


where  may  be  taken  zero  if  £  is  zero,  but  otherwise  is  an 

arbitrary  (fixed)  positive  number  less  than  1. 


Similarly,  with  t^  and  t  points  on  IP  in  the  order  £,  t^ , 


z  we  define 


(A.26)  711m  =  sup  | 


-ni 


m-1 


>>: 


(t)[ui1  +  £  log  t  U10  ]  + 

llm  t  ^hp'm+l-TJi  /  ,  1  ll,m-sv  /L  11s  ^  °  12s 

^  s=o 


+  b12,m-S(C)  [u21s  +  «  lo8  <=  usoJ)l 
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21m 


sup 


4m  <*2 


tl>t2  eiP  t2 


7i-m  _  ^7i 


-m 


+ 


t-li  r  q21(t2)-q21(ti) 


m+l-7i 


d  f  P21^tl^  \  m+2 

dti  V  ’si^i5  )  1  + 


m-1 


s=o 


b21,m-s(tl>  +  *  log 


ll,m-sv  1 


(t,)l 


Ulls  +  S  108  *1  u22s 


b22,m-s^tl^  +  ^  log  (t^  b12,m-s^tl^l 


U21s  +  «  108  C1  u22s 


m-2  _d_  f  P21^tl^  V 
1  dtl  V  **21^1^  ' 


* 


where  it  is  of  course  presumed  that  the  terms  / ^  and  p^^t^q^Ct) 
are  missing  if  q^(t)  =  q^t).  Conversely  if  q^(t)  7^0  t^en  £  is 
to  be  taken  zero.  If  I-  =  0  then  7X  can  be  taken  to  be  zero;  other¬ 
wise  is  an  arbitrary  positive  number  less  than  1. 


We  also  define  a,  3,  7  and  5  by 


( 1-r) )a  =  suP|t  V  (t)|;  (l-q)3=  sup  \t~\  (t)| 

t  €p  U  t  €  f  12 


(A. 27)  (1-7)7  =  sup  |t. 

ti , t2  £  ft 


“7 


b21(tl>  +  1  108  (If)  bll(tl>l 


d-^S  =  sup  +  e  log  (”) 

ti»t2  e  <y  L 


(i2l(t2)-q2l(t1) 


q21(t2)-ci2i(t:i) 


where  t^  and  t^  are  taken  as  in  (A. 26)  and  7  may  be  taken  zero  if 
|  is  zero  but  otherwise  7  is  an  arbitrary  (fixed)  positive  number 


less  than  1. 


128  - 


By  the  analysis  of  Appendix  B  the  solution  to  equation  (A.19) 
exists  and  is  a  unique  vector  of  functions  holomorphic  in 


Theorem  A. 1  If  corresponding  to  the  formal  vector  solution  W^(z) 
of  equation  (A. 3)  we  can  define  a  region  oQ(z,£)  as  described  by  the 
conditions  on  page  123  then  the  equation  (A. 3)  possesses  an  actual  solution 
vector 


of  functions  holomorphic  in  ( z , £ ) ,  such  that 


(A. 29) 


<  exp 


l£21m(z)l 


u*  Y (tT'i"n,ur  y^tni'm‘h 


(X  3  U1  1  ( t 

:  v, 


for  all  z  in  o 8  (z»5).  The  solution  W  (z)  depends  on  £  and  an 
arbitrary  positive  integer  m.  The  function  q^(z)  is  defined  as  in 
(AA);  the  numbers  u^^  are  defined  by  either  (A#9)  or  (A.IO).  The 
numbers  u.Qm»  ^ilm  an<^  a»  3>  ? *  5  are  defined  by  (A. 25),  (A. 26)  and 
(A. 27)  respectively.  If  £  =  0,  both  q  and  q!  may  be  taken  zero; 
otherwise  q  and  q!  are  arbitrary  positive  numb er s  less  than  one. 


By  the  analysis  of  Section  1.2  we  can  use 

a- 6 


a  p 


exp 


\ 


7  & 


sinh  x  /  ^ 


X 


6  cosh  y  0 

m  : 

0  cosh  x 


6 -a- 


(A. JO) 


to  explicitly  evaluate  the  right  hand  side  of  (A. 29). 
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APPENDIX  B 


EXISTENCE  AND  UNIQUENESS  OF  SOLUTIONS  OF  INTEGRAL  EQUATIONS 


In  this  Appendix  we  establish  the  existence  of  solutions  of  integral 
equations  used  in  the  thesis.  Among  other  works  Erdely's  [3?]  development  of 
a  real  variable  theory  applicable  to  single  Volterra  vector  integral  equations 
is  probably  related  most  closely  to  the  analysis  of  this  appendix.  Generaliza¬ 
tions  such  as  the  contraction  mapping  principle  (see  e.g.  [9])  can  also  be  made 
applicable.  We  shall  however  construct  a  proof  along  the  lines  of  the  usual 
differential  or  integral  equations  proofs  which  at  the  same  time  applies 
directly  to  our  problem  and  conditions. 

As  in  [20]  we  first  establish  the  existence  of  a  domain,  we  then 
proceed  to  construct  a  sequence  of  vectors  which  converges  to  a  solution.  In- 
order  to  avoid  unnecessary  repetition  we  shall  establish  existence  and 
uniqueness  only  for  the  general  singular  system.  Two  different  cases  arise 
however s  (a)  the* case  of  a  single  Volterra  (vector)  integral  equation,  and 
(b)  the  case  of  a'  simultaneous  pair  of  Volterra  (vector)  integral  equations. 

We  shall  consider  each  of  these  cases  separately. 


B . 1  One  End-Point  of  Integration 

Lemma  B . 1  If  under  the  generalized  extreme  eigenvalue  conditions  of 

Section  3.4.2,  <90  (z,  £)  contains  a  point  z  ,  with  z  bounded  and  distinct 
from  £,  then  $9(z,  0  is.  an  open  connected  doma in. 

Proof.  <7(9  defined  on  page  85  is  certainly  a  domain.  By  assumption 
there  is  a  path  0J  in  $0  connecting  z  with  £  satisfying  the  four 
generalized  extreme  eigenvalue  conditions.  Similarly  there  is  such  a  path 
(P  connecting  an  arbitrary  point  z^  (zq  bounded)  on  with  £.  We 

can  clearly  choose  a  radius  p  >  0  sufficiently  small  such  that  when  we 
join  an  arbitrary  point  z^  in  the  interior  of  the  disc  S(p)  (S(p)  is 
entirely  in  c&)  of  radius  p  about  zq  with  Zq  by  a  straight  line  tf9,, 
then  conditions  l),  2)  and  4)  are  satisfies  on  ru  ip  It  remains  to  be 
shown  that  condition  3)  is  also  satisfied. 
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To  achieve,  this  we  define  matrices 


(B.l) 


Ki(jp )  =  T>tUejp  lexP^i(T)‘Q1(,:)-I4  (<lj(T)-qj(t))]|  (1-1,2,...,  if  i  /  j) 


K1(*V  =  TS?t  €  y  le*PEQi(T)-Q1(t1)-I  (qj(T  J-qjO: 

1  1  U  1  1 


))] 


of  non-negative  elements  where,  in  equation  (B.l),  t,  t,  t  ,  t  are  points 

l  l 

on  PUP,  in  the  order  z,T»t,z,T,t,£  and  where  the  limit 
1  iii° 

superiors  are  taken  for  each  individual  element  of  the  matrices. 


Since  the  element  of  the  matrices  Q^(t)  (i  =  1,2,..., i)  are 
holomorphic  functions  bounded  in  0& ,  each  of  the  matrices  consists  of 
non-negative  bounded  elements.  In  fact  the  elements  of  K^(  )  (i  =  1,2,. ..,i; 

i  /  j)  are  bounded  by  hypothesis  of  this  proof.  We  now  choose  r  and  t 
on  pup  ,  in  the  order  z  ,  r,  t,  £.  We  then  consider 

1  l 

(B.2)  =  T  n  f(,JpJexp[VT>-Qi<t)-i .  (qj(T)-qj(t))]  |. 

(i  =  1,2 . t-,  1  /  j) 

If  t  (in  B.2))  is  on  jp  so  is  t  ,  and  in  this  case,  by  hypothesis  the 
matrix  on  the  right  consists  of  bounded  elements. 


On  the  other  hand  if  both  t  and  t  are  nn  ,  the  matrix 
on  the  right  is  similarly  bounded  (K^ (  )  of  equation  (B.l)).  Suppose 

then,  that  r  is  on  ^  and  t  is  on  .  In  this  case 


t\l  |J  .....s.l-i)  |[((>:  !-(»>  P-(t>  Pjiv.)  »  ( SO.,* 


Tedw  bn» 


.« 3  b»bwfod  avlSageo-oon 


. ;  rq  r>  r  »  .5  eJ  >;  Jot;^  r.6  babciuod  &im  (t  \  t 


■ 
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(B.3) 


exp  (Q.(t)-Q.(i:)  - 


(qj(r)-qj(t))] 


< 


|exp[Qi(zo)-Qi(t)-I 


(qj(z0)-qj(t))]|  lexp[Qi(r)-Qi(zo)-I^  (qj(r)-qj(zo))]| 


so  that  (SWJ  <  K  («>)  k'(£),  which  shows  that  the  elements  of 
K.(SV©  (i  =  1.2,  •••  ,Z‘,  i  /  j)  are  bounded,  and  simultaneously  establishes 
that  O$(z>0  is  an  open  connected  domain. 


We  next  establish  the  existence  and  uniqueness  of  the  solution  of 
the  integral  equation  (3,4.22),  To  achieve  this  goal,  we  first  observe  that 
all  the  known  functions  (i.e.  all  functions  except  the  solution  vector)  in 
(3.4.22)  are  holomorphic  and  single-valued  in  every  simply  connected  subset 
of  o$(z,f;).  If  o$(z,£)  jLs  multiply  -connected,  we  interpret  q6(z,0  (or 
any  open  subset  of  0$(z,£))  to  b<e  a  Riemann  surface  on  which  all  known  functions 
of  equation  (3*4.  22)  are  holomorphic  and  single-valued. 


We  need  only  establish  existence  and  uniqueness  of  the  solution 
of  (3.4.22)  for  each  t  in  an  arbitrary  open  bounded  subset  O#*(z,0  Of 
(z,£)  for,  from  this,  existence  and  uniqueness  follows  for  all  of  p$(z,£). 


We 


assume  £  is  in  0$*(z,C)  if  £  is  bounded.  If  £  is  unbounded  we 


choose  an  arbitrary  subfamily  {  (J  )  of  the  family  {(^)  (see  Section  3.4.4) 
such  that  if  is  in  {  $*)  then  =$  (Jtf  \  here  (f  is  a  fixed  path 

from  £  to  a  boundary  point  of  o$*(z,£),  while  lies  in  p$*(z,£) 

and  joins  with  an  arbitrary  point  t  of  $  (z,£).  With  rj  defined  as 


•* 


♦  *o5  (23  )  lo 
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in  the  definition  of  the  extreme  eigenvalue  case  (Section  3*4.2)  we  assume 
in  any  case  —  i.e.  whether  £  bounded  or  unbounded  —  that  ^ )  is 
to  remain  uniformly  bounded.^ 


As  is  usual  procedure,  we  define  a  sequence  {h^(r))  (k  =  0,1,***) 
of  n  x  1  vectors  of  holomorphic  functions  by  (see  equation  (3.4.22)) 


(B.4) 


ho(T)  =  0 


\+1(t 


)-r 

Ji 


D.(t)-D  (t) 
e  J  J  [t: 


B(t)hk(t)  -  JR;m(t)  .0%  ] 


dt 


s=v 


where  t  is  any  point  on  (p  ,  The  integral  is  taken  over  (P  where  ^P 
satisfies  the  generalized  extreme  eigenvalue  conditions  in  Section  3*4.2, 
together  with  the  remarks  above. 


Let  us  next  show  that  the  sequence  (h^(r))  is  well  defined,  and 
that  each  element  of  this  sequence  of  vectors  is  bounded.  By  equation 
(B.4)  we  have 


i.(t)  =  -  r  e 

1  Ji 


f  rV  (t)  dt 

/_,  jm'  (s  -  v); 


(B.5) 


s=v 


h,  i(T)  "  h.  (t)  -  ^  e  ^ 
k+1'  '  k'  '  ,/j 


t  D,(r)-D4(t)  _2 


t  B(t)[hk(t)-hk_1(t)]dt  (k  =  1,2, 


). 


+ 


In  our  proof  of  existence  and  uniqueness  all 
"independent  of  the  path",  and  hence  there 
the  restriction  from  {ft7}  such  that 

V(t 


such  that 


)  is  uniformly  bounded. 


our  contour  integrals  will  be 
:here  is1no  loss  of  generality  in 

(t1*  )  is  bounded  to  {  }  C[jp) 


. 


■ 


»  -3  i*v* 
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Clearly  h^(r)  is  a  well  defined  vector  of  holomorphic  functions  of 
r  ,  for  each  t  in  aQ*( z,£).  Further 


(B.6) 


VT>I  5  + 


where  the  vector  on  the  right  of  (B.6)  is  defined  as  in  Theorem  5.4.1. 
Suppose  that  we  have  verified  existence  and  holomorphy  for  h^(r), 


k  >  1,  together  with  the  inequalities 


q-lv.k-1 


{b yr  (f-1)) 

(b-7)  lhk^T^'hk-i^T)l  =  (k  .  i): 


where  B  and  r\1  are  defined  as  in  Theorem  5*4.1.  Adding  the 
inequalities  we  get 

k 

(B.8)  |\(t)|  ^  ^  lhs('«')-hs_1(T)| 

s=l 


S=1 


(s  -  i): 


jm  'p 


exp  {U;:^(tBi-”)  +  /^^(t^-1)} 


This  implies  by  (B.4)  that  h^^T)  also  weH  defined  as  a  vector  of 
holomorphic  functions  of  r  .  Further,  from  the  second  of  (B.5)» 
(5.4.25)  and  (B.7)  we  get 


*  r*,>$  >  ^  iw/i 


20I  y  fcp  70 Xr  f  I  »  &  3  8 1*  i>«i22isv  svorf  nr  3n/f3  >i  oq?»8 


'■  ;  i  9  9-id  ,f.  1-up  8  j.  :  <  'v 


.  J  I(t)x_8*M  •)/!  J  S  IW/I  (8.5) 


3:^8  W  8»1  JlJ  ■  >9fi 


c  1  iV  1  *•*  >4  t  >6  o>.  l  (T),  3arf*  (4.8)  aaliqj  x  alrfT 
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(B*9)  lhk+i(T)"hk(T) I 


<  JT |(H)M^l(^(fl)!1  (W1-)  + 


jm  T 


jm  **p 


f  T>  /  T^-lv^k-1 

<  r  _ _ L.vJ. _ _ _ drBy"  rtT1'1nfuv*YftT1l*m)4-7v  y^t-Bi-m-iv 

=  J,  (k  -  X):  dlBT?,t(t  )JlUjmV  +7jmV  ' 


} 


k! 


so  that  ( B . 7 )  holds  for  all  k.  The  same  is  true  for  (B.8). 


Thus  we  see  that  the  sequence  {h^(r)}  is  well  defined  for  r 
in  <fri*  ,£)•  Further  (B.7)  implies  that  (see  e.g.  [18]  pages  I66-I67) 


(B.10) 


W*  (t)  -  ' 11111  \ (t) 

jm  =  k  — ►  00  k  ' 


exists  uniformly  with  respect  to  r  in  p$*(z, £).  Thus  W^^t)  is 
holomorphic  in  this  domain.  Moreover  from  (B.4)  and  the  uniform  conver¬ 
gence  it  follows  that  ^jm(z)  satisfies  (3.4.22) ,  and  consequently  it  is 
a  desired  solution  of  this  equation. 


It  remains  to  show  uniqueness.  Suppose  W(z)  ji  W^(z) 
satisfies  (3. 4.22),  where  the  elements  of  W(z)  are  uniformly  bounded 
on  a  path  P  described  above.  Then 

pt  D.(r)-D.(t)  0 

(B.ll)  W(t)  -  WV  (t)  =  /  e  J  3  B(t)t~2(W(t)-W^(t))dt 


and  hence,  setting  a(t)  =  SUp ^  |W(t)  -  (t)|,  iterating  (B.ll) 

t  €  Qr  jm 

and  using  (3.4.23)  we  get 


jm 


i{(w-^o^r  Jt  ♦  r^,v  ^ 


■ 

\  .  f"-xr95'Y*"fHi  —  ' _ iJ—  « 


•  r£j-i  X  Jguq  [°-  ]  .[.3  !J<  ;  u  J  e;  JU<j(  i  (Y.f?)  rs  f3i;/r 


/  \ 

1  ^  OP  ’  •]  '  '  M  W 


03  339q^31  ff3Xv  ^XatoUfW  83h1X9 


-  J;>\  3  loo  -tnu  9d.  U  (u,v{)  «•«'•  «>VOI  3).  .ft!  .  fc/.'fj  Oi  >2rfC  tC  / 


al  ii  ^I3n*upr  r?oo  fan-  1.  .  )  tol  h  L  Jfte  (.*■  ,  .  M  *:rf3  awoIIo5  32 


oi38L-  a  i.U  io  (,'oJcJuIoa  btj/  ah  8 


(  I  <  .  o  •  L-  a  r  ■ .  •  ■  i  3*. 


-  (t)„>-{t)W 


JO),  -  (i)’-l  c/5  -  ■  (t)&  *«***»«  t»OC9rf  ()0B 
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{err  (tT1-1)}k 
_ 

kl 


(B.  12) 


|W(r)  - 


Wjm(T)l  = 


A(t) 


where  k  is  an  arbitrary  positive  integer.  However,  the  right  hand  side 
of  (B.12)  tends  to  the  zero  vector  as  k  -►  00,  so  that  the  assumption 
W(z)  ^  W^(z)  is  untenable.  This  proves  that  the  solution  (B.lo)  of 
equation  (3.4.22)  is  unique  in  the  space  of  n  x  1  vectors  of  functions 
holomorphic  in  0$  { z,  £)  which  remain  bounded  as  z  £  along  any  path 
P  specified  above. 


Theorem  B , 1  The  solution  to  equation  (3.4.22)  exists  and  is  a 
unique  vector  of  functions  holomorphic  on  the  Riemann  surface  qS  ( z ,  £  ) 
defined  by  the  extreme  eigenvalue  conditions  in  Section  Moreover , 

the  inequality  (B.8)  constitutes  a  bound  for  this  solut ion. 


B .2  Two  End-Points  of  Integration 

Lemma  B . 2  If  under  the  generalized  interior  eigenvalue  con- 
ditions  of  Section  3-4.2  (page  86),  c&  (^,  "H )  contains  z  with 

z  bounded  and  distinct  from  C,  and  £  ,  then  rl)  is 

12  1  2 

an  open  connected  domain. 

w 

Proof.  By  assumption  there  is  a  path  (P  in  c&"  is  defined  on 

page  85)  connecting  and  ^  ,  passing  through  z  and  satisfying  the 

four  generalized  interior  eigenvalue,  conditions  of  Section  3*4.2.  The  same 
path  P  connects  any  arbitrary  point  z  in  this  manner.  By  an  argu¬ 

ment  similar  to  that  used  in  the  proof  of  Lemma  B.l  we  can  clearly  choose 
a  radius  p  >  0  sufficiently  small  such  that  when  we  alter  P  to  pass 


*Ul~ri)  *Ysi 


■ 


>bi  •  bi  r><  3/fji  i  Arfl  v  '-  .inasinl  :■  Jl«oq  '  ei  ri  atari* 


o  7  j  s  orf*  j*  (<5l.T)  to 

al 


1><  np^oodT 


.S  .<t.£ 


?■  ss2  (r  «g  «(  '  .. "StListsS  ~  ~  ii 

i  , '  _.2>  >«!/>.  •.  : 


■  ,4.i-  j->  /9*j  lo  'no)  -.:bi  <  .'li.vrtdjjfcd  7ol  a  ftSiJ  '09j  it  ^ 
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through  the  point  z^  instead  of  through  z^  as  indicated  in  Figure  3,  where 
z  is  an  arbitrary  point  interior  to  the  disc  S(p)  then  conditions  l), 

2)  and  b)  of  the  generalized  interior  eigenvalue  conditions  remain  satis¬ 
fied.  Condition  3)  can  however 


FIGURE  3 


D  (x)-D  (t) 

:  Jl  Jl  a(t)t  1 


also  be  made  satisfied  by 

taking  p  sufficiently  small. 

To  prove  this  we  need  merely 

show  that  as  p  -+  0  and  the 

altered  path  in  Figure  3 

tends  to  <P  , 

D.,(t)-D  (t) 

:  Jl  Jl  a(t)t  n|| 


t,T  €  $> 


sup  ||  e 
t,T  € 


where  t  and  r  are  points  located  on  the  paths  and  tP  in  the  order 

t,  t,  £  .  A  more  convenient  way  of  achieving  this  is  to  put 

D,p(T1)-D1P(t:1) 

(3 . 1J! )  A(t,t1,T,T1)  =  e  J  J  a^tx)ti  ^ 


VT,'V0 

i  J  J  a(t}t 


where  t  and  t  are  located  as  described  above,  and  t  =  t,  ~  r 

outside  S(p)  while  in  the  interior  of  S(p),  t^  and  are 

arbitrarily  located  on  .  Let  us  denote  the  norm  on  the  left  in  (B.13) 
by  a*  and  that  on  the  right  by  a.  Then  it  follows  that 


(3.15) 


a'  <  a  +  sup  ||A(t,t 

tx  ,tx  e  I?'  ;  t , t  e 


' 


•  3 .  1  f  .  T  1  •  . ;  :■  d  •  '  ■  • 

J  Sf  ru  >1  sf.dftJiiu*  -J  Otflfiv  g  '3  io  btsll  six  T9}^  9ri3  o  (*t  bn*  (S 


9iu;;il  fix 

v  55 


(i)r.a-(T)r.a 

■  fj&q  rfj  .ro  b  :J  -  >1  >  trrJ  >q  s»i&  T  brrfc  J  9i*rfw 


j  c  1  ij  ! )  :  •  "V  ■  o  r  rc-..i  A  t-  , 

=  •  bod  i  ?r>  9  si  ~  oiia  •  xMw 

>r  i  j  1  •  ■;  •  1  ■  ('  ;• 
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The  norm  on  the  right  of  (B.15)  is  defined  as  a  sum  of  absolute  values  of 
holomorphic  functions.  With  f(z)  holomorphic,  log  |f(z)j  is  harmonic; 
the  maximum  modulus  principle  (see  Whitney  [3])  applies  and  the  norm  on 
the  right  of  (B.15)  tends  to  zero  as  p  -*0.  A  similar  argument  applies  to 
3,  7  and  S.  Moreover  the  difference  -'Tjp(tri  ^)|  (with  rj 

chosen  as  in  Theorem  3.4.3)  can  t>e  made  as  small  as  we  please  by  taking  p 
sufficiently  small.  Lemma  B.2  follows. 

Before  proceeding  to  establish  existence  and  uniqueness  of  the 
solution  of  the  equation  (3. 4.28)  we  again  observe  that  if  (jQ 
is  simply  connected  then  all  known  functions  of  equation  (3.4.28)  are 
single-valued  on  (fi(£vZ2>T\)'  li  on  the  other  hand  1»C2»rl)  is  multiply- 
connected  then  we  interpret  (jQ  {£ (or  any  open  subset  of  ,5 2,ti)) 

to  be  a  Riemann  surface  on  which  all  known  functions  of  equation  (3.4.28)  are 
holomorphic  and  single-valued. 


We  again  need  only  establish  existence  and  uniqueness  of  the  solution 
of  (3.4.28)  for  each  t  in  an  arbitrary  open  bounded  subset  d 9  (£1}£0,q) 

1  c. 

of  o0<Sl  for  from  this,  existence  and  uniqueness  follows  for  all  of 

$  (£x.£a.n). 


Let  us  choose  an  arbitrary  but  fixed  path  (P  belonging  to  the  family 
of  paths  {!?}  which  satisfies  the  conditions  of  Section  3*4.2.  We  next 
pick  any  two  distinct  points  and  ^  located  on  }P  in  the  order 

£>2*  Further,  we  assume  that  and  are  bounded.  Assume 


moreover  that  with  a,  3,  7  and  5  defined  as  for  Condition  4)  (Interior 
eigenvalue  conditions  —  Section  3. 4.2)  we  have  ^{cw- 5+\/r(a-62+437}^( tP  ^)  =  l-€ 


’ 


.  f.i  .  ■'  •  : . l  j 


2l&  t  2l  .  .  ~  L°./J 


.  i  '  '  '  ‘■f'JV  ■-  'U. 


where  e  is  a  positive  number  less  than  1,  We  now  form  a  class  [y  }  C  {  y } 

of  paths  such  that  each  P  coincides  with  the  fixed  path  P  everywhere 

except  on  that  part  of  P  which  joins  and  l  here  we  allow  & 

to  vary  such  that  -^{04- 5-t-s/~(a-  5 )  ^  +  VpV  t ^  <  1  -  -g-e.  Thus  we  define 

oO*Uvlz. n)  -  Up  where  ^ ^  is  that  part  of  ^  which  is  left  after 

we  remove  the  portions  from  (P  to  £'  and  from  C‘  to  (P. 

1  bl  2  2 


We  again  define  a  sequence  of  n  x  1  vectors  of  holomorphic 


functions 


by 


(B.16) 


h,  -  (  t  ) 

(hk(T))  =  {kl  > 
\2(t) 


(k  =  0,1,2,.*.) 


ho(r)  =  ( 


-hol_(T)- 


0 


ho2(T)  0 


hk2(r) 


r 

'Js, 


0 


0 


p 

J  t 

->  o 


0 


0  e 


-2 


V(t)  6(t)  h 


k-1,1  (t\ 

R.  - (t) 

/  Jml 

1  dt  +  ( 

k-1,2  ^ 

\  *  ,  * 
R.  0(t) 
jm2  * 

J  e  J  (k  >  1). 


where  r  is  any  point  on  ^  ,  where  the  first  integral  is  over  a  path  P 
from  to  t  and  the  second  is  over  ^  *  from  r  to  £  . 


The  sequence 


h  1  ( r  ) 
fki  ^ 


(k  =  0,1,...) 


is  well  defined. 


•  . .  lnn.i  J.  .  cr{?cq  o 
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For,  by  equation  (B.16)  we  have 


(B.17) 


/•hll^T^  -qj(r) 

\2(t)  Rjm2(T) 


,Vi,i(t)'\i(t)n  L  0  e  Jl'  7  Jl  0 


Djl(r)-DJ1(t) 


s. 


hk+l,2^'hk2^  0  j ^  0  e 


\  „Dj2(T)'Dj2(t) 


a(t)  p(t) 
7(t)  6(t) 


h,  „(t)-h,  ,  „(t)  C 
k2  *  k-1 ,2V  ' 


dt  . 
2  ’ 


by  the  definition  of  the  vectors  Rmjp(T)  (k  =  1>2)j  h^(r)  is  clearly  a 
well-defined  vector  of  functions  holomorphic  in  ($  (£^»£p>rl).  Further,  if 

we  define  functions  >  )  (k  =  1,2)  by  equation  then 


(B.18) 


< 


ii»>12(T)ir 


Suppose  that  we  have  verified  existence  and  holomorphy  for  h^(r),  k  >  1, 
together  with  the  inequalities  (B.18)  and 


(B.19) 


IVi,i(t)'Yi(t 


llhk+l,2^T^'hk,2^T 


< 


(k  >  1) 


tlrfqiooiG/vxf  i  I  i  >nu  o  iqiosv  haallsb-J  . 


P  .  (e‘« 
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This  implies  by  the  second  of  (B.l6)  that  *'S  a^so  weH  defined  as 

a  vector  of  holomorphic  functions  of  r,  for  each  t  in 


Let  us  prove  the  inequality  (B.19)  for  k  =  1,  By  equations 
(B. 18) ,  (B.I7)  and  (3.4.31)  we  find  that 


(3.20)  (  )<(  1 

'  h2,2<T>-v2(T>u  ^VT*s2»K52(tV1)' 


so  that  (B.20)  is  true  for  k  =  1.  Assuming  (B.20)  true  for  some  fixed  positive 
integer  k,  we  substitute  this  inequality  into  the  second  of  (B.I7)  to  obtain 


(B.21) 


hk+2,l^"hk+l,l^T^L  fy  °  \  /a  P\  Jhk+l,l^’hk,l^ 


^k+2,2^T^"hk+l,2^T 


51  „V 

J  T 


7  5  Hhk+l>2^t^“hk,2^t^ 


0--n)dt 


£7+2 


r 


k+l 


0  a  3  ^ 

0  /  2  7  6 


T 


( l-n)dt 


£7  +  2 


y  y  Vy,<!ry' 

0  <^tn'ly  7  5  VSx-«2)' 


The  first  of  the  above  inequalities  is  obtained  directly  from  the  second 
of  (B.I7)  together  with  (3*4.31)}  the  second  of  (B.2l)  is  obtained  using 
(B.20).  Thus  (B.I9)  holds  for  each  integer  k  >  1.  Moreover  summing  (B.19) 
and  (B.18)  we  find  that  under  the  conditions  imposed  on  , 


. 


■ 
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o 


(S1’T\  X>T(tI1  ^  “\  /  -  .  \  /ml 

)+(  )e  f  (Wj't''-1) 

WT>'  *m2^2>  °  ^,?2(tn'1) 


/-hk,l^T4  /ml 

(B.22)  (  )< 


(£i_.52) 

*»2«l-£8> 


for  all  integers  k  >  2,  where  the  matrices  E  and  f(  )  are  defined  as 
in  Theorem  3.4.3*  Consequently  (see  e.g.  [18]  pages  166 -67) 


(B.25) 


,Wjml<T> 


( 


hkl(x). 


j  —  lim 

\  v  *  .  /  k— >oc  .  ./ 

^jm2  T  hte<T> 


exists  uniformly  with  respect  to  r  in  oQ  Thus  (t)  is 

holomorphic  in  this  domain.  Moreover  from  (B.16)  and  the  uniform  convergence 
it  follows  that  (t)  satisfies  (3.4.28),  and  consequently  it  is  a  desired 

v 

solution  of  this  equation. 


In  order  to  show  uniqueness  we  suppose  that  there  exist  vectors  W^(z) 

and  W~(z)  of  sqme  dimension  as  ,(z)  and  wV  ^(z)  respectively,  such  that 

2  jmlv  '  jm2 

the  elements  of  these  vectors  are  holomorphic  in  ykr  ,  bounded  on  JP  , 

and  such  that 

W1(Z\ 

(B.2U)  (  V( 


W  ( z)  W.  _(z) 
2V  7  jm2v  7 


also  satisfies  (3.4.28).  Then 


Wi(t)-W^i(t) 


■J 


r  0  eDJl(T)'DJl(t) 


0 


1 


rrJ\  Di2(T)’Di2(t:) 

0  /  0  e  J  J 

J  z 
s2 


Wi(t)-w^mi(t). 


W.(t)-W^  _(t)  t 
2  '  jm2  * 


dt 

2 


a(t)  3(t) 

V(t)  6(t) 


(B.25) 


|<  S™)”  '  1  »•;  )  <■•■> 


oo  ri  o linw  »H3  L/je  (di«H)  <^Oi  1  *vo»loHf  .nX&mob  »i/(J  nX  oXdu'iotooloiS 


.aroi  >  *  io  i1  1  ul  • 
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fll 


* 
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Setting 


(B.26) 


A(t)n  /rB"ps>  *llwi(T>-wjmi(T)ll. 


A2(t)  TS“Pp  *HW2(T)-Wjn.2(T 


and  iterating  (B.25)  we  obtain 


(B.27) 


) 


where  k  is  an  arbitrary  positive  integer.  However  the  right  hand  side  of 
(B.27)  tends  to  the  zero  vector  as  k  -»-oo}  so  that  the  assumption  (B.2U)  is 
absurd.  This  proves  that  the  solution  of  the  equation  (3.H.28)  is  unique  in 
the  class  of  n  x  1  vectors  of  functions  holomorphic  in  qQ  which 

remain  bounded  on  any  path  (p  specified  above. 


Theorem  B.2  The  solution  of  equation  (3.I+.28)  exists  and  is  a  unique 
vector  of  functions  holomorphic  on  the  Riemann  surface  defined  by  the  interior 
eigenvalue  condition  in  Section  3. 4.2.  Moreover,  the  inequality  (B.22) 
constitute  a  norm  bound  for  this  solution. 


It  is  noteworthy  that  once  a  holomorphic  solution  vector  is  known 
at  some  point  t,  this  same  solution  vector  can  be  uniquely  extended  to 
every  bounded  subset  of  ()l§  by  a  method  in  [19].  We  shall  not  prove  this 
statement  here,  since  in  this  thesis  we  have  not  exploited  this  possibility 
with  respect  to  obtaining  sharper  error  bounds  near  the  origin. 
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APPENDIX  C 


SOLUTION  OF  INTEGRAL  INEQUALITIES  AND  PROOFS  OF  THEOREMS  IN  SECTION  3.4.5 


In  this  appendix  we  state  and  prove  some  extensions  of  Gronwall's 
[I7]  (or  Bellman's  [2])  inequality;  we  then  use  these  extensions  to  carry  out 
the  proofs  of  Theorems  in  Sections  3.3.5  and  3.4.5*  For  extensions  of  this 
inequality  to  finite  difference  equations  and  to  (scalar)  Volterra  integral 
inequalities  for  which  the  kernel  may  be  more  general  see  the  papers  [ 43 » 44 ] . 


C.l  Extensions  to  Matrices 


Lemma  C ,  1  For  k  a  positive  integer t  let  cp,  G  and  H  be  k  x  1 

vectors  and  let  E  and  F  be  k  x  k  matrices  of  piecewise  continuous  non¬ 
negative  functions  of  (7  in  a  <  a  <  b  such  that  F(cr)  E(cr)  and  J^F ( s )E( s )d 
commute .  Then  if  for  each  <r  in  a  <  a  <  b 

(C.l)  cp(ar)  <  E(<t)  r  [F(s)  cp( s )  +  G(s)]ds  +  H(cr) 

^  a 

the  following  inequality  also  holds  in  a  <  cr  <  b: 

cr  p(j 

a 


Proof.  From  the  statement  of  the  Lemma  there  exists  a  k  x  1  vector 

^(a)  of  non-negative  piecewise  continuous  functions  such  that  we  can  replace 
the  inequality  (C.l)  by  the  equality,  i.e. 

(C.3)  <p(cr)  =  E(ct)  f  [F(s)  <]>(s)  +  G(s)]ds  +  H(cr)  -  tfr(er). 

a 


exp  (  /  F(t)  E(r)dr){G(s)  +  F(s)H(s))ds 


(C. 2)  cp(a)  <  H(or)  +  E(a)  / 


The  equation  (C.3)  is  however  a  Volterra  integral  equation  in  cp  with 


”>  i  n»  •  a  -or*  a  .  .•  bno  i»3«  <>••;  *rt;  n 

IQ  \  A3  r  J  HrtO..  8.'l9*X9  9..3rfJ  d«0  ll  9*  [  yJi  •  t  ([S1  a  '  ir-ttHdCi  To)  ‘  V.*  ] 

*  IoV  (- ■  !  ■  )  r.1  k  ■  ■  i  V  it  •  E>  9J  > n ••  ;  03  yjjl  O  >r: 

«?9<2oq  ortj  393  J[6*T9rr9s  9‘JOri  ad  ys<o  Ismo;!  3/fj  riol/fw  lot  8'->Jt32I eoprini 


(c)If  ♦-  a!  («)0  +  («)<5  (e)3l  (tj)3  >  (d)<p  (1.0) 


salsfi  ?/  s?Haeg<>."*  aai«ali2).  ■  . 


ob((  +  (»)J)}{Tj  (t );  (r)  }  qxs  '»  (t>)3  +  (x>)K  >  (u)q>  (s  0) 


h  t  B9*I  9tl3  lo  Jrr»ra  4R3vi  &si3  '.toi*!  ,  LoovJ 

>w  lot!  f.  .  jo  roi  •  U  .rri  ir  e  tw.9:vj2rj  9Vl  t*>  (‘dJ'V 


.9.,  ,^H(vop»  srfl  yc/  .0)  y  ;  01/;  rr fc  arf  ! 


;oK  -  (v)H'  4-  sr>[ (e)»  +  (a)*  (o)3)  1  (t>):  «  (-o)a:  (£.0) 


'  o i  yp9  Ic|;  (i  a*is:Jl©V  a  iov&woW  8i  ?.0)  nci^aopo  arf 
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unique  solution  (see  e.g. 


[16]). 


Decreasing  the  magnitude  of  the  elements  ^(cr)  increases  the  magnitude  of 
the  elements  on  the  right  hand  side  of  (C.4).  The  inequality  (C.2)  follows. 


Lemma  C.2  Let  E  ^(cr)  exist  for 
exist .  Then  along  with  (C.2)  the  ( less 
inequality 


<T  i£  [ a , b  ]  and  let  ^  !*■ ) 

sharp  but  easier  to  evaluate) 


(C.5)  ®(ct)  <  E(ct)  exp  {  f'  F(t)E(t)(1t) 


[E_1(a)H(a) 


Y 


a 


a,cr 


(e'!h) 


r 


a 


G(s)ds ] 


also  holds . 


Proof .  The  inequality  (C.5)  follows  directly  if  we  integrate  the  right 
of  (C.2)  by  parts. 

It  is  noteworthy  that  if  P  dE  X(s)  H(s)  >  0  over  each  arbitrary 

A  _  i  i  ^ 

sub-interval  A  on  [a,b],  then  (E  H.)  =  E  (<t)H(ct)  -  E  (a)H(a); 

cl  j  U 

in  this  case  the  inequality  (C.5)  is  further  simplified.  If  in  addition  E(s) 
is  a  scalar  then  (C.5)  reduces  to 


In  this  thesis  we  have  not  attempted  to  make  full  use  of  this  property, 


* 


Unique  in  the  class  of  k  x  1  vectors  of  piecewise  continuous  functions. 


' 


r  ?3  0  U  I  yi  ">ffl  1  oj  *  |«  :  '.f-  Jon  ti  ;■  8JC»ar i$  8i/f  J  III 
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particularly  with  respect  to  our  analysis  concerning  the  case  of  the  regular 
singular  point.  The  derivation  of  the  inequality  (C,6)  leads  us  to  conjecture 
that  it  is  possible  (by  making  the  transformation  eYm(z)  =  zU  q^m(z)  in  the 
integral  equation  (3.2.9)  and  appropriately  altering  the  conditions  on  the 
paths  of  integration)  to  obtain  error  bounds  for  the  regular  singular  (and 
the  general)  case  without  modifying  the  usual  canonical  forms. 

C . 2  Extensions  to  Fredholm  Integral  Equations 

Lemma  C ,  3  For  lc  a  positive  integer ,  let  G  (1c  x  k)  II  (k  x  k), 
q>  (k  x  l)  and  0  (k  x  l)  be  matrices  of  non-negative  piecewise  continuous 
functions  of  s  in  ( a ,  b  )  .  JJi  in  addition  we  are  given  that  the  inequality 

rh 

(C.7)  cp(cr)  <  G(cr)  H(s )  cp(  s)ds  +  0(a) 

^  a 

holds  for  each  a  in  (a,b),  and  that  the  eigenvalues  of  the  matrix 

P  b 

(C.8)  F  =  /  G(s)  H(s)  ds 

'  a 

are  less  than  one  in  magnitude,  then  the  inequality 

1  Pb 

(C.9)  cp(cr)  <  0(a)  +  G(a)  {i  -  F]  I  H(s)  0(s)ds 

'a 

also  holds . 


Proof .  From  the  statement  of  the  lemma  there  exists  a  unique  k  x  1 
vector  ^(a)  of  non-negative  piecewise  continuous  functions  sue’:  that  ue 
can  replace  the  inequality  (C.7)  by  the  equality 


rr'tl  oi  uox/  1  Jr;  s  jiwvDsJtq  svll  -Md  ugh  %o  (o)ty  soi oov 


(d  (Y  )  Jjl  «;■  r?l  o  » jiiJ  f  r 
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r  b 

(C.io)  Cp((j)  =  '  G(cr)  H(s)  cp(s)ds  +  0(a)  -  ^(c j) 

a 

The  equation  (C.IO)  is  however  a  Fredholm  integral  equation  in  cp;  if  no 
eigenvalue  of  the  matrix  F  (in  (C.8))  is  one,  then  (C.IO)  has  the  unique 
solution  (see  e.g.  [16]) 

(C.ll)  cp(cr )  =  0(a)  -  +  G(a)  { I  -  F }  1  f  H(s)  {©(s)  -  f  (s))ds 

a 

If  in  addition  all  the  eigenvalues  of  F  are  less  than  one  in  magnitude,  then 

00 

(C.12)  [I  -  F)"1  =  y  Fv 

v=o 


and  the  sum  of  the  power  series  on  the  right  of  (C.12)  is  a  matrix  of  non¬ 
negative  elements.  Thus  decreasing  the  magnitude  of  the  elements  ^(a) 

(and  Tjr(s))  on  the  right  of  (C.ll)  increases  the  elements  of  cp(cr).  The 
inequality  (C.9)  follows. 

C. 3  Proofs  of  Theorems  in  Sections  3.3.6  and  5.4.3 

For  sake  of  avoiding  unnecessary  repetition  we  shall  only  prove  the 
theorems  in  Section  3.4. 5j  since  the  Theorems  in  Section  5*3*5  are  special 
cases  of  those  in  Section  3*4.5 

* 


Unique  in  the  class  of  k  x  1  vectors  of  piecewise  continuous  functions. 


. 
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The  vector  functions  i)  ,(£_,z)  and  il/  defined  by 

ml'  1  '  7  m2  2 

(3.4.44)  or  (3.4.45)  depending  on  whether  £  is  unbounded  or  bounded  res- 

* 

pectively  are  integrals  of  vectors  of  absolute  values,  To  simplify  the 
analysis  of  this  section  we  write  these  vector  functions  in  the  form 


(c.13) 


^  cnl(t)  I dt I 

0(z,t  )  =  r  2  C  (t)  I dt I 
rm2v  ,b2'  /  m2x  '  1  1 

v  z 


Using  the  definitions  (3.4.43)  for  a,  (3,  7  and  5  together 
with  the  results  of  the  previous  section,  we  obtain  from  (3.4.34) 


(C. 14) 


< 

■  r  {[a^ml(t)|+P|wVm2(t)|]|^|  +  Cnl(OJ|dt| 

- 

J  2ndw;mdOM|w;m2(t)|]|^i  +  cra2(t))|dt|J 

Z  L 

By  Section  3*4.2,  the  contour  with  end  points  ^  and  £ 
is  mapped  onto  the  real  interval  [a,b],  the. points  t  =  z,  and  £ 

corresponding  to  the  points  s  =  a,  cr  and  b  respectively.  Hence  with 

u)(s)ds  =  ,  0k(s)ds  =  Cmk(t)  |  dt  j  ,  |Wjmk(z)|  =  cpk(cr) ,  (k  =  1,2), 

the  equations  (C.14)  may  be  written 


(C.15) 


ra 

CP1(0')  <  1  { [a  cp^  ( s )  +  pq>2(s)]  w(s)  +  e^s)}  ds 


pb 

q>p(<0  <  /  {[7  cp^s)  +  &cpp(s)]  w(s)  +  0  ( s ) }  ds 

-  cr  ^ 


# 

Variations  of  vectors  over  a  contour. 


1*331:1  t  dd  y«jzk  (;U.O)  uaoldfci/p©  od3 


—  ■■"  — ""  ■-  . - . — * 
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On  applying  Lemma  C.l  to  each  equation  in  (C„15)  we  obtain 


cp 


pCT  p(J 

(ct)  <  /  exp  {  /  a  w(T)dT]  {3  cp  (s)  u>(s)  +  0  (s)}  ds 

a  li  a  ^ 


(G. 16) 


pb  nS 

cpp(^)  5  /  exp  {  /  &  0)(T)dT}  {r  CP  (s)  co(s)  +  6  (s))  ds 

^  1  cr  cr  1  ^ 


Substituting  the  bound  for  cp^(s)  in  the  second  of  (C„l6)  for 


cp^s)  in  the  first  we  obtain 


p(J  pc  p  b  p  u 

(C.I7)  cp  (cr)  <  /  exp  {a  '  w(r)dT)[3  j  exp  [5  '  w(r)dT}{7  cp^u)  co(u)  + 

1  “  1 '  a  s  ' '  s  *•'  s 


+  @2(u)}  du  w(s)  +  0^(s)]ds 


The  repeated  integral  (C.I7)  may  be 
interpreted  as  a  double  integral  over  the 
shaded  region  indicated  in  Figure  4. 
Splitting  up  the  integral  as  an  integral 
over  a  triangle  plus  an  integral  over  a 
rectangle  and  interchanging  the  order  of 
integration  in  each  case  we  obtain 


(C.  18)  cpi(cr)  <  0  (o’) 


n cr  p  r  po  pT 

j  j  exp{a  '  w(u)du}3  exp{&  '  w(v)dv}{7cp  (t)u>(t)  +  0p( t) )w( s )ds  dT 

'1  3  >'  Cl  '  ■  o  o  t- 


a  a 


pb  pff  p  cr 

+  '  '  expfa  '  co(u)du 


r 


T 


cp{a  ’  w(u)du}3  exp{5  oj(v)dv}(7cp  (t)w(t)  +  0  (  t)  }<d(s  )ds  dT 
ct  ’  a  s  s  1  ^ 


■ 


■ 
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where  we  have  put 


(C.19) 


r a 

0  (cr)  =  '  exp 
'  a 


r  a 

{a  co(u)du]  0  (s)ds 


On  evaluating  the  inner  integrals  in  (C.lB)  we  find  that 


(C.20) 


x2(^>t)}{7  cP1(t)co(t)  +  0g(r )  ]dT 


where  X^(ct,t)  and  X^cTjt) 


are  given  by  the  solutions  of 


a  X1(ct,t)  +  X1(o-,t)5 

p  (J  p  T  p  T 

(C.21)  =  exp{a  /  w(s)ds}[(3  exp{5  /  to(s)ds}  -  exp  {-a  (s)ds)$] 


a  X2(ct,t)  +  X2(o-,t)6 


p  0"  p  0"  p  T 

=  [exp{a  /  td(s)ds}(3  -  3  exp{-  5  '  oo(  s  )ds }  ]exp{5  J  u)(s)ds} 


a 


According  to  Bellman  [11]  both  of  the  above  equations  have  unique 
solutions  provided  +  v ^  /  0  (i  =  1,2,  ••*,«;  j  =  1 ,2,  •  •  •  ,n-K ) ,  and 

Vj  being  the  eigen  values  of  a  and  5  respectively,  and  where  a  is  a 
k  x  k  matrix.  By  another  theorem  of  Bellman  [2],  if  this  condition  is 
not  satisfied,  it  can  always  be  made  satisfied  by  an  arbitrarily  small  change 
in  the  elements  of  the  matrices  a  and  5.  Nevertheless  an  exact  explicit 
expression  of  the  right  hand  side  of  (C.20)  is  not  practical  for  obtaining 
error  bounds,  and  we  make  a  slight  sacrifice  of  sharpness  of  error  bound  in 
order  to  achieve  a  matrix  bound  which  is  much  simpler  to  evaluate.  Several 
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simplifications  of  the  bounds  are  possible  at  this  point;  one  of  these  is 
as  follows: 


p  T  p  (J  p  T 

(C.22)  exp  (a  bo(u)du}(3  exp  {&  '  oo(v)dv}  u)(s)ds 


s 


J 


< 

‘  1  a 


r,T  p  ( J  r\  T 

/  exp  (a  w(u)du}3  exp  (5  '  w(  v)dv}  w(s)ds 


cl 


-1  r0”  r*  fT 

=  a  [exp  [a  '  co(u)du)  -  exp  [a  '  co(u)du}]3  exp  [5  I  w(v)dv} 

1  a  t  a 

-1  ra  rT 

<  a  [exp  [a  I  u)(u)du}  -  I  ]|3  exp  [5  '  u)(v)dv] 

=  Jo  K  J*. 


* 


Similarly 


(c.23) 


ra 

[a  j  co(u)du]  3  exp  [5 
^  s 


rT 

1  oj(v)dv]  u)(s)ds 

y'  s 


< 


-1 

a 


[exp  [a 


ra 

'  u)(u)du]  -  I  ]  3  exP 
J  a  K 


fB 


rT 

J  oo(v)dv] 
a 


=  G(o-)  H(t) 


where 


(C. 24) 


-1  r 

G(cr)  =  a  [exp  [a  '  u)(u)du]  -  I  ] 

a  K 

H(t)  =  exp  [6  f  co(v)dv] 

'■  a 


The  possible  difficulty  arising  as  a  result  of  the  vanishing  of  the  deter¬ 
minant  of  a  is  avoided  if  we  remember  that  the  exponentials  are  defined 
by  power  series  expansions. 


. 
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Thus  combining  the  above,  we  have 


(C.25) 


q>(cr)  <  0  (ct)  + 


r 


G(a)  H(t)  {7cp(t)  w(t)  +  02(t)}c1t 


On  applying  Lemma  C.3  to  equation  (C.25)  and  transform  back  to  complex 
variable  notation  we  obtain  the  results  of  Theorem  3.4.4  (and  (3.3.4). 


Let  us  now  proceed  to  prove  Theorem  3.4.5  (and  3.3.5).  For  this 
purpose  we  let  norms  a,  (3,  7  and  S  be  defined  as  in  Theorem  3*4.3*  If 
we  define  non-negative  scalar  continuous  functions  cp^,  cp^,  oj,  Q ^  and  0^, 

and  0 ^  hy  cPk((J)  =  llWjmk(z)ll»  cP](s)cls  =  cP2^S)dS  =  I  d  Vyp  ( t  >  C  g )  | 

(equation  (3-4.42))  w(s)ds  =  =  1-^-i-dIlJ —  — 4(fJ- 1  c{s .  and  6  is 

t2^  t(s)2"Tl  05  3 

defined  using  the  definitions  of  this  paragraph  together  with  (C.19).  Again 
a  s  <  b . 


Starting  with  equations  (C.15)  where,  instead  of  having  vector 
inequalities  we  now  have  scalar  inequalities,  the  analysis  proceeds  exactly 
as  above  to  equation  (C,2l).  Since  however,  the  elements  a,  (3,  7  and  5 
along  with  X^(ct,t)  and  X0(cr,r)  are  now  scalar  quantities,  we  are  able 
to  explicitly  solve  the  equations  (C.2l),  to  obtain 

rb 

(C.26)  cp(cr)  "S  ^^(Cr)  +  X(t,Ct)  [7CP(t)  w(t)  +  (  T )  )  dT 

^  a 


where 


(C.27)  X(T,cr)  = 


p  p  cr  p  t  p  cr 

“  ^  [exp  [a  ■  co(  s  )  ds  +  S  w(s)ds)  -  exp  [a  w( s  )ds  }  ]  ,a<T<cr 
a  a  t  — 


0W-& 


p  cr  p  t  p  t  r 

[exp  [a;  u)(s)ds  +  5  '  U)(s)ds]  -  cxp[&  '  oo(  s )ds }  ]  ,o<T<b  . 
a  a  cr 


i  n:  ;i  bn  ( •)  i  q:  >  •  ^aly;I  a  nO 


j  .  i-r 
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We  over-estimate  the  bound  on  the  elements  X(t,<j)  in  order  to 
obtain  an  integral  equation  which  satisfies  Lemma  C.3.  Over  the  whole  range 
of  (a,b)  each  of  the  functions  in  (C.27)  satisfies 

(C.28)  X(t,ct)  <  X*(cr)  =  exp  |a  P  w(s)ds  + 

a  +  J  a 

+  |  |  P  oo(s)dsT  sinh  P  u>(s)dsl-  . 

d  '■  a  1  d  -'a 

where  0,  =  max  (a, 5). 

On  substituting  the  right  of  (C.28)  into  (C.26)  we  again  obtain  an 
integral  inequality  of  the  form  (C.25);  applying  Lemma  C.3  to  this  inequality, 
the  results  of  Theorem  3.4.5  follow. 


I 

’ 


